AFRL-RV-PS- 

TR-2015-0116 


AFRL-RV-PS- 

TR-2015-0116 


RADIO  INTERFERENCE  MODELING  AND 
PREDICTION  FOR  SATELLITE  OPERATION 
APPLICATIONS 


Tien  M.  Nguyen,  et  al. 


Department  of  Electrical  Engineering  and  Computer  Science 
The  Catholic  University  of  America 
Washington,  DC  20064 


25  Aug  2015 


Final  Report 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


AIR  FORCE  RESEARCH  LABORATORY 

Space  Vehicles  Directorate 

3550  Aberdeen  Ave  SE 

AIR  FORCE  MATERIEL  COMMAND 

KIRTLAND  AIR  FORCE  BASE,  NM  87117-5776 


DTIC  COPY 

NOTICE  AND  SIGNATURE  PAGE 


Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for 
any  purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U.S. 
Government.  The  fact  that  the  Government  formulated  or  supplied  the  drawings, 
specifications,  or  other  data  does  not  license  the  holder  or  any  other  person  or  corporation; 
or  convey  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that 
may  relate  to  them. 


This  report  is  the  result  of  contracted  fundamental  research  deemed  exempt  from  public  affairs 
security  and  policy  review  in  accordance  with  SAF/AQR  memorandum  dated  10  Dec  08  and 
AFRL/CA  policy  clarification  memorandum  dated  16  Jan  09.  This  report  is  available  to  the  general 
public,  including  foreign  nationals.  Copies  may  be  obtained  from  the  Defense  Technical 
Information  Center  (DTIC)  (http://www.dtic.mil). 


AFRL-RV-PS-TR-20 15-01 16  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR 
PUBLICATION  IN  ACCORDANCE  WITH  ASSIGNED  DISTRIBUTION  STATEMENT. 


//SIGNED// 

PAUL  HAUSGEN,  Ph.D. 

Technical  Advisor,  Space  Based  Advanced  Sensing 
and  Protection 


//SIGNED// 

JOHN  BEAUCHEMIN 

Chief  Engineer,  Spacecraft  Technology  Division 
Space  Vehicles  Directorate 


//SIGNED// 
STEVEN  A.  LANE 
Program  Manager 


This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange,  and  its 
publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its  ideas  or  findings. 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

25-08-2015  Final  Report  23  Sep  2014  to  24  Aug  2015 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 

FA9  453-14-1-0259 _ 

Radio  Interference  Modeling  and  Prediction  for  Satellite  5b. GRANT  NUMBER 

Operation  Applications 

5c.  PROGRAM  ELEMENT  NUMBER 

62  60  IF 


5d.  PROJECT  NUMBER 

8809 


5e.  TASK  NUMBER 

PPM00019115 


5f.  WORK  UNIT  NUMBER 

EF122748 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


6.  AUTHOR(S) 

Tien  M.  Nguyen,  Hien  Tran,  Zhonghai  Wang,  Amanda  Coons,  ZiJian 
Mo,  Genshe  Chen,  Charles  C.  Nguyen,  Gang  Wang 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Electrical  Engineering  and  Computer 
Science  The  Catholic  University  of  America 
Washington,  DC  20064 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 
Space  Vehicles  Directorate 
3550  Aberdeen  Ave  SE 
Kirtland  AFB,  NM  87117-5776 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL/RVSV 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AFRL-RV-PS-TR-2 01 5-0116 


14.  ABSTRACT 

This  report  describes  an  innovative  radio  frequency  interference  evaluation  and  prediction 
framework.  This  report  presents  associated  modeling  and  prediction  tools  to:  (1)  detect  the 
presence  of  radio  frequency  interference  using  receiver's  synchronizer,  (2)  estimate  the 
radio  frequency  interference  characteristics,  (3)  evaluate  effectiveness  of  the  existing 
unified  S-band  command  waveforms  employed  by  civil,  commercial  and  military  satellite 
operations  ground  stations,  and  (4)  predict  the  impacts  of  radio  frequency  interference  on 
unified  S-band  command  systems.  The  proposed  framework  and  radio  frequency  interference 
prediction  tools  allow  a  communication  designer  to  estimate  the  optimum  transmitted  signal 
power  to  maintain  a  required  satellite  operations  quality-of-service  in  the  presence  of  both 
friendly  and  unfriendly  radio  frequency  interference  sources. 


15.  SUBJECT  TERMS 

space  communication;  satellite  communication,  satellite  operations,  radio 
frequency  interference 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Steven  A.  Lane 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

Unlimited 

138 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


(This  page  intentionally  left  blank) 


Approved  for  public  release;  distribution  is  unlimited. 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES . iii 

LIST  OF  TABLES . vi 

1  SUMMARY . 1 

2  INTRODUCTION . 1 

2. 1  Description  of  Problems  for  RFI  Detection  and  Prediction . 3 

2.2  Impacts  of  RFI  on  Synchronization  Loops  Performance . 6 

2.3  Scope  of  Effort . 8 

2.4  Organization  of  Report . 9 

3  METHODS,  ASSUMPTIONS,  AND  PROCEDURES . 10 

3.1  Assessment  of  Existing  RFI  Tools  and  RFI  Tool  Selection . 10 

3.2  Framework  for  RFI  Tool  Development  and  RFI  Assessment . 12 

3.2.1  Framework  for  RFI  Detection  and  Prediction  Tool  Development . 12 

3.2.2  Framework  for  RFI  Detection  and  Prediction . 13 

3.2.3  Framework  for  RFI  Assessment . 15 

3.3  RFI  Analytical  and  Simulation  Model  Development . 16 

3.3.1  Analytical  Models  for  Unified  S-band  Wavefonns . 16 

3.3.2  Analytical  Models  for  RFI  Signals . 23 

3.3.3  Simulation  Models  for  the  Second-Order  Phase  Locked  Loop . 24 

3.3.4  Existing  Carrier  Acquisition  Approaches . 29 

3.3.5  Analytical  Models  for  Assessing  PLL  Acquisition  Perfonnance  in  the  Presence  of 

RFI  . 37 

3.3.6  Analytical  Models  for  Assessing  PLL  Lock  Detector . 42 

3.3.7  Analytical  Models  for  Assessing  PLL  Tracking  Performance . 45 

3.3.8  RFI  Detection  Models . 50 

3.3.9  USB  Signal  Estimation  Models . 53 

3.3. 10  RFI  Signal  Estimation  and  Prediction  Models . 66 

3.3. 1 1  USB  Link  Budget  Models  for  Assessing  SATOPS  Performance . 95 

4  RESULTS  AND  DISCUSSION . 100 

4.1  RFI  Detection  and  Prediction  Model  Integration . 100 

4.2  IFT  SATCOM  Tool  Integration . 101 

4.3  Verification  and  Validation  of  RFI  Detection  and  Prediction  Tools . 106 

5  CONCLUSIONS . 114 

Approved  for  public  release;  distribution  is  unlimited. 


1 


6  RECOMMENDATIONS . 116 

REFERENCES . 118 

LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS . 119 


Approved  for  public  release;  distribution  is  unlimited. 

ii 


LIST  OF  FIGURES 


Figure  1.  A  Typical  Satellite  Control  Network . 2 

Figure  2.  Typical  Friendly  RFI  Scenarios . 4 

Figure  3.  Definition  of  Interference  Bandwidth  and  Corresponding  CCI  and  ACI  Scenarios . 5 

Figure  4.  Notional  Subcarrier  Tracking  Jitter  Performance  of  PCM/PSM/PM  without  RFI . 7 

Figure  5.  Notional  Subcarrier  Tracking  Jitter  Performance  of  PCM/PSM/PM  with  RFI . 7 

Figure  6.  RFI  Detection-Prediction  Tool  Development  Framework . 13 

Figure  7.  RFI  Detection-Prediction  Models  Development  Framework . 14 

Figure  8.  BER  Performance  Model  for  Predicting  the  Impacts  of  RFI  on  SATOPS  Perfonnance 

. 15 

Figure  9.  RFI  Assessment  Framework . 16 

Figure  10.  PCM/PSK/PM  Modulator-Transmitter  Model . 17 

Figure  11.  PCM/Bi-Phase  Modulator-Transmitter  Model . 17 

Figure  12.  NRZ  Data  Format  for  PCM/PSK/PM-Sinewave  Subcarrier . 19 

Figure  13.  Power  Spectral  Density  for  PCM/PSK/PM-Sinewave  Subcarrier  Waveform . 19 

Figure  14.  Receiver-Demodulator  for  PCM/PSK/PM-Sinewave  Subcarrier  Waveform . 20 

Figure  15.  Bi-Phase  Data  Fonnat  Used  By  PCM/PM-Bi-Phase  Waveform . 21 

Figure  16.  Power  Spectral  Density  of  PCM/PM/Bi-Phase  Waveform . 21 

Figure  17.  Receiver-Demodulator  for  PCM/PM/Bi-Phase  Wavefonn . 22 

Figure  18.  Advanced  Second-Order  PLL  Employed  by  NASA  Spacecraft  and  Satellites . 26 

Figure  19.  Recursive  Implementation  of  F(z),  B(z),  K(z)  and  V(z)  Using  Bilinear  Transfonnation 

. 26 

Figure  20.  Typical  Second-Order  PLL  Employed  by  Civil  and  Commercial  Satellites . 29 

Figure  21.  Carrier  Acquisition  Approaches  for  USB  SATOPS  Systems . 30 

Figure  22.  Carrier  Frequency  Acquisition  Using  Frequency  Sweep  Technique . 32 

Figure  23.  Minimum  Sweep  Rate  and  Probability  of  Lock  as  a  Function  of  Receive  Power . 33 

Figure  24.  Phase  Acquisition  Process  of  a  PLL . 35 

Figure  25.  Typical  Maximum  Likelihood  (ML)  Phase  Estimator . 36 

Figure  26.  Typical  PLL  Lock  Detector  Employed  By  Civil,  Commercial  and  Military  SATOPS 
Satellites . 43 

Approved  for  public  release;  distribution  is  unlimited. 

iii 


Figure  27.  Proposed  Framework  and  Algorithm  for  RFI  Detection . 53 

Figure  28.  Linearized  Model  of  Advanced  NASA  PLL  Shown  in  Figure  18  for  Tracking  Mode56 
Figure  29.  Mathematical  Framework  and  Algorithm  to  Incorporate  PLL  Tracking  Error  into  USB 

Signal  Estimation  Model . 56 

Figure  30.  Simulation  Model  for  Incorporating  Carrier  Tracking  Error  into  NASA  Advanced 

Second-Order  PLL . 57 

Figure  31.  Simulation  Model  for  Incorporating  Carrier  Tracking  Error  into  Typical  Type-I 

Imperfect  Integrator  Second-Order  PLL . 60 

Figure  32.  Tracking  Phase  Error  as  a  Function  of  Time . 64 

Figure  33.  Difference  of  Tracking  Error  as  a  function  of  Time . 64 

Figure  34.  Phase  Plane  Plot  for  a  Typical  PLL  Using  Type-I  Imperfect  Integrator  for  ISR  =  -15 

dB . 65 

Figure  35.  Phase  Plane  Plot  for  a  Typical  PLL  Using  Type-I  Imperfect  Integrator  for  ISR  =  25 

dB . 65 

Figure  36.  Block  Diagram  of  the  USB  Carrier  Signal  Component  Estimator . 66 

Figure  37.  Block  Diagram  for  Extracting  the  CW  RFI  Signal  from  the  Received  Signal,  R(t)  ...  68 

Figure  38.  Block  Diagram  for  CW  RFI  Estimation,  Prediction  and  Detection  Model . 69 

Figure  39.  Block  Diagram  for  CW  RFI  Signal  Power  Estimation . 70 

Figure  40.  Block  Diagram  for  CW  RFI  Frequency  Estimation . 75 

Figure  41.  Block  Diagram  for  CW  RFI  Phase  Estimation  and  RFI  Detection . 77 

Figure  42.  Receiver  Operating  Curve  for  the  CW  RFI  Signal  Detector  Shown  in  Figure  41 . 78 

Figure  43.  Block  Diagram  for  Extracting  the  WB  RFI  Signal  from  the  Received  Signal,  RwbO)  79 

Figure  44.  Block  Diagram  for  the  WB  RFI  Estimation,  Prediction  and  Detection  Model . 81 

Figure  45.  Block  Diagram  for  the  WB  RFI  Power  Estimation . 82 

Figure  46.  Block  Diagram  for  the  WB  RFI  Signal  Frequency  Estimator . 85 

Figure  47.  Block  Diagram  for  WB  RFI  Phase  Estimation  and  RFI  Detection . 88 

Figure  48.  Block  Diagram  for  WB  RFI  Data  Stream  and  Rata  Rate  Estimation . 90 

Figure  49.  Block  Diagram  for  WB  RFI  Synchronizer . 92 

Figure  50.  Block  Diagram  for  WB  RFI  PSD  and  Data . 93 

Figure  51.  RFI  Modeling  and  Prediction  Tool  Integration . 102 

Approved  for  public  release;  distribution  is  unlimited. 


IV 


Figure  52.  IFT  SatComm  Tool  Integrated  With  RFI  Prediction  Results . 103 

Figure  53.  BER  Performance  With  CW  RFI  or  WB  RFI  and  Imperfect  Carrier  Tracking  Loop  104 

Figure  54.  Tracking  Error  Results  Under  No  RFI  and  CW  RFI  or  WB  RFI  Scenario . 104 

Figure  55.  Acquisition  Time  Result  Under  CW  RFI . 105 

Figure  56.  Acquisition  Time  Result  Under  WB  RFI . 106 

Figure  57.  Plots  of  Total  Carrier  Acquisition  Times  Versus  Loop  SNR  Without  RFI  and  With 

CW  RFI . 107 

Figure  58.  Plots  of  Total  Carrier  Acquisition  Times  as  Functions  of  Loop  SNR  Without  RFI  and 

With  WB  RFI  at  ISR  =  -10  dB . 108 

Figure  59.  Plots  of  Total  Carrier  Acquisition  Times  as  Functions  of  Loop  SNR  Without  RFI  and 

With  CW  RFI . 109 

Figure  60.  Plots  of  Total  Carrier  Acquisition  Times  as  Functions  of  Loop  SNR  Without  RFI  and 

With  WB  RFI  at  ISR  =  -40  dB . 110 

Figure  61.  Plots  of  PLL  Tracking  Jitter  in  the  Absence  of  RFI  and  Presence  of  CW  and  WB  RFI 

Signals  as  Functions  of  Loop  SNR  With  A/rfi  =  5  Hz . Ill 

Figure  62.  Plots  of  PLL  Tracking  Jitter  in  the  Absence  of  RFI  and  Presence  of  CW  and  WB  RFI 

Signals  as  Functions  of  Loop  SNR  With  A/rfi  =  10  Hz . 112 

Figure  63.  Plots  of  BER  due  to  RFI  for  Both  CW  RFI  and  WB  RFI  Signals  as  Functions  of  Loop 

SNR  With/RFi  =  5  Hz . 113 

Figure  64.  Plots  of  BER  due  to  RFI  for  Both  CW  RFI  and  WB  RFI  Signals  as  Functions  of  Loop 
SNR  With/RFi  =  10  Hz . 113 


Approved  for  public  release;  distribution  is  unlimited. 


V 


LIST  OF  TABLES 


Table  1.  Other  Aerospace  RFI  Tools . 12 

Table  2.  Typical  Parameters  for  USB  Waveform  Acquisition  and  Tracking  Modes . 22 

Table  3.  Carrier  Modulation  Mode  (CMM)  for  SATOPS  USB  Acquisition . 29 

Table  4.  Typical  Type-I  Imperfect  Integrator  Employed  By  NASA  PLL . 63 

Table  5.  Rain  Attenuation  Loss  in  dB  Calculated  for  2GHz  and  3  GHz . 98 

Table  6.  Signal  Losses  Due  to  Scintillation  for  15°  and  10°  Elevation  Angle  at  2.06  GHz . 98 

Table  7.  Signal  Losses  Due  to  Scintillation  for  15°  and  10°  Elevation  Angle  at  1.8  GHz . 98 


Approved  for  public  release;  distribution  is  unlimited. 


ACKNOWLEDGMENTS 


This  material  is  based  on  research  sponsored  by  Air  Force  Research  Laboratory  under 
agreement  number  FA9453-14-1-0259.  The  U.S.  Government  is  authorized  to  reproduce  and 
distribute  reprints  for  Governmental  purposes  notwithstanding  any  copyright  notation  thereon. 

DISCLAIMER 

The  views  and  conclusions  contained  herein  are  those  of  the  authors  and  should  not  be 
interpreted  as  necessarily  representing  the  official  policies  or  endorsements,  either  expressed  or 
implied,  of  Air  Force  Research  Laboratory  or  the  U.S.  Government. 


Approved  for  public  release;  distribution  is  unlimited. 


(This  page  intentionally  left  blank) 


Approved  for  public  release;  distribution  is  unlimited. 


1  SUMMARY 


This  report  describes  innovative  radio  frequency  interference  (RFI)  evaluation  and 
prediction  frameworks  that  incorporate  effects  of  RFI  on  the  satellite  receiver  synchronizers. 

This  report  also  presents  associated  RFI  modeling  and  prediction  tools  to:  (i)  detect  the  presence 
of  RFI  using  the  receiver’s  synchronizer,  (ii)  estimate  RFI  characteristics,  (iii)  evaluate 
effectiveness  of  the  existing  unified  S-band  (USB)  command  wavefonns  employed  by  civil, 
commercial  and  military  satellite  operations  ground  stations,  and  (iv)  predict  the  impacts  of  RFI 
on  USB  command  systems.  The  approach  presented  here  would  allow  a  communications 
designer  to  characterize  both  friendly  and  unfriendly  RFI  sources,  and  evaluate  the  impacts  of 
RFI  on  civil,  commercial,  and  military  USB  satellite  operations  systems.  In  addition,  the 
proposed  frameworks  and  RFI  prediction  models  would  allow  a  designer  to  estimate  the 
optimum  transmitted  signal  power  required  to  maintain  a  USB  satellite  operational  quality-of- 
service  (QoS)  in  the  presence  of  both  friendly  and  unfriendly  RFI  sources. 

2  INTRODUCTION 

SATOPS,  which  is  defined  as  a  satellite  operation  service  provided  by  a  satellite  control 
network  (SCN),  includes  tracking,  telemetry,  and  command  (TT&C)  services.  Figure  1  illustrates 
a  typical  civil,  commercial,  and  military  satellite  control  network.  It  contains  the  following 
components: 

•  Ground  Tracking  Station  (GTS):  Usually  located  at  a  remote  place  with  a  clear  view  of 
sky; 

•  Satellite  Control  Center  (SCC):  Provides  support  for  the  operation,  control,  and 
maintenance  of  satellites  -  SCC  executes  TT&C  operations; 

•  Network  Control  Center  (NCC):  Provides  support  for  the  scheduling,  satellite  network 
configuration,  and  operation  of  SCC,  GTS,  and  satellites;  and 

•  Ground  Communications  Network  (GCN):  Supports  data  transfers  among  SCC,  GTS, 
and  the  NCC. 

SATOPS  is  usually  operating  in  L-band  (1-2  GHz)  or  S-band  (2-4  GHz)  frequencies,  where 
there  are  many  potential  sources  of  radio  frequency  interference  (RFI)  to  TT&C  data  links.  We 
have  conducted  a  survey  of  RFI  analysis  tools  and  found  that  there  are  many  RFI  analysis  tools 
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available  to  analyze  the  impacts  of  RFI  on  SATOPS  systems.  However,  most  of  RFI  tools  do  not 
consider  the  time-factor  (TF)  when  calculating  bit  or  carrier  signal-to-noise  ratio  (SNR) 
degradation  due  to  interference.  Some  RFI  tools  incorporate  the  TF  by  incorporating  the 
SATOPS  schedules  along  with  interference  protection  criteria  (IPC1 2)  to  identify  the  blockage  due 
to  RFI.  These  tools  assume  that  RFI  increases  a  receiver’s  noise  floor.  None  of  these  tools 
consider  the  effects  of  RFI  on  the  receiver’s  synchronization  loops. 


Figure  1.  A  Typical  Satellite  Control  Network 


This  report  describes  an  approach  that  will  incorporate:  (i)  the  receiver  synchronization 
loops,  including  carrier  and  subcarrier  tracking,  and  symbol  timing  loops,  and  thereby  evaluating 
the  interfering  TF  at  the  synchronization  loop-level,  (ii)  statistical  technique  to  estimate  the  type 
of  interference  and  predict  its  impact  on  the  synchronization  loops,  and  (iii)  statistical  techniques 
to  optimize  the  SNR  of  the  transmitted  TT&C  signal  to  overcome  the  synchronization  errors 
caused  by  RFI  sources.  The  focus  is  on  the  SATOPS  using  USB  wavefonns  . 

Potential  radio  frequency  (RF)  sources  that  can  cause  interference  to  existing  civil, 
commercial,  or  military  USB  SATOPS  systems  can  be  classified  into  two  categories: 


1  IPC  is  Interference  Protection  Criteria  recommended  by  National  Telecommunication  and  Information 
Administration  (NTIA)  and  International  Telecommunication  Union  (ITU). 

2  USB  waveforms  are  used  by  civil  and  commercial  SATOPS  systems  operating  in  Unified  S-Band  frequencies. 
Recently,  U.S.  Air  Force  has  been  making  progress  to  migrate  its  SATOPS  operations  to  USB. 
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•  Category  I  -  friendly  RFI  sources:  The  sources  are  from  friendly,  unintentional 
interference  from  satellite  data  links  from  neighboring  orbiting  satellites.  This  category 
also  includes  friendly  SATOPS  tracking  stations  within  the  field-of-view  (FOV)  of  the 
victim’s  SATCOM  links.  Friendly  RFI  sources  can  cause  co-channel  interference  (CCI) 
or  adjacent  channel  interference  (ACI),  depending  on  the  transmitting  carrier 
frequencies.  Potential  friendly  RFI  sources  that  are  common  for  civil  and  commercial 
personal  communication  systems  include  3G  wireless,  4G  wireless,  WiFi  and  WIMAX. 

•  Category  II  -  unfriendly  RFI  sources:  These  are  caused  by  intentional  jammers, 
including: 

o  Classical  jammer:  A  high  power  jamming  signal  transmitted  with  the  intention  to 
shut  down  the  entire  USB  SATOPS  system  operation  or  to  significantly  reduce 
the  efficiency  of  USB  SATOPS  data  link  throughput;  and, 
o  Smart  jammer:  Purposed  to  monitor  the  USB  SATOPS  data  links  and  inject  weak 
jamming  signals  in  the  channel  with  the  intention  to  degrade  the  USB  SATOPS 
quality-of-service  (QoS). 

2.1  Description  of  Problems  for  RFI  Detection  and  Prediction 

Figure  2  illustrates  typical  friendly  RFI  scenarios  for  uplink  and  downlink,  respectively. 
The  figure  shows  that  potential  RFI  occurs  when  the  victim  receiver  (on-board  a  satellite  in  the 
case  of  an  uplink  or  at  the  ground  station  in  the  case  of  a  downlink)  simultaneously  collect  power 
from  a  desired  signal  source  and  other  undesired  interferers.  Existing  RFI  analysis  tools  evaluate 
the  extent  to  which  the  victim  receiver  collects  sufficient  power  from  one  or  more  interfering 
links  to  degrade  the  performance  below  an  acceptable  interference  protection  criteria  (IPC)  level, 
at  which  blockage  occurs.  The  problems  addressed  in  our  work  deal  with  the  detection  and 
prediction  of  the  RFI  signals  coming  from  both  friendly  and  unfriendly  sources. 
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Figure  2.  Typical  Friendly  RFI  Scenarios 


Our  research  developed  a  new  RFI  detection  algorithm  that  would  monitor  the  SATOPS 
receiver’s  synchronizer  performance.  The  RFI  detection  algorithm  assumes  that  the  interfering 
link  must  be  operating  and  satisfy  the  following  considerations: 

i.  Geometry  Consideration:  Victim  receiver  must  be  within  the  beamwidth  (or  field-of- 
view)  of  an  interferer’s  transmitting  antenna  and,  at  the  same  time,  the  interferer  must  be 
within  the  beamwidth  of  the  victim  receiver’s  antenna. 

ii.  Frequency  Consideration:  Any  interfering  link  must  have  part  or  its  entire  signal  power 
spectrum  within  the  victim  receiver’s  “interference  bandwidth”  (see  Figure  3).  Therefore, 
the  frequency  separation  between  the  victim  and  interfering  receivers,  and  the  spectral 
spreading  of  the  interfering  signal,  determines  the  contribution  to  the  interference  power. 
Likewise,  the  victim  receiver’s  “interference  bandwidth”  is  determined  based  on  the 
power  spectral  density  (PSD)  of  the  victim  signal.  The  spectral  spreading  of  the 
interfering  signal  depends  on  its  PSD,  which  in  turn  depends  on  its  waveform,  i.e.,  data 
rate,  coding  and  modulation  scheme. 
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iii.  Interfering  Power  Consideration:  From  the  total  interference  power  point-of-view,  the 
RFI  effects  from  all  the  interfering  links  (friendly  or  unfriendly)  must  exceed  the  IPC  to 
declare  that  a  potential  RFI  event  is  to  occur.  The  IPC  is  usually  derived  from  the 
victim’s  receiver  operating  threshold  or  obtained  from  national  or  international  regulatory 
bodies,  such  as  NTIA  (National  Telecommunications  and  Information  Administration), 
FCC  (Federal  Communications  Commission)  or  ITU  (International  Telecommunication 
Union).  Figure  3  describes  the  interference  bandwidth  and  the  corresponding  CCI  and 
ACI  scenarios. 

iv.  Interfering  Time  Consideration:  The  effects  from  all  the  interfering  links  must  be 
sufficient  from  the  time  duration  (i.e.,  time-factor  (TF))  point-of-view  to  cause  an 
unacceptable  degradation  in  victim  receiver  performance  in  order  to  declare  an  RFI  event 
has  been  detected.  The  TF  will  also  be  incorporated  into  the  receiver’s  synchronization 
loop  to  ensure  that  the  detected  event  is  truly  caused  by  an  RFI  source. 

The  RFI  detection  tool  development  and  analysis  framework  described  in  this  report  will 
incorporate  the  above  four  considerations  for  the  interference  links. 


Figure  3.  Definition  of  Interference  Bandwidth  and  Corresponding  CCI  and  ACI  Scenarios 


To  predict  RFI  occurrence,  the  RFI  analysis  framework  must  also  incorporate  the  four 
considerations  for  the  interference  links.  The  total  power  coming  from  the  interfering  links  must 
exceed  the  IPC  and  stay  within  victim’s  carrier  tracking  loop  and  symbol  synchronizer  loop 
bandwidth  long  enough  to  force  the  synchronization  loops  out-of-lock.  If  the  loops  can  re¬ 
acquire  and  track  the  desired  signal  quickly,  the  operator  may  not  declare  an  RFI  occurrence. 
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The  interfering  TF  consideration  is  a  complex  issue.  It  deals  with  the  carrier  tracking  loop 
and  the  symbol  synchronizer  perfonnance  (also  referred  to  as  timing  synchronization  loop, 
timing  synchronizer,  or  timing  recovery  loop).  In  addition  to  the  four  considerations  for  the 
detection  problem,  the  following  considerations  must  be  incorporated  in  the  RFI  framework  for 
the  prediction  problem: 

•  Carrier  Tracking  Performance  Consideration:  The  interference  signals  have  to  exceed  the 
carrier  threshold  power  level  and  stay  within  the  carrier  tracking  loop  bandwidth  for  a 
certain  time  duration  before  the  carrier  loop  drops  lock  and  the  carrier  loop  is  unable  to 
re-acquire  the  carrier  within  a  specified  time  window. 

•  Sub-Carrier  Tracking  Perfonnance  Consideration:  The  interference  signals  have  to 
exceed  the  subcarrier  threshold  power  level  and  stay  within  the  subcarrier  tracking  loop 
bandwidth  for  a  certain  time  duration  before  the  subcarrier  loop  drops  lock  and  the 
subcarrier  loop  is  unable  to  re-acquire  the  subcarrier  within  a  specified  time  window. 

•  Symbol  Timing  Synchronization  Perfonnance  Consideration:  The  interference  signals 
have  to  exceed  the  threshold  symbol  signal-to-noise  ratio  (SNR)  level  and  stay  within  the 
symbol  synchronizer’s  loop  bandwidth  for  a  certain  time  duration  before  the  symbol 
synchronizer  drops  lock  and  is  not  able  to  re-acquire  the  timing  within  a  specified  time 
window. 

1 .  Hop  /  De-Hop  /  Spread  /  De-Spread  Consideration:  For  a  victim  receiver  using  a 

wavefonn  employing  spread-spectrum  or  frequency-hopping  techniques,  the  interfering 
signals  must  stay  within  the  spreading  or  hopping  bandwidth  long  enough  to  cause  the 
de-spread  or  de-hop  synchronizer  to  drop  lock. 

2.2  Impacts  of  RFI  on  Synchronization  Loops  Performance 

This  section  elaborates  on  the  synchronization  loop  performance  considerations  described 
previously.  Figure  4  and  Figure  5  show  a  notional  second-order  Costas  loop  subcarrier  tracking 
the  jitter  performance  of  a  unified  S-band  (USB)  wavefonn,  namely  PCM  /  PSK  /  PM  (Pulse 
Code  Modulation  /  Phase  Shift  Keying  /  Phase  Modulation),  in  the  absence  of  and  presence  of 
RFI.  Figure  4  shows  that,  in  the  absence  of  RFI,  the  subcarrier  tracking  jitter  does  not  have  a 
break-point.  It  reaches  a  steady-state  value  around  1.2  degrees.  On  the  other  hand,  Figure  5 
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shows  that  the  tracking  jitter  grows  as  RFI  power  increases.  At  a  specific  interference  power  that 
exceeds  the  subcarrier  power  threshold  level  that  corresponds  to  a  threshold  subcarrier  tracking 
jitter,  e.g.,  15  degrees,  the  tracking  jitter  grows  logarithmically.  When  this  happens,  the 
subcarrier  tracking  loop  declares  that  an  RFI  event  has  occurred  and  the  loop  would  be  switched 
from  tracking  mode  to  acquisition  mode.  The  RFI  framework  must  be  developed  to  incorporate 
this  event. 


-  Second  Order  Costas  Loop 

-  One  sided  Loop  Bandwidth  =  7Hz 

Jitter  =  1.21° 
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Figure  4.  Notional  Subcarrier  Tracking  Jitter  Performance  of  PCM/PSM/PM  without  RFI 


Figure  5.  Notional  Subcarrier  Tracking  Jitter  Performance  of  PCM/PSM/PM  with  RFI 
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2.3  Scope  of  Effort 

The  objective  for  this  research  effort  was  three-fold: 

1 .  Investigate  radio  frequency  interference  and  the  effects  that  it  has  on  wireless 
communication  systems,  both  private  and  commercial,  both  terrestrial  and  for 
satellite  communications; 

2.  Develop  RFI  tools  that  can  detect  sources  of  interference,  predict  the  effects 
of  multiple  sources  at  different  frequencies,  and  allow  evaluation  of  the 
impact  to  operational  systems;  and 

3.  Asa  first-order  exercise  of  the  tools,  consider  satellite  operation  (SATOPS) 
systems  employed  by  civil  and  commercial  unified  S-band  waveforms. 

Our  work  focused  on  USB  waveforms  employed  by  private  and  commercial  service  providers  to 
support  military  and  commercial  SATOPS  applications.  To  achieve  this  objective,  the  following 
tasks  were  performed  and  completed  during  the  period-of-perfonnance  for  this  effort: 

i.  Task  #1:  The  Catholic  University  of  America  (CUA)  conducted  an  assessment  of  existing 
RFI  tools  related  to  SATCOM  and  SATOPS  and  then  selected  an  RFI  tool  that  could 
accurately  address  geometry  considerations,  frequency  considerations,  and  power 
considerations.  North  Carolina  State  University  (NCSU),  CUA’s  subcontractor,  provided 
support  for  this  task. 

ii.  Task  #2:  CUA  developed  a  framework  for  RFI  tool  development  applicable  to  both 
satellite  communications  and  satellite  operations. 

iii.  Task  #3:  CUA  developed  analytical  and  simulation  models  to  accurately  address 
interfering  time  duration  for  satellite  operation  applications.  The  focus  of  this  task  was  on 
the  civil  and  commercial  SATOPS  USB  wavefonns. 

iv.  Task  #4:  CUA  integrated  the  analytical  and  simulation  models  developed  in  Task  3  with 
the  RFI  tool  selected  in  Task  1.  A  beta  version  of  the  RFI  tool  for  USB  wavefonns  was 
completed. 

v.  Task  #5:  CUA  integrated  the  model  developed  in  Task  4  into  existing  Intelligent  Fusion 
Technology  (IFT,  a  sub -contractor)  SATCOM  tools  to  display  the  RFI  prediction  and 
detection  results.  IFT  provided  integration  and  Matlab  simulation  support  for  this  task. 
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vi.  Task  #6:  CUA  investigated  verification  and  validation  approaches  for  the  RFI  detection 

and  prediction  tools,  and  provided  recommendations  for  experimentally  validating  the 

modeling  and  assessment  tools. 

2.4  Organization  of  Report 

This  report  is  organized  as  follow: 

•  Section  3.1  provides  a  summary  of  the  results  found  in  Task  #1 .  Key  results  from  the 
survey  of  existing  RFI  analysis  tools  are  presented. 

•  Section  3.2  describes  the  framework  for  the  RFI  tools  developed  in  Task  #2.  This 
section  presents  the  proposed  innovative  framework  for  developing  analytical  and 
simulation  tools  to:  (a)  predict  the  characteristics  of  RFI  sources,  and  (b)  evaluate 
impacts  of  RFI  to  USB  SATOPS  systems.  The  framework  is  applicable  to  both 
satellite  communications  and  satellite  operations. 

•  Section  3.3  presents  the  analytical  and  simulation  models  developed  in  Task  #3. 

These  models  can  accurately  predict  the  interfering  time  duration  for  the  SATOPS 
applications.  As  mentioned  earlier,  the  focus  of  this  section  is  on  the  civil  and 
commercial  SATOPS  USB  waveforms. 

•  Section  4. 1  discusses  the  integration  results  obtained  in  Task  #4.  This  section 
describes  how  the  analytical  and  simulation  models  presented  in  Section  3.3  were 
integrated  with  the  RFI  tool  selected  in  Section  3.1. 

•  Section  4.2  describes  the  integration  results  obtained  in  Task  #5.  This  section 
discusses  how  the  model  was  ported  into  existing  IFT  (Intelligent  Fusion  Technology, 
Inc)  SATCOM  tools  to  display  the  RFI  detection  and  prediction  results. 

•  Section  4.3  presents  the  investigation  results  found  in  Task  #6.  This  section  presents 
the  numerical  results  obtained  from  the  analytical  and  simulation  models  presented  in 
Section  3.3,  and  discusses  the  verification  and  validation  approaches  of  these  models 
for  the  RFI  detection  and  prediction  tools  presented  in  Section  4.1. 

•  Section  5  and  Section  6  summarize  the  work  and  discusses  future  studies.  In  addition, 
this  section  also  provides  recommendations  for  experimentally  validating  the 
modeling  and  assessment  tools. 

Approved  for  public  release;  distribution  is  unlimited. 


9 


3  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 


3.1  Assessment  of  Existing  RFI  Tools  and  RFI  Tool  Selection 

A  thorough  survey  was  conducted  of  existing  RFI  analysis  tools  available  in  the  public 
domain.  As  listed  below,  many  RFI  analysis  software  tools  are  available  to  analyze  the  impacts 
of  RFI  on  digital  communication  systems.  As  pointed  out  earlier,  most  of  these  RFI  analysis 
tools  do  not  consider  the  time-factor  (TF)  when  calculating  the  bit  or  carrier  signal  to  noise 
degradation  due  to  interference.  These  tools  assume  that  the  RFI  power  simply  increases  the 
receiver’s  noise  floor  and  adds  to  the  amount  of  interference  power  collected  by  a  victim’s 
receiver.  However,  some  RFI  analysis  tools  did  consider  the  TF  by  incorporating  the  SATCOM 
and  SATOPS  schedules  along  with  the  Interference  Protection  Criteria  (IPC)  (recommended  by 
NTIA  and/or  ITU)  to  identify  the  blockage  (i.e.,  outage)  due  to  RFI,  but  these  tools  also  assumed 
that  the  RFI  power  causes  an  increase  in  the  receiver  noise  floor.  None  of  these  tools  considered 
the  effects  of  the  RFI  on  the  synchronization  loops,  including  de-spreading  loop,  de-hop  loop, 
carrier  tracking,  and  timing  loop,  and  incorporating  the  TF  at  the  synchronization  loop-level. 
Below  is  a  list  of  the  existing  RFI  tools  and  their  associated  capabilities. 

•  JPL  Enhanced  Spectral  Analysis  Tool  for  RFI  Analysis  and  Spectrum  Management: 
Provides  predictors  of  signal  /  interference  /  fdters  /  spikes,  power  spectrum  plots, 
calculations  for  spectral  power,  telecommunication  link  budget  tables  and  interference 
analysis  tables,  configuration  menu  for  the  analysis  table  that  contains  options  on 
antenna,  amplifier,  and  receiver  for  ground  terminal  characterization. 

•  NTIA  Spectrum  /  Orbit  Utilization  Program  for  Interference  Analysis  of  Geostationary 
Communication  Satellites:  This  tool  calculates  interference  and  computes  the 
communication  impacts  of  that  interference,  including:  (i)  the  feeder  link  carrier-to- 
interference  ratio  (C/I)  is  calculated  at  each  satellite,  and  (ii)  at  each  earth  station  the 
downlink  C/I  is  calculated,  taking  into  account  interference  from  all  other  satellites  in  the 
geostationary  orbit. 

•  AGI’s  Systems  Tool  Kit  (STK)  for  RFI  Analysis:  STK  models  satellite  dynamics  and  RF 
propagation,  communication  link  budget,  avoidance,  interference  and  jamming  in 
dynamic  “real-world”  scenarios. 
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•  JPL  Interference  Analysis  Program:  The  tool  is  used  to  select  the  radio  frequencies  that 
would  be  used  for  uplink  and  downlink  communications  with  the  objective  of  avoiding  or 
minimizing  the  possibility  of  RFI  between  existing  and  planned  missions. 

•  Lincom  LinCsim  (acquired  by  Titan  on  January  31,  2000;  Titan  was  later  acquired  by  L-3 
Communications  on  June  3,  2005):  Analytical  simulation  package  provides  modeling  of 
RFI  effects  on  bit  error  rate  perfonnance  including  detailed  models  for  the  space  shuttle 
and  Tracking  and  Data  Relay  Satellite  System  (TDRSS)  links. 

•  The  Aerospace  Corporation  Web-based  Tool  for  Analyzing  Effects  of  RFI  on  Satellite 
Communication  Systems:  The  tool  can  calculate  the:  (i)  interfering  frequency  factor  (the 
fraction  of  the  wave  that  is  within  the  victim’s  beamwidth)  on  a  victim  satellite,  (ii) 
interfering  satellite  attributes  including  altitude,  inclination,  earth  center  inertial  longitude 
and  orbit  phase  for  multiple  constellations  with  multiple  satellites,  and  (iii)  mean  and 
variance  of  RFI  durations  and  outputs  the  times  at  which  RFI  occurs  and  which  satellite 
and  ground  stations  are  involved. 

•  The  Aerospace  Corporation  RFI  Tool  for  Analyzing  IMT-2000  Interference  With  POD 
Satellites:  The  tool  uses  Satellite  Orbit  Analysis  Program  (SOAP)  to  generate  temporal 
latitude/  longitude  profiles  for  orbiting  satellites  and  calculate  interference  power  density 
using:  (i)  different  contributions  from  base  stations  in  rural  areas  and  cities,  (ii)  ITU 
propagation  path  loss  and  base  station  antenna  model,  (iii)  NTIA  traffic  models  for  IMT- 
2000  (derived  from  ITU  models),  (iv)  parameterized  percentage  of  IMT-2000  land 
coverage,  region  classification,  and  traffic  environment  mixes  for  cities  and  rural  areas, 
(v)  communications  link  calculation  models  for  CDMA-2000  and  W-CDMA,  and  (vi) 
actual  geometry  to  be  considered  in  computing  base  station  effective  isotropic  radiated 
power  (EIRP),  space  loss,  and  atmospheric  loss. 

In  addition  to  the  RFI  analysis  tools  described  above,  we  are  also  aware  of  other  RFI  tools  that 
were  developed  by  the  Aerospace  Corporation  in  the  mid-2000’s  for  both  commercial  and 
military  applications.  Table  1  summarizes  these  RFI  tools.  The  RFI  approaches  used  by  these 
tools  were  similar  to  the  above  tools  and  did  not  include  the  time-factor  using  synchronization 
loops.  This  effort  leveraged  the  Aerospace  Corporation  Web-based  Tool  described  above  for 
incorporating  the  RFI  detection  and  prediction  framework  and  models  to  assess  the  impacts  of 
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RFI  on  SATOPS  USB  systems.  The  newly  developed  RFI  analysis  tool  incorporats  the  effects  of 
the  carrier  synchronizer  loop  in  assessing  the  impacts  of  RFI  on  SATOPS  USB  systems. 


Table  1.  Other  Aerospace  RFI  Tools 


SW  Tool 

Analysis  Performed 

Task/Program 

GMSK/QPSK  with  pulsed  RFI 

ICO/Teledesic 

ACI  for  combinations  of  QPSK,  8PSK 

Gapfiller 

SPW  (Signal 

and  16QAM 

Processing 

Workstation) 

GMSK  with  a  variety  of  jammers 

AEHF 

QPSK  with  single-tone  jammer 

GBS 

Baseband  filtering  of  SGLS  to  reduce 
out-of-band  emissions,  and  SGLS 
interference  to  PCS 

AFSCN 

FORTRAN 

and  C 

Uplink  interference  in  USB 

SGLS  Migration  to 

Programs 

USB  Study 

3.2  Framework  for  RFI  Tool  Development  and  RFI  Assessment 

This  section  describes  the  framework  for:  (i)  RFI  tool  development,  (ii)  RFI  detection  and 
prediction,  and  (iii)  RFI  assessment,  respectively. 

3.2.1  Framework  for  RFI  Detection  and  Prediction  Tool  Development 

Figure  6  illustrates  a  framework  for  the  development  of  required  analytical  and 
simulation  tools  for  detecting  and  predicting  RFI  sources.  The  framework  is  color-coded  to  show 
the  models  and  tools  developed  under  this  effort  and  the  status  of  the  tool  development.  Each 
block  in  Figure  6  represents  a  model,  which  can  be  an  analytical  model  or  a  simulation  model, 
required  for  the  detection  and  prediction  of  RFI  sources.  The  analytical  and  simulation  models 
were  developed  and  integrated  according  to  the  proposed  framework,  and  implemented  in 
Matlab  to  provide  an  RFI  detection  and  prediction  tool  for  assessing  the  impacts  of  RFI  on 
SATOPS  systems.  The  blocks  labeled  “Completed”  means  that  a  model  was  developed  and 
implemented  in  Matlab.  The  lighter  blue  boxes  indicate  models  of  interest  for  this  phase  of  the 
study.  The  darker  blue  boxes  indicate  that  the  models  would  be  developed  when  time  permitted, 
and  the  darkest  blue  boxes  indicate  the  models  are  to  be  developed  in  the  next  phase  of  this 
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effort.  As  illustrated  in  Figure  6,  the  focus  of  this  study  was  on  two  USB  command  waveforms, 
namely,  PCM/PSK/PM  (pulse  code  modulation,  phase  shift  keying,  phase  modulation)  and 
PCM/PM/Bi-Phase  (pulse  code  modulation,  phase  modulation,  bi-phase)  [1,2].  Section  3.3 
describes  the  models  labeled  “Completed”. 


Start 


Completed 


Completed 


PCM/PSK/PM 
Signal  Model 


Completed 


PCM/PM/Bi- 
Phase  Signal 
Model 


|~  |  Focus  of  this  Phase 

□  When  Time  Permits 
M  Next  Phase 


Completed 


Develop 
AWGN  Model 


Develop  RFI 
Signal  Models 


Completed 


DTTL  =  Digital  Transition  Tracking  loop 
AWGN  =  Additive  White  Gaussian  Noise 


Completed 


Develop  USB 
Carrier  Loop 
Model 


Develop  USB 
Subcarrier 
Loop  Model 


Develop  USB 
Timing 
Recovery 
Loop  Model 


f  RFI  Detection  &  N 
Prediction 
Models 

Completed 


2nd-Order 
Costas  Loop 


2nd-Order 

DTTL 


BER  Performance  Model  Completed 


Imperfect 

Carrier 

Model 


Figure  6.  RFI  Detection-Prediction  Tool  Development  Framework 


3.2.2  Framework  for  RFI  Detection  and  Prediction 

As  mentioned,  the  focus  of  this  effort  was  modeling  of  RFI  detection  and  estimation  for 
the  carrier  synchronization  loop.  Figure  7  illustrates  a  framework  for  the  development  of  RFI 
detection  and  prediction  models.  The  framework  describes  the  analytical  and  simulation  models 
required  for  the  detection  and  prediction  of  the  RFI  sources  using  existing  USB  carrier  tracking 
loops.  As  indicated  in  this  figure,  the  carrier  synchronizer  loop  is  a  second-order  phase  locked 
loop  (PLL).  The  framework  also  describes  an  algorithm  to  integrate  these  analytical  and 
simulation  models  for  detecting  and  characterizing  RFI  sources  that  can  be  used  for  mitigating 
the  interference  effects  on  the  bit  error  rate  (BER)  performance  of  the  USB  SATOPS  systems. 
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Section  3.3  provides  a  detailed  description  of  the  required  analytical  and  simulation  models 
labeled  as  “Completed”. 


Start 


Completed 


£s  *(<> 


Assume 
PLL  Starts 
at  Tracking 
Mode 


Com  )leted 


BER  Performance  Model 


RFI 

Characteristics 

Characterized 


RFI 

Estimation  & 
Prediction 
Model 


RFI  Detection  &  Prediction  Models  for  USB  Carrier  Loo 


[End 


Figure  7.  RFI  Detection-Prediction  Models  Development  Framework 


The  block  entitled  “BER  Performance  Model”,  shown  in  Figure  7,  is  used  to  predict  the 
BER  performance  degradation  due  to  the  presence  of  RFI  and  is  further  expanded  and  illustrated 
in  Figure  8.  This  model  was  used  for  predicting  the  bit  SNR  performance  degradation  due  to  RFI 
for  the  specified  quality-of-service  (QoS)  in  terms  of  the  BER.  The  predicted  SNR  degradation 
included  the  effects  of  the  imperfect  carrier  tracking  due  to  additive  white  Gaussian  noise 
(AWGN)  and  RFI,  and  the  receiver  noise  floor  raised  by  the  RFI  power. 
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Figure  8.  BER  Performance  Model  for  Predicting  the  Impacts  of  RFI  on  SATOPS 

Performance 


3.2.3  Framework  for  RFI  Assessment 

Finally,  a  framework  for  assessing  RFI  impacts  on  overall  performance  of  the  SATOPS 
USB  system  is  presented  in  Figure  9.  Using  the  bit  SNR  degradation  found  from  the  “BER 
Performance  Model”,  this  framework  describes  a  process  to  overcome  the  effects  of  RFI  by 
increasing  the  transmitted  SATOPS  signal  power. 
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Figure  9.  RFI  Assessment  Framework 


3.3  RFI  Analytical  and  Simulation  Model  Development 

This  section  provides  a  detailed  description  of  the  required  analytical/simulation  models 
shown  in  Figure  6,  Figure  7,  Figure  8  and  Figure  9. 

3.3.1  Analytical  Models  for  Unified  S-band  Waveforms 

As  described  earlier,  the  PCM/PSK/PM  and  PCM/PM/Bi-Phase  waveforms  were  the 
focus  of  this  study.  Simplified  modulator  and  transmitter  models  for  these  two  waveforms  are 
presented  in  Figure  10  and  Figure  11,  respectively.  For  PCM/PSK/PM-Sinewave  Subcarrier 
waveform,  the  NRZ  (Non-Return-to-Zero)  command  data  is  PSK  modulated  on  the  subcarrier, 
with  subcarrier  frequency,  fsc,  set  at  8  kHz  when  the  command  bit  rate  is  less  than  2  kbps,  and  at 
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16  kHz  when  the  bit  rate  is  between  2  kbps  and  4  kbps.  The  modulated  command  subcarrier  is 
then  phase-modulated  on  the  carrier  frequency,^.  The  uplink  command  (or  tele-command)  S- 
band  carrier  frequency  is  selected  between  2025  MHz  and  2125  MHz.  The  PCM/PM/Bi-Phase 
wavefonn  is  used  when  the  command  data  rate  is  greater  than  4  kbps.  For  this  waveform,  the 
NRZ  command  data  is  converted  to  bi-phase  data  format  and  then  phase  modulated  on  the 
command  S-band  carrier  frequency.  The  modulation  index,  m,  for  these  two  waveforms  is 
selected  to  ensure  proper  power  is  allocated  to  the  carrier  signal  component  and  command  signal 
component  to  achieve  predefined  carrier  and  command  margins  [3].  Margin  is  defined  as  the 
difference  between  the  required  signal  power  and  received  signal  power.  A  detailed  description 
of  this  waveform  is  given  below. 


Simplified  PCM/PSK/PM  Modulator/Transmitter  Model  R(0 


Rs  =  Rb  =  Bit  rate  in  bps 


NRZ  Command 

(SA  (SA 

Phase 

High  Power 

Data  Source 

Modulator 

Amplifier  (HPA) 

Rs  <2  kbps  for  8  kHz 

Rs  <  4  kbps  for  1 6  kHz  sin(27l/sct))  m  yj2P^sm(2nfct)) 

0.1rad<m<1.4Rad 


fsc  =  8  kHz  or  1 6  kHz 


Uplink  (Earth  -  to  -  Space) :  2025 MHz  <  fc  <2125  MHz 


Figure  10.  PCM/PSK/PM  Modulator-Transmitter  Model 


Simplified  PCM/PM/Bi-Phase  Modulator/Transmitter  Model  /(/)  R(t) 


Rs  =  Rb  =  Bit  rate  in  bps 


NRZ  Command 

NRZ-to-Bi-Phase 

Phase 

High  Power 

Data  Source 

Converter 

Modulator 

Amplifier  (HPA) 

i?  >  4  kbps 

m 

0. 1  rad  <  m  <  1 .0  Rad 


Note  that  for  the  USB  Downlink  (Space  -  to  -  Earth) : 
2200  MHz  <y  <2290  MHz 


A/2I\sin(2jr/’t)) 


Figure  11.  PCM/Bi-Phase  Modulator-Transmitter  Model 
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Analytical  and  simulation  models  and  tools  are  required  to  detect  and  predict  the  RFI 
signal  characteristics  in  the  presence  of  these  two  USB  waveforms,  and  to  assess  the  impacts  of 
RFI  on  these  two  waveforms.  Using  the  USB  command  modulator/transmitter  models  presented 
in  Figure  10  and  Figure  1 1,  a  high-level  framework  for  the  RFI  detection-prediction  tools 
development  is  described  in  Figure  6.  Figure  6  describes  the  required  USB  signal,  noise  and 
interference  models  to  detect  and  predict  the  RFI  signals  in  the  presence  of  the  two  USB  signals. 
As  indicated  earlier,  Figure  6  is  color-coded  to  show  that  the  “USB  Carrier  Loop  Model”  and 
“Imperfect  Carrier  Model”  are  the  focus  for  this  phase  of  the  research. 

3. 3. 1.1  PCM/PS K/PM  Sine  Wave  Subcarrier  Waveform 

The  mathematical  expression  for  the  PCM/PSK/PM-sine-wave  subcarrier  is  defined 
by: 

S(t)  -  yjZP.  cos (2  T  f(t+  nxi(t)sm(2jrfxt))s  (1 

where 

PT  =  total  transmitted  power, 

fc  =  transmitted  carrier  frequency  (Hz), 

m  =  command  modulation  index  in  radians, 

d(t)  =  NRZ  command  data  with  symbol  rate  Rs,  and 

fsc  =  command  subcarrier  frequency  (Hz). 

Figure  12  illustrates  the  NRZ  data  format  used  by  the  PCM/PSK/PM-sine-wave 
subcarrier  waveform.  As  shown  in  this  figure,  the  data  format  is  a  non-return-to-zero  (NRZ) 
format.  Figure  13  shows  the  power  spectral  density  (PSD)  of  the  PCM/PSK/PM-sine-wave 
subcarrier  waveform.  The  subcarrier  is  used  to  separate  the  carrier  component  and  the  data 
component.  When  the  subcarrier  frequency  to  bit  rate  ratio,  n,  is  larger  (i.e.,  n  >  4),  the 
interference  from  the  data  component  to  the  carrier  tracking  loop  can  be  neglected. 
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NRZ-Level  or  NRZ  Data  Format 


Rs  =  Rb  =  Bit  rate  in  bps 
Ts  =  Tb  =  Bit  duration  in  second 

Figure  12.  NRZ  Data  Format  for  PCM/PS  K/PM-Sinewave  Subcarrier 


PCM/PSK/PM  Power  Spectral  Density 


/;,  =  9rs 


Figure  13.  Power  Spectral  Density  for  PCM/PSK/PM-Sinewave  Subcarrier  Waveform 

The  receiver-demodulator  for  the  PCM/PSK/PM-sine-wave  subcarrier  wavefonn  is 
shown  in  Figure  14.  This  figure  is  also  color-coded  to  show  that  the  focus  of  this  study  is  on  the 
carrier  synchronizer  loop.  This  report  focused  on  the  acquisition  and  tracking  perfonnance  of  the 
carrier  PLL  in  the  presence  of  RFI  sources. 
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PCM/PSK/PM  Sinewave  subcarrier 


Figure  14.  Receiver-Demodulator  for  PCM/PSK/PM-Sinewave  Subcarrier  Waveform 


3. 3. 1.2  PCM/PM/Bi-phase  Waveform 

Figure  15  illustrates  the  bi-phase  data  format  used  by  the  PCM/PM/Bi-Phase  waveform. 
The  mathematical  expression  for  the  PCM/PM/Bi-Phase  is  defined  by  [1]: 

Sft)  =  ftp  cos(2^;,/  +  md/t))  ^2) 

where  dft)  is  the  bi-phase  command  data  with  symbol  rate,  Rs,  or  bit  rate,  Rh-  Figure  16  shows 
the  PSD  of  the  PCM/PM-Bi-Phase  wavefonn.  The  spectral  null  created  by  using  the  bi-phase 
data  format  allows  the  communication  designer  to  insert  the  carrier  component  into  the 
wavefonn.  When  the  command  data  rate  is  greater  than  or  equal  to  4  kbps,  the  spectral  null 
becomes  “flat”  and  causes  negligible  degradation  to  the  carrier  PLL  performance.  The  receiver- 
demodulator  for  the  PCM/PM-Bi-Phase  waveform  is  shown  in  Figure  17.  Similarly,  this  figure  is 
color-coded  to  illustrate  that  the  focus  of  this  study  is  on  the  carrier  synchronizer  loop.  Table  2 
shows  acquisition  and  tracking  typical  parameters  used  by  the  SATOPS  USB  receiver- 
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demodulator  for  acquisition  and  tracking  of  the  carrier,  subcarrier,  and  timing  of  the  received 
waveforms. 


Symbol: 


level  A 

level  B 


level  A 

level  B 


10  0  10  110 

NRZ-Level  or  NRZ  Data  Format 


Bi-Phase  or  Split-Phase  Data  Format 


Rs  =  Rb  =  Bit  rate  in  bps 

Ts  =  Tb  =  Bit  duration  in  second 


Figure  15.  Bi-Phase  Data  Format  Used  By  PCM/PM-Bi-Phase  Waveform 


PCM/PSK/Bi-Phase  Power  Spectral  Density 


Figure  16.  Power  Spectral  Density  of  PCM/PM/Bi-Phase  Waveform 
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PCM/PM/Bi-phase 


Color  Code 


Focus  of  this  Phase 
When  Time  Permits 
Next  Phase 


Acronym: 

LNA  =  Low  Noise  Amplifire 

RF  =  Radio  Frequency 

BPF  =  Band-Pass  Filter 

IF  =  Intermediate  Frequency 

DTTL  =  Digital  Data  Transition  Tracking  Loop 


Figure  17.  Receiver-Demodulator  for  PCM/PM/Bi-Phase  Waveform 


Table  2.  Typical  Parameters  for  USB  Waveform  Acquisition  and  Tracking  Modes 


Trackinq  Mode 

Acquisition  Mode 

•  Carrier  Recovery  Loop  (PLL) 

•  Carrier 

-  Loop  Bandwidth  (BL):  10  -  400  Hz 

-  Acquisition  Loop  Bandwidth: 

B,Tu  <  0  1 

•  BLAcq  >  Doppler  Frequency 

-  Linear  Sweep  of  +  55  kHz 

•  SubCarrier  Recovery  Loop  (Costas) 

•  12.18  seconds  (duration) 

-  Loop  Bandwidth:  100  -  400  Hz 

-  Lock  Detector  Integration  time 

-  Damping  factor:  l/sqrt(2) 

~  300  symbols 

-  BLTb  <  0.1 

•  Sub-Carrier 

•  Timing  Recovery  Loop  (DTTL) 

—  Lock  Detector  Integration  time 

—  Loop  Bandwidth:  10  -  20  Hz 

~  256  symbols 

-  Damping  factor:  l/sqrt(2) 

•  Timing 

B,T\  <  0  1 

-  Lock  Detector  Integration  time 

UL1  b  —  1 

~  400  symbols 
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3.3.2  Analytical  Models  for  RFI  Signals 


The  received  signal  at  the  SATOPS  USB  receiver-demodulator  can  be  characterized  by: 

R(t)  —  S(t)  +  I(t)  +  n(t),  (3) 

where  n(t)  is  the  average  white  Gaussian  noise  with  one-sided  power  spectral  density,  No.  S(t)  is 
the  received  USB  signal  of  interest.  For  PCM/PSK/PM-sine-wave  subcarrier,  it  is  given  by: 

S(t)  =  a^2PT  cos(2 n{fc  +  fd )t  +  md(t)  sm(2fx,t)  +  0c )) 

S  (t)  =  a  ^2^  [exp  {j{2n(fc  +  fd)t+ 0C )).  exp  (j(md  (t)  sm(2rfsct))] 

S(t)  =  Real]  S~(t)  \ 

1  J  (4) 


where 


fd  =  Doppler  shift  due  to  satellite  movement, 

0C  =  received  carrier  phase  shift  of  the  desired  signal,  and 
a  =  channel  gain  factor. 

I(t)  is  the  RFI  signal  received  at  the  SATOPS  USB  receiver-demodulator,  and  it  is 
defined  as: 

7(0  =  ^2 fd,  (7)  cos(2/r(/;  +  A f)t  +  G, ) 

oo 

d,(t)=  Yj^rP^-nT,-  r), 

n=— oo 

where 


(5) 


Pi  =  RFI  power, 

A fc  =  RFI  frequency  offset  from  the  desired  carrier  frequency,^.,  and 
0i  =  RFI  phase. 

Here,  /„  and  p(t)  characterize  the  type  of  RFI.  The  RFI  signal  can  be  rewritten  in  the  following 
form: 

7(0  =  y/2P^dft)  Re[exp(j2n(fc  +  A fc)t  +  0;)],  (6) 

with/,,  the  RFI  data  sequence,/!^  the  RFI  unit  amplitude  rectangular  pulse  of  duration,  7/,  and  t 
the  time  asynchronization  of  the  RFI  with  respect  to  the  desired  signal,  s(t).  RFI  signals  can  be 
classified  into:  (1)  narrowband  RFI  or  continuous  wave  (CW)  RFI,  (2)  wideband  RFI,  and  (3) 
pulsed  RFI.  For  narrowband  RFI,  /„  =  1  and  7/  =  0.  For  wideband  RFI,  /„  is  an  independent  and 
identically  distributed  sequence  taking  on  equiprobable  values  of  ±1,  with 
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which  is  the  bit  duration  of  the  desired  signal  in  seconds.  For  pulsed  RFI,  In  =  1  during  77  and  In 
=  0  during  TPRF  -  77,  where  T\  is  the  RFI  pulse  duration  in  seconds,  and  PRF  is  the  pulse 
repetition  frequency. 

3.3.3  Simulation  Models  for  the  Second-Order  Phase  Locked  Loop 

A  typical  analog  PLL  model  used  by  advanced  NASA  spacecraft  and  satellite  is  shown  in 
Figure  18.  This  PLL  has  advanced  features  to  cut-out  the  loop  noise  and  smoothing  the  tracking 
phase  error  for  improved  carrier  tracking  perfonnance.  These  two  features  are  characterized  by 
the  low-pass  filter  B(s)  and  phase  smoother  V (s)  shown  in  Figure  18,  in  addition  to  existing  loop 
filter  F(s)  and  VCO  filter  K(s).  Figure  18  also  describes  the  characteristics  of  the  filters.  Figure 
19  describes  the  digital  modeling  of  the  advanced  second-order  PLL  shown  in  Figure  18  [4].  The 
figure  shows  detailed  description  of  the  filter  coefficients.  The  digital  closed-loop  transfer 
function  of  the  digital  filters  shown  in  Figure  19  is  defined  as  [4]: 

_  AI<  [B(Z)F(Z)V(Z)K(Z)\ 

-  1  +  AK  [B(Z)F(Z)V(Z)I<(Z)]  (g) 

The  one-sided  loop  noise  bandwidth  is  given  by: 

!*'=1  (9) 

Substituting  (8)  into  (9),  one  obtains: 
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where 


with 


Bdl  = 


1 


2TsH\l)\ 


aoboQo  ~  ooBiQi  +  a3B2Q2  —  oq B3Q3  +  B4Q4 


“o{(ao  “  al)Qo  ~  (“oai  -  “3“4)Qi  +  (aoa2  -  a2aA)Q2  -  (a0a3  -  aia4)<33} 


Bo  =  io  +  bf  +  b\  4-  63  +  b\;  B\  =  2(6o£>i  +  b\b2  +  62A3  +  6364) 

B2  —  2(b0b2  +  6163  4-  b2bi)\  B3  —  2(bob3  4-  6164);  B4  —  2iof> 4 
Qo  =  0061^4  —  a0a3e2  4-  04(0162  —  e3e4);  Q 1  =  000164  —  aoa2a3  4-  04(0102  —  o3e4) 

Q2  =  000162  —  aoa2ei  4-  04(0263  —  0362);  Q3  —  01(0162  —  e3e4)  —  02(0161  —  a3e3)  4-  a3(eie4  —  0362) 

Q4  =  a0  [e2(aia4  -  aoa3)  4-  e5(oo  -  oj)]  4-  (e2  -  e§)[ai(oi  -  a3)  +  (a0  -  a4)(e4  -  o2)] 

C\  =  ao  +  fl2>  e2  =  al  +  fl3>  e3  =  a2  +  fl4 
€4  =  do  +  fl4  >  €5  =  Oo  +  <*2  +  G4 

ao  =  2Fsj4ioj4h  +  AKAo;  ai  =  2Fs(j4ioi4oo-4i  +  ■Aoi-^ii^i  -  AioAnBi  -  -4ii>lio-4i)  +  AK(ZAo  -  Bo) 
a2  =  2Fs^i(^oo^oi  -  hoihn)  +  2FsBi(4ioj4ii  -  ^10^00  -  ^01^11)  +  3AK(Ao  -  Bo) 
a3  =  2FsB1(AooA10  -  A01A00  +  Aoi-An)  -  2Fs>li A0i>l0o  +  AK(A0  -  3 B0);  a4  =  2FsB1A10A0o  -  AKB0 
b0  =  AKA0;  61  =  j4A'(3j40  -  Bo);  h  =  3AK(A0  -  3 B0) 


(10) 


63  =  AK{A0  -  3B0);  b4  =  -AKBq 


(ID 
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INPUT 
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*(*)  =  £ 


V(S)  = 


1 

(1  +  tvS) 


Figure  18.  Advanced  Second-Order  PLL  Employed  by  NASA  Spacecraft  and  Satellites 


F(Z)  = 


(ApZ  -  Ba) 


K(Z)  = 


Ts(Z  +  1) 
2(2-1) 


Ts  =  Sampling  time  in  sec  =  — 
Fs 

Fs  =  Sampling  frequency  in  Hz 


B(z)  =  (1  +  z~') 

(AqqZ~1  4*  >4 ii ) 


Aq  —  1  4*  &o i  -^l  —  1  “I-  &o,  B0  —  do  —  1,  B\  —  bo  —  1 
Aqq  =  1  —  Co,  A\\  —  1  4-  Co 


Uo  Z-'+An) 


„  _  2trc  _  2r2  2 r, 

Co  Ts  ’  Ts  '  0  Ts 


A  —\  —  C 
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C  = 

'-'Ol 
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^10  —  1  '-01 

2r 


Figure  19.  Recursive  Implementation  of  F(z),  B(z),  K(z)  and  V(z)  Using  Bilinear 

T  ransformation 


The  sampling  frequency,  fs,  for  the  digital  filters  must  be  selected  to  satisfy:  fs  >  1  MHz. 
The  loop  bandwidth,  BDL,  for  tracking  mode  must  be  selected  to  satisfy:  BD[  TS  =  BmJb  <0.01. 
where  Tb  is  the  symbol  time  duration. 
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The  carrier  loop  SNR,  or  LSNR,  for  tracking  mode  in  the  absence  of  RFI  is  found  to  be: 

r2  / 


LSNR  = 


Jp(m)PT 

No  Bdl 


(12) 


The  carrier  LSNR  for  acquisition  mode  in  the  absence  of  RFI  is  defined  as: 


LSNRAcq  = 


J 0  (m)PT 

Nq  ^Acq-DL 


(13) 


A  typical  second-order  PLL  employed  by  civil  and  commercial  satellites  is  shown  in  Figure  20. 
As  shown  in  this  figure,  B(s)  =  V(s)  =  1,  and  the  loop  filter,  F(s),  is  usually  selected  as  a  Type-II 
(imperfect  integrator)  or  a  Type-Ill  (perfect  integrator)  filter.  Note  that  the  figure  also  defines  the 
voltage  controlled  oscillator  (VCO)  gain  in  terms  of  the  clock  frequency  and  the  number  of  bits 
in  the  accumulator.  For  imperfect  loop  filter  Type-II,  the  open-loop  response  G(s)  is  defined  as: 


G(s)  =  AK.K , 


F(s)  AK.K , 


VCO 


vco- 


1  +  sr. 


ST, 


(14) 


For  perfect  loop  filter,  the  closed-loop  response  H(s)  is  defined  as: 

G(s)  ( A  +  Bs)e~sT° 


H(s) 


where 


1  +  G(s)  s2  +  (A  +  Bs) 


a  =  tvAK.Kvco 
b  =  t2.AK.Kvco 


(15) 


(16) 


When  the  time  delay  is  zero,  the  closed-loop  filter  becomes: 

H(S>-  ,2,0 _ 2 


s  +2  gmNs  +  ml 


(17) 


The  loop  noise  bandwidth  is  found  as: 

Bl  =  Loop  noise  bandwidth  (BW)  in  Hz 


Bl  =  J|H(f)|V/  =  S7, 


1 

g  +  *g. 


(18) 


where 
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ujn  =  Natural  frequency  of  the  loop  in  rad/sec  =  Va  =  ^r].AK.KI 
g  =  Loop  damping  factor  = 


vco 


Lrco 


2tuN  2^jrl.AK.Kl 
For  perfect  loop  filter,  the  open-loop  response  G(s )  is  defined  as: 

m  _  ak.kvco 


G(s )  =  AK.K 


vco ■ 


- 1 -T-y 


,-sTn 


For  perfect  loop  filter,  the  closed-loop  response  H(s)  becomes: 

G(s)  ( A  +  Bs ) 


HAs) 


where 


1  +  G(s)  s2  +  (A  +  Bs) 


a_AK.Kvco 


B  =  —.AK.K 


vco 


Similarly,  the  loop  bandwidth  for  a  perfect  loop  filter  is  found  to  be  the  same  as  before: 

1 


BLP  =  JlH(f)l  2pdf  =  mA 


Sp 


4  ^ 


and 


or Np  =  Natural  frequency  of  the  loop  in  rad/sec  =  Va  = 

B  B 


AK.K 


vco 


C,,  =  Loop  damping  factor  = 


2ot 


n  2  \AK.Kyco 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 
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Imperfect  Integrator  loop  filter  Type-ll: 


K(s)  =  *hvco_ 
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V(s)=  1 


t2  =  Loop  filter  parameter 
Td  =  Pure  time  delay  of  the  loop 


fcLK  =  System  clock 

N  =  Number  of  bits  in  the  acummulator  =  32bits  (typical  accumulator) 


Figure  20.  Typical  Second-Order  PLL  Employed  by  Civil  and  Commercial  Satellites 


3.3.4  Existing  Carrier  Acquisition  Approaches 

Most  USB  Command  Detector  Unit  (CDU)  performs  carrier  acquisition  at  the  satellite  by 
going  through  four  carrier  modulation  mode  (CMM)  recommended  by  the  Consultative 
Committee  for  Space  Data  System  (CCSDS)  listed  in  Table  3. 


Table  3.  Carrier  Modulation  Mode  (CMM)  for  SATOPS  USB  Acquisition 

CARRIER  MODULATION 

MODE  (CMM)  DESCRIPTION 

CMM-1  UNMODULATED  CARRIER  ONLY 

CMM-2  CARRIER  MODULATED*  WITH  "ACQUISITION 

SEQUENCE" 

CMM-3  CARRIER  MODULATED  WITH  COMMAND  DATA 

CMM-4  CARRIER  MODULATED  WITH  "IDLE  SEQUENCE" 
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For  CMM-2,  the  “acquisition  sequence”  is  usually  an  alternating  sequence  (+1)  with  a 
predefined  length,  L,  in  terms  of  number  of  command  bits  N.  The  length  L  is  usually  selected 
using  the  following  relationship: 

L  =  2N.  (25) 

The  “idle  sequence”  is  usually  a  non-alternating  sequence. 

CCSDS  also  recommends  two  approaches  for  the  carrier  frequency.  The  approaches  are 
described  in  Figure  21.  The  main  difference  between  the  two  approaches  is  the  approach  for 
reacquiring  the  carrier  when  the  subcarrier  and/or  data  timing  synchronizers  drop  lock.  Approach 
1  assumes  the  subcarrier  and/or  data  timing  synchronizers  drop  lock,  but  the  carrier  synchronizer 
is  still  “in  lock”.  Approach  2  assumes  the  subcarrier  and/or  data  timing  synchronizers  drop  lock 
and  the  carrier  synchronizer  is  also  “drop  lock”.  The  carrier  acquisition  mode  includes  carrier 
frequency  acquisition  and  carrier  phase  acquisition  modes.  The  subsections  below  describe  these 
two  modes. 


Carrier  Acquisition  Approach  1 


Carrier  Acquisition  Approach  2 


Figure  21.  Carrier  Acquisition  Approaches  for  USB  SATOPS  Systems 
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3.3.4. 1  Carrier  Frequency  Acquisition  Mode 

The  most  commonly  used  technique  for  frequency  acquisition  mode  is  the  frequency 
sweep  technique  shown  in  Figure  22.  The  technique  consists  of  the  following  steps: 

•  Step  1:  Select  the  frequency  sweep  rate  based  on  the  estimated  received  carrier 
signal-to-noise  power  ratio,  or  CNR,  or  carrier  SNR. 

o  The  frequency  sweep  rate  as  a  function  of  received  signal  power  is  derived 
based  on  simulation  and  actual  operational  data. 

•  Step  2:  Uplink  command  carrier  frequency  is  swept  during  carrier  acquisition 
using  the  selected  sweep  rate. 

•  Step  3 :  When  the  command  signal  carrier  frequency  and  the  receiver  PLL  VCO 
frequency  coincide,  the  PLL  will  lock  onto  the  carrier  signal. 

o  Assuming  that  the  selected  sweep  rate  is  not  too  large. 

In  practice,  the  SATOPS  USB  receiver  utilizes  a  hard  limiter  in  front  of  the  PLL  to  control  the 
tracking  loop.  The  second-order  PLL  is  typically  employed  by  military,  civil  and  commercial 
satellites  for  tracking  the  carrier  signal  component.  The  minimum  sweep  rate  and  probability  of 
lock  for  the  loop  are  usually  obtained  by  simulation  using  actual  operating  conditions.  When  the 
carrier  is  acquired,  the  PLL  will  track  it  with  certainty,  i.e.,  with  probability  of  lock  greater  than 
90%.  Figure  23  shows  the  minimum  sweep  rate  and  probability  of  lock  as  a  function  of  the 
received  power. 
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Figure  22.  Carrier  Frequency  Acquisition  Using  Frequency  Sweep  Technique 
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Received  Power,  dBm 


Figure  23.  Minimum  Sweep  Rate  and  Probability  of  Lock  as  a  Function  of  Receive  Power 


33.4.2  Carrier  Frequency  Acquisition  Time  Using  Sweeping  Technique 

For  the  frequency  sweep  technique,  the  carrier  frequency  acquisition  time,  Tacq,  is 
developed  using  the  minimum  sweep  rate  specified  in  Figure  23  and  simulation  results  using 
actual  operational  conditions  from  practical  SATOPS  missions,  and  it  is  found  to  be: 

_  A a>un  _  2nfun 
—  Agj  _  A<i)  ’ 

where  fUn  is  the  frequency  uncertainty  due  to  Doppler  and  drift  of  the  voltage  controlled 
oscillator  and, 
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A  co  =  Sweep  rate  in  Hz/sec  =  co 


1- 


1 


Loop  SNRAcq  -2  j 


,  for  3  <  Loop  SNRAcq  <4.75 


Aa>  =  < 


co: 


co 


1- 


Loop  SNR  Acq  -  4  y 


,  for  6  <  Loop  SNR  <  9.5 


,  for  Loop  SNR  Acq  >  9 . 5 


2  2BLAcq 
= 


1  +  ^T 


,4"  =  0.707 


BLAcq  =  Acquisition  loop  bandwidth 


Note  that  during  acquisition,  the  loop  SNR  is  defined  as: 

■J  (i  ( >n )  Pr 


Loop  SNR  Acq  = 


NnB 


0  L Acq 


(27) 


(28) 


(29) 


(30) 


Note  that  without  loss  of  generality,  one  can  assume  the  channel  gain  factor  is  equal  to  unity. 


33.4.3  Carrier  Phase  Acquisition  Mode 

For  carrier  phase  acquisition  mode,  the  most  commonly  used  techniques  in  practice  are 
carrier  phase  sweep  and  synchronization  (or  sync)  word.  The  carrier  phase  acquisition  process  is 
described  in  Figure  24  and  is  briefly  described  below: 

To  =  Time  required  for  the  PLL  to  estimate  the  initial  phase  Go,  and 

Te  =  Time  it  takes  the  PLL  to  reach  steady  state  phase  error  9e. 

The  carrier  phase  acquisition  time,  TpAcq,  is  calculated  using  the  following  relationship: 

TpAcq  =  T0  +  Te  .  (31) 

The  carrier  phase  acquisition  time  improves  significantly  with  phase-aided  acquisition  using  the 
sync  word,  or  phase  sweeping  technique,  to  estimate  initial  phase.  The  discussion  of  these 
techniques  is  provided  below. 
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Figure  24.  Phase  Acquisition  Process  of  a  PLL 


3.3. 4.4  Carrier  Phase  Acquisition  Using  Sync  Word  Technique 


Typical  carrier  phase  acquisition  technique  is  to  use  sync  word  (or  alternating  sequence 
or  acquisition  sequence)  to  estimate  initial  phase.  This  subsection  assumes  that  the  Doppler 
frequency  and  timing  are  known  in  advance  and  that  the  command  data  is  formatted  as  an  NRZ 
rectangular  data  pulse,  i.e., 


Pit)  = 


J  1,  0  <  t  <  Tb 
[0,  elsewhere 


(32) 


It  is  also  assumed  that  the  sync  word  or  acquisition  sequence  with  length  L  is  known,  and  the 
sync  word  is  transmitted  during  CMM-2  (see  Figure  20  and  Figure  21),  using  direct  modulation 
with  word  sequence.  Mathematically,  the  carrier  signal  modulated  with  an  acquisition  sequence 
can  be  expressed  as: 
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^CMM- 2  (0  =  yf^dAcq  (t  )[exp  (j{ln(fc  +fd)t  +  6C ))] 
ScMM-2  (0  =  ^eal|“^CMM-2  (0 '  j 

d  Acq  (0  =  YjCn.p{t-nTb), 

n=-L/ 2 

where 


(33) 


/d  =  Doppler  shift  due  to  satellite  movement, 

6C  =  Received  carrier  phase  shift  of  the  desired  signal,  and 
Cn  =  Acquisition  sequence  with  length  L. 

The  carrier  signal  modulated  with  an  acquisition  sequence  or  sync  word  sequence  is  received  at 
the  SATOPS  USB  receiver-demodulator  and  processed  by  a  typical  Maximum  Likelihood  (ML) 
phase  estimator  shown  in  Figure  25. 


Figure  25.  Typical  Maximum  Likelihood  (ML)  Phase  Estimator 

The  performance  of  the  ML  phase  estimator  shown  in  Figure  25  is  expressed  in  terms  of 
the  variance  of  the  estimated  carrier  phase,  which  is  given  by  [5]: 


cr  =var  \6C-0C\  = 


1  N0 

rn 

1 

2 L  (2J;(m)PTTb) 

U  L) 

Bit  SNRAcq 

where 


Bit  SNR 


Acq 


2  Jl(m)PTTb 
Nn 


(34) 


The  above  equation  as  a  function  of  bit  SNR  is  only  for  acquisition  mode,  and  without  loss  of 
generality,  one  can  assume  the  channel  gain  factor  is  equal  to  unity.  For  a  given  bit  SNR  and  a 
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specified  variance  of  the  estimated  carrier  phase,  the  length  of  the  sync  word  sequence  is 
detennined  using  Equation  (34). 


3.3. 4.5  Carrier  Phase  Acquisition  Time  Using  Sync  Word 


The  length  of  sync  word  must  be  selected  to  ensure  the  variance  of  the  carrier  phase 
estimate  error  is  small  such  that  the  time  it  takes  the  PLL  to  reach  steady-state  phase  error,  0e ,  is 
about  l/2f?Acq-DL-  The  acquisition  models  presented  here  ensure  the  probability  of 
synchronization  failure  at  1 0s.  The  phase  sync  work  can  be  computed  from: 


where 


T 

Phase  Sync 


2  B 


■  +  LT. 


LAcq 


(35) 


L  =  length  or  number  of  bits  in  the  sync  word, 
Tb  =  sync  work  bit  duration  in  seconds,  and 
5LAcq  =  acquisition  loop  bandwidth  in  Hz. 


3. 3. 4.6  Carrier  Phase  Acquisition  Time  Using  Phase  Sweeping  Technique 


Based  on  the  simulation  results  using  the  various  sources,  the  carrier  phase  acquisition 
time  using  phase  sweeping  technique  is  found  to  be: 


T  =  < 

J  Phase_Swp  ' 


1  +  ,  for  Loop  SNR.  =  10  dB 


25 LAcq  45 i^cq 


25 LAcq  45 LAcq 


,  for  Loop  SNRAcq  =  14  dB 


The  total  carrier  acquisition  time,  Tacq,  is  calculated  using  the  following  relationship: 


T  =T  +T 

±  Acq  Freq  ±  Phase 


(36) 

(37) 


where  rFreq  is  defined  in  (26),  and  rPhase  is  defined  as  in  (35)  or  (36),  depending  on  the  phase 
sync  technique  used  by  the  SATOPS  USB  receiver. 


3.3.5  Analytical  Models  for  Assessing  PLL  Acquisition  Performance  in  the  Presence  of 
RFI 

This  section  presents  analytical  models  for  evaluating  the  impacts  of  CW  and  wideband 
RFI  signals  on  the  second  order  PLL  acquisition  perfonnance. 
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3. 3. 5.1  Analytical  Models  ofPLL  Acquisition  with  CW  RFI 

Recall  from  previous  section,  the  CW  RFI  signal  is  defined  as: 

7(0  =  a/2^  cos(2tt(/c  +  A  fc)t  +  0,),  (38) 

where 

Pi  =  RFI  power, 

A fc  =  RFI  frequency  offset  from  the  desired  carrier  frequency,/;,  and 
9i  -  RFI  phase. 

For  continuous  wave  (CW)  RFI,  the  effects  of  RFI  depend  on  the  RFI  carrier  frequency,  namely, 
case  1  for  out-of-band  interference,  and  case  2  for  in-band  interference.  The  in-band  interference 
occurs  when  the  RFI  carrier  frequency  is  within  the  carrier  acquisition  loop  bandwidth,  RLACq, 
and  the  out-of-band  interference  happens  when  RFI  carrier  frequency  is  outside  the  carrier 
tracking  loop  bandwidth. 


3. 3. 5. 1.1  Case  1:  A  fc>  Bum 

For  case  1,  the  loop  SNR  in  the  presence  of  CW  RFI  is  given  by: 

LSNR. 

LSNRr  ~  1 


VCW- 1 


where 


LSNRAcq  = 


N0BLAcq 


*  Acq- 1 


and  A*ri  = 1 


(39) 


(40) 


For  case  1,  the  carrier  acquisition  time,  Tacq,  remains  the  same  as  the  non-interference  case 
presented  in  the  previous  section. 


3.3. 5. 1.2  Case  2:  Afc<  Bum 


For  case  2,  the  loop  SNR  in  the  presence  of  CW  RFI  is  given  by: 


LSNR, 


LSNR 


‘ Acq 


CW- 2 


A 


Acq -2 


where 


^Acq—2  —  [1+INR] 


(41) 

(42) 
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and  INR  is  the  ratio  of  Pit  No,  the  interference  power-to-noise  ratio.  For  case  2,  the  carrier 
acquisition  time,  racq.Cw-2,  will  be  longer  due  to  the  loop  SNR  degradation  by  a  factor  AAcq-2-  The 
frequency  sweep  rate  will  be  calculated  using  LSNRcw-i,  i.e.. 


A(yCPF_2  —  coN 


1- 


1 


,1/2 


Loop  SNR 


_  9 

CW-2  J 


for  3  <  Loop  SNR cw  2  <4.75 


(43) 


A  cocw_2  —  < 


co : 


,  1/2 


1- 


v  LSNRCW_2  -4 j 


,  for  6  <  LSNRcw.2  <9.5 


^,for  LSNRcw_2  >9.5 


(44) 


2  ^^LAcq 

coN  = 


l+x 

c2. 


,  4"  =  0.707 


BLAcq  =  Acquisition  loop  bandwidth 


(45) 


Again,  without  loss  of  generality,  one  can  assume  the  channel  gain  factor  is  equal  to  unity.  The 
carrier  frequency  acquisition  time  becomes: 

_  Jn 


L  Freq-CW-2 


A  co, 


CW-2 


(46) 


Similarly,  for  case  2,  the  carrier  phase  acquisition  time  becomes: 


T  -  < 

1  Phase_Swp-CW-2  ' 


1  +  12-2  ,  for  LSNR  cw_2  =  1 0  dB 


Acq-DL  4S Acq-DL 


Acq-DL  4S Acq-DL 


,  for  LSNR  cw_,  =14  dB 


(47) 


where 


1 


[  Phase_Sync-CW—2 


2  B 


+  7  T 

^  ^ CW-2 1  b 


LAcq 


Lcw  2  =  Required  length  or  number  of  bits  in  the  sync  word  in  the  presence  of  CW  RFI 
Tb  =  Sync  word  bit  duration  in  second 

For  CW  RFI  case  2,  the  performance  of  the  ML  phase  estimator  becomes: 


.  (48) 


cr 


f  1  1 

(l  +  INR) 

-CW-2 

IT 

\  ^cw-2  y 

'  BitSNRAcq_ 

(49) 
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where  bit  SNR  is  defined  as  before,  i.e.,: 


Bit  SNRAcq  = 


2J-(m)PTTh 


To  compensate  for  the  CW  RFI,  one  can  increase  the  length  of  the  sync  word  or  the  bit  SNR.  For 
the  phase  sweep  technique,  the  carrier  acquisition  time  in  the  presence  of  CW  RFI  becomes: 


T  =  T  +T 

±  Acq-Sweep-CW—2  Freq—CW—2  ±  Phase_Swp-CW—2 


For  sync  word  technique,  the  carrier  acquisition  time  is  given  by: 


T  —  T  +T 

Acq—  Phase-Sync—  CW— 2  Freq—CW—2  Phase_Sync-CW—2 


33.5.2  Analytical  Model  of  PLL  Acquisition  with  WB  RFI 

Recall  from  the  previous  section,  the  wideband  (WB)  RFI  signal  is  defined  as: 
Ht)  =  ^[2P1  d,(t)cos(2n(fc  +  A  fc)t  +  0,), 


where 


di(t)  =  £“=_oo 4  ■  p(t  -  nf;  -  r), 


Pi  =  RFI  power, 


Ti=  RFI  bit  duration  in  seconds,  and  7)  <  T/h 
A fc  =  RFI  frequency  offset  from  the  desired  carrier  frequency,/;,  and 
T\  =  RFI  phase. 

The  loop  SNR  in  the  presence  of  wideband  (WB)  RFI  can  be  shown  to  have  the  following  fonn: 

LSNR4 

LSNRu/d  = - ^ 


where 


LSNRAcq  = 


*  Acq -WB 


■/)('»)  4 

N0BLAcq 


^Acq-WB  ~  1+  $ rfI-WB^BdL’  fun)  ’  ~  INR 

^ 0  J  ^0 

sm.„  (/,/»)  =  [s,  (/  -  fu. )  +  S,  (/  +  fu. )] 


W)=r,  j 


BAcq-DL  V  ■,  _jrp  N  1 


sin (nfT,) 


df,  where  Tj  =  l/7?; 
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Again,  fjn  is  defined  as  the  frequency  uncertainty  due  to  Doppler  and  VCO  drift. 

For  the  phase  sweep  technique,  the  carrier  acquisition  time  in  the  presence  of  WB  RFI 
becomes: 


T  —  T  +  T 

±  Acq—  Sweep—  WB  Freq—WB  ±  Phase_Swp—WB 


For  sync  word  technique,  the  carrier  acquisition  time  is: 

T  =  T  +T 

±  Acq—  Phase- Sync-WB  Freq—WB  ±  Phase_Sync-WB 


where: 

/  _  --rflJn 

1  Freq-WB 

A  coWB 

n 

where: 


co : 


co. 


,1/2 


1- 


vLSNRWB-4y 


J,  for  6  <  LSNR  WB  <  9.5 


-y,for  LSNR  WB>  9.5 


T 

1  Phase_Swp-WB 


for  LSNR  WB  =  10  dB 

2  BLAcq  4BLAcq 

— - —  +  — ,  for  LSNR  WB  =  14  dB 

2  BLAcq  4  BLAcq 


1 


1  Phase  Sync-WB 


2  B 


.+  T  T 

^  ^  WB 1  b 


LAcq 


Lwb  =  Required  length  or  number  of  bits  in  the  sync  word  in  the  presence  of  WB  RFI 
Tb  =  Sync  word  bit  duration  in  second 


For  WB  RFI,  the  performance  of  the  ML  phase  estimator  becomes: 


{  1  1 

^  Acq -WB 

■WB 

IT 

\  wb  y 

BitSNRAcq_ 

<7" 


^  Acq— WB  ~  [l  +  INR.S RFI  _WB  (B Acq 

-DL  ’  fun  )] 

Sm  fun  )  =  [St  (/  -  fun  ) +  S ,  (/  +  fUn  )] 


W)=^  J 


Acq-DL 

f 

sin (ytfTi) 

J 

Acq-DL 

xfT, 

df ,  where  T:  =  \  /  R, 


(57) 

(58) 


(59) 


(60) 


(61) 


(62) 


(63) 


(64) 
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For  the  ML  phase  estimator,  the  acquisition  loop  bandwidth,  BAcq_DL,  is  identical  to  IWb-  To 
compensate  for  the  WB  RFI,  one  can  increase  the  length  of  the  sync  word  or  the  bit  SNR. 

3.3.6  Analytical  Models  for  Assessing  PLL  Lock  Detector 

Figure  26  shows  how  a  PLL  lock  detector  is  integrated  into  the  PLL  operation  for 
acquiring  and  tracking  the  carrier  signal.  Note  that  the  parameters  BDL  and  BAcq.DL  shown  in 
Figure  26  denote  the  digital  loop  bandwidth  for  tracking  and  acquisition  mode,  respectively, 
assuming  that  the  second-order  PLL  is  implemented  in  the  digital  domain.  Regardless  of  whether 
the  loop  is  implemented  in  digital  domain  or  analog  domain,  the  loop  bandwidth  should  be  set 
the  same  for  both  cases,  i.e.,  BD l  =  B L  and  BAcq  =  BAcq.DL,  where  B\  and  BAcq  are  the  analog  loop 
bandwidth  for  tracking  and  acquisition  modes,  respectively. 

As  shown  in  Figure  26,  the  PLL  initially  will  be  in  the  acquisition  mode,  acquiring  the 
carrier  signal  component.  In  this  mode,  the  loop  bandwidth  will  be  open  wide  enough  to  cover 
the  range  of  the  expected  frequency  uncertainty  of  the  carrier  signal,  normally  110  kHz  for 
SATOPS  applications.  During  this  mode,  the  phase  lock  loop  will  monitor  the  carrier  phase  error 
over  an  observation  time  If.  If  the  phase  error  is  less  than  a  threshold  values  a-rhcar,  then  the  PLL 
will  switch  to  the  tracking  mode.  For  the  tracking  mode,  the  PLL  tracking  loop  bandwidth  will 
be  set  according  to  the  loop  bandwidth  product  f?LIb  <0.01,  where  It,  is  the  command  data  bit 
duration.  This  section  describes  how  the  threshold  values  axhcar  and  II  are  set  for  the  optimum 
lock  detector  performance. 
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Figure  26.  Typical  PLL  Lock  Detector  Employed  By  Civil,  Commercial  and  Military 

SATOPS  Satellites 


The  value  a-rhcar  is  selected  to  ensure  the  probability  of  lock  is  greater  than  90%. 
Theoretically,  it  is  chosen  by  using  the  following  relationship: 


®Thcar 


L JLSNREff 


+  ®Margir 


(Radians'). 


(65) 


The  effective  loop  SNR  (LSNReg)  shown  in  (65)  is  defined  as  the  loop  SNR  in  the  presence  of 
noisy  carrier  reference  for  a  specified  BER  performance.  The  LSNRe ff  can  be  calculated  using 
(70)  shown  below,  and  aMargin  is  the  phase  error  margin  that  is  specified  by  the  communications 
designer  to  take  into  account  the  phase  jitter  caused  by  VCO  drift,  residual  Doppler,  phase  noise, 
etc.  In  practice,  the  phase  error  margin  can  be  5%  to  10%  of  the  value  of  \N(LSNRes). 

The  value  a-rhcar  is  selected  using  the  effective  loop  SNR  that  is  evaluated  in  the  absence 
of  RFI.  It  is  assumed  that  the  PLL  tracking  error  follows  the  Tikhonov  distribution  and  that  the 
carrier  tracking  phase  error  is  small.  By  using  Taylor  series  expansion  on  the  average  uncoded 
BER  equation,  one  can  approximate  the  uncoded  BER  as  [6]: 
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beR*l  =  - 


where 


=  ^erfc{jBSNR0}+±j 


-J-BSNR()cry{BSm} 
n 


BSNR0  =2 

N n 


1 


<ya  = 


(66) 


As  before,  if  one  defines  the  PLL  tracking  loop  bandwidth  as  BL  and  the  carrier  tracking  loop 
SNR,  LSNR0  as: 

LSNR0=J20(m )  Pr 


N0  Bl 


(67) 


then  the  uncoded  BER,  in  the  absence  of  RFI  with  small  carrier  tracking  error,  can  be  rewritten 
as  a  function  of  the  carrier  tracking  loop  SNR,  LSNRq: 


BERfti  = 


\  erfc^aiBjjLSNR,  }- 


1  ^ a{BLTb).LSNR0 


-{aiBJ-J.LSNS,} 


LSNRn 


where 


a  —  2 


'  Jo  (Pi) 


(68) 


(69) 


The  first  tenn  of  (68)  represents  the  ideal  BER  perfonnance  with  perfect  carrier  tracking,  and  the 
second  tenn  is  the  BER  degradation  due  to  imperfect  canier  tracking  due  to  the  presence  of 
AWGN.  The  imperfect  carrier  tracking  is  represented  by  the  carrier  tracking  jitter,  Oq.  This  term 
is  also  defined  as  the  “Radio  Loss”  (RL).  From  (68),  the  effective  loop  SNR  is  defined  as: 

Effective  Loop  SNR  =  LSNREff  =  a{B,Tb).LSNR() 

Substituting  (70)  into  (68),  the  uncoded  BER  becomes: 


(70) 


BERm  -  — 

m.  2 


LSNREff  1+ 


a 


{BLTb\ 


LSNR 


71 


Eff 


-{. LSN 


LSNR, 


v« r  .  (71) 

The  effective  loop  SNR  must  be  chosen  for  the  uncoded  BER  to  meet  the  desired  threshold  BER 
value,  e.g.,  BERThrcshokl  =  10"9  or  10'6. 
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The  effective  loop  SNR  must  also  be  chosen  to  ensure  that  the  mean  time  between  cycle 
slips,  Jmbcs,  is  sufficiently  large  that  probability  due  to  cycle  slip  occurring  due  to  noise  can  be 
ignored: 

T  n*  ^  {x-LSNRej  } 

1MTBCS  ~  R  ^ 

.  (72) 

To  ensure  a  small  probability  of  false  lock  at  approximately  10'4  [6],  the  integration  time  T\  of 
the  PLL  detector  can  be  set  as  follows: 


0.018  *LSNREff*BL 

Effective  loop  SNR  must  be  set  to  meet  the  specified  uncoded  BER  performance  and  mean  time 
between  cycle  slips.  The  above  equations  are  used  to  determine  the  required  carrier-to-noise 
power  spectral  density  (C/N0). 


3.3.7  Analytical  Models  for  Assessing  PLL  Tracking  Performance 

This  section  describes  analytical  models  for  evaluating  the  impacts  of  CW  and  wideband 
RFI  signals  on  the  PLL  tracking  performance. 


3.3. 7.1  Analytical  Model  of  PLL  Tracking  Jitter  in  the  Absence  of  RFI 


The  performance  of  the  PLL  is  characterized  by  the  variance  of  the  tracking  phase  error, 
erg  ,  and  in  the  absence  of  RLI,  it  is  given  by: 


where 


A 

LSNR 


(74) 
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LSNR  =  Loop  SNR  = 

Pc  =Jq  (m)Pr 

i+^sCDm 


c  N ote  that  Pc  /  N0  is  the  carrier  SNR 


n0bl 


A  = 


Nn 


sCD(k,f )  =  -  Wsc)  +  W  +  *fsc)\ 


k  =  \,k=odd 


Br 


SAf)  =  Tt\ 


sin  (7frh) 
*fTh 


df, where  Th  =  1  /  Rh 


(75) 


Note  that  without  loss  of  generality,  one  can  assume  the  channel  gain  factor  is  equal  to  unity. 
Note  that  whcn  /^c  >  4i?b,  (Joc  becomes: 

<y2  ~  1 

~  LSNR  .  (76) 

In  practice,  for  reliable  tracking,  the  loop  SNR  is  usually  set  at  about  1 1  dB  or  more  to  ensure 
reliable  tracking  performance. 


33.7.2  Analytical  Model  of  PLL  Tracking  Jitter  with  CW  RFI 

A  continuous  wave  (CW)  RFI  signal  is  defined  as: 

/(t)  =  cos(2n(fc  +  A fc)t  +  0,),  (77) 

where 

Pi  =  RFI  power, 

A fc  =  RFI  frequency  offset  from  the  desired  carrier  frequency,  fc,  and 
0[  =  RFI  phase. 


33.7.2.1  Case  1:  Afr  >  Bnr. 

For  Case  1,  the  variance  of  the  tracking  phase  error,  Oqc ,  in  the  presence  of  CW  RFI  can 
be  shown  to  have  the  following  fonn: 

<j2  =  — 

LSNR  .  (78) 

Here,  A  is  defined  as  in  (74),  and  for  fsc  >  4 Rb,  A  =  1 . 
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3. 3. 7. 2. 2  Case  2:  Afr  <  Bnr 


For  Case  2,  the  variance  of  the  tracking  phase  error,  <Jqc,  in  the  presence  of  CW  RFI  can 
be  shown  to  have  the  following  fonn: 


„2  _  Acw 

CT  /)  — 

'  LSNR 
ACPF=[l  +  INR] 
p 

INR  =  — —  =  Interference  -  to  -  noise  p  ower  ratio 

N°  .  (79) 


Here,  one  can  assume  that  /C  >  4i?b-  Note  that  when  INR  is  greater  than  the  carrier  SNR,  the  PLL 
will  drop  lock  on  the  command  signal  and  lock  onto  the  CW  RFI  signal.  Note  that  for  this  case, 
the  effective  loop  SNR,  LSNREff,  is  calculated  from  (79),  using:  LSNREff  —  Using  this 

J J  11  A cw 

effective  loop  SNR,  one  can  calculate  the  SATOPS  BER  perfonnance  in  the  presence  of  AWGN 
and  CW  RFI  taking  into  account  of  imperfect  carrier  synchronization  using  (71). 


3.3. 7.3  Analytical  Model  of  PLL  Tracking  Jitter  with  WB  RFI 

Recall  from  the  previous  sections,  the  WB  RFI  signal  is  defined  as: 

/ (t)  =  JUfdft)  cos(2  n{fc  +  A fc)t  +  6l ) 

oo 

d,{t)=  Yjh-P^-nTi  ~t), 

n=— co 

where 

P,  =  RFI  power 

T,  =  RFI  bit  duration  in  seconds,  where  T)  <  Tb 
A fc  =  RFI  frequency  offset  from  the  desired  carrier  frequency  fc 
Oj  =  RFI  phase 

The  variance  of  the  tracking  phase  error,  Oqc,  in  the  presence  of  CW  RFI  can  be  shown  to  have 
the  following  fonn: 

2  _  ^ WB 

6c  LSNR .  (81) 

Here,  one  can  assume  that  fsc  >  4 Rb,  and  the  degradation  factor  AWB  is  found  to  be: 
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where 


(82) 


A„=[\  +  INR.S„(BL,T„Afcj\ 

P 

INR  =  -L- 

SWB  (/,  7) ,  Afc )  =  ^  [SRn  (/  -  A/c  )  +  Sm  (/  +  A/c  )] 


sinC^Ty) 


dj 


here 


TI=l/i?/<Tb=l/^ 


Note  that  for  this  case,  the  effective  loop  SNR,  LSNREff,  is  calculated  from  (81),  using: 
LSNREff  —  Again,  using  this  effective  loop  SNR,  one  can  calculate  the  SATOPS  BER 

1  1  hWB 

perfonnance  in  the  presence  of  AWGN  and  WB  RFI  taking  into  account  of  imperfect  carrier 
synchronization  using  (71). 


3.3. 7.4  Mean  Value  of  Tracking  Error  in  the  Presence  of  RFI 

This  section  derives  a  simulation  model  to  characterize  the  mean  value  of  the  phase  error 
in  the  presence  of  RFI  signal.  The  model  presented  here  assumes  that  the  PLL  is  in  acquisition 
mode  at  Carrier  Modulation  Mode  2  (CMM-2).  This  means  the  transmitted  SATOPS  signal  uses 
direct  modulation  with  a  known  acquisition  sequence,  dACq(t).  For  this  case,  the  received  signal  at 
the  satellite  is  given  by: 

^CMM-2^)  =  ^CMM-2^f)'^d(f)  +  n{f) 

where 

S  CMM-2  (t)  =  ^2Pr  dAcq  (t)  cos(27T(/;.  +  fd  )t  +  ec )  (83) 

where  n(t)  is  the  AWGN  defined  earlier. 


3. 3. 7. 4.1  CW RFI  Case:  Afr  >  B, 

As  defined  previously,  the  CW  RFI  signal  I(t)  is  given  by: 

Hf)  =  V2 Pi  cos(2yr(/c  +  Af)t  +  0, ) 

The  received  signal  becomes: 
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RcMM-2(t)  =  y^PrdAcq(t)cos{27T{fc+fd)t  +  ec )+ cos(2;r(/c  +  tyc)t  +  9I)+  n{t)  ^ 

The  VCO  signal  is  defined  as: 

VCO(t)  =  sin(  2  7v(f,  +fd)t  +  9rco  (t))  ^ 

Assume  that  the  acquisition  sequence  is  perfectly  removed  by  the  correlation  processing  at  the 

satellite  receiver.  The  tracking  error  signal,  e(t),  is  formed  by: 

(0  =  VCO(t).R(t).dAcq (I)  «  J. 1PT  sin(0vco ( t )  -0C)+  ^2P,  sin (0VCO (0 “  dcw 9))-dAcq (0  +  nLPF (1) 
where 

dcw(t)  =  2 K(fd  -  A fc)t  +  Qj 

n,PF  (t)  =  Low  pass  filter  version  of  theAWGN 

The  mean  value  of  the  tracking  error  signal  is  found  to  be: 

E[ecw(t)\~  j2PTE[9vCO(t)-9c] 

Assuming  that  (9vco ( t )  —  9C)  is  small  ^gg^ 

For  CW  RFI,  the  mean  value  of  the  interference  signal  has  been  removed  during  the  sync  word 
correlation  processing  as  shown  in  (87),  and  the  mean  value  of  the  tracking  error  is  the  average 
of  the  difference  between  the  incoming  carrier  phase  and  the  VCO’s  phase. 

3.3.7. 4.2  WB  RFI  Case 

As  defined  previously,  the  WB  RFI  signal  I(t)  is  given  by: 

7(0  =  di  (0  cos(2?r(/c  +  A fc)t  +  9d)  (g9) 

Assume  that  the  acquisition  sequence  is  perfectly  removed  by  the  correlation  processing,  the 
tracking  error  signal  e(t)  is  found  to  be: 


’,U)  =  VCOU).R{t).d  Acq  (t)  ~  sm(6^co  (t)  —  Oc)+4  2  A  d  a«,  0)d,  O)  sin  (Ovco  (/)  -  0CW  (/))  -+-  n,,.h.  O) 


(90) 


Assume  the  acquisition  sequence  and  clAcq(t)  are  independent.  We  get: 

E[ewM^WrE[9vco{t)-9L] 

Assuming  that  (9vco (0  -  9C )  is  small  ^  j \ 

Equation  (91)  shows  that  the  mean  value  of  the  PLL  tracking  error  is  found  to  be  the  same  as  the 
CW  RFI  case. 
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3.3.8  RFI  Detection  Models 


Currently,  there  are  three  approaches  for  the  RFI  detection,  namely: 

•  Approach  #1:  Employ  ITU  and  NTIA  definitions  for  “harmful”  interference  and  use  the 
maximum  allowable  PFD  (Power  Flux  Density)  expressed  in  dB/m  /Hz  to  determine 
whether  harmful  interference  has  occurred. 

•  Approach  #2:  Monitor  the  received  BER  and  when  the  received  BER  falls  below  a 
threshold  value  for  7rh  second,  one  can  declare  an  RFI  event  has  occurred. 

•  Approach  #3 :  Monitor  the  sync  errors,  if  they  exceed  the  threshold  values  for  Trfi 
seconds,  then  one  can  declare  an  RFI  event  has  occurred. 

As  indicated  in  the  introduction  and  our  survey  results  presented  in  Sections  3.1  and  Section  3.2, 
Approach  #1  and  Approach  #2  are  the  ones  that  are  currently  adopted  by  many  existing  RFI 
analysis  tools.  Our  research  effort  considered  a  combined  Approach  #  2  and  Approach  #3  to 
detect  an  RFI  event.  The  combined  approach  takes  into  consideration  of  the  synchronization 
loops  and  BER  performance  of  the  SATOPS  system.  The  combined  approach  includes  the 
following  steps: 

•  Step  1:  Assume  the  PLL  is  in  tracking  mode,  monitor  the  carrier  PLL  tracking  error, 


^Tracking* 

o  RFI  is  flagged  when  the  PLL  carrier  tracking  phase  error  aTrackmg  exceeds  a 
predefined  threshold  value  aThresoid  for  Trfi  second,  i.e.,  aTracking  >  oThreshoid  for 
Trfi  seconds. 

o  From  Section  3.3.6,  the  carrier  PLL  tracking  error,  aTracking,  is  evaluated  using  the 
following  relationships: 


cr 


Tracking 


_  ^RFl  -i 

LSNR^ 


i  =  1, 2 


Recall  that  LSNRq  is  the  loop  SNR,  which  is  defined  as: 


(92) 


LSNR()  =  Jl(m)  Pj 


n0bl 


(93) 


For  i  =  l,  the  loop  receives  the  desired  carrier  signal  component  in  the  presence  of 
CW  RFI,  and  Arfm  is  given  by,  from  Section  3.3.6: 
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Am-,  =[!  +  «*] 

where 

p 

inr  =  — 

Nn 


(94) 


where  INR  is  defined  as  the  Interference  power-to-Noise  Ratio.  For  i  =  2,  the  loop 
receives  the  desired  carrier  signal  component  in  the  presence  of  WB  RFI,  and 
Arfi-2  is  given  by,  from  Section  3.3.7: 


*  Acq-WB 


=  [l  +(INR).S„_„(BDL,f„)\ 


S fu, )  =  [ W  "  fm  >  +  S,  (/  +  /„ )] 


s,  (f)  =  T,  j 


sin  (nfT,) 

xfT, 


df,  where  T,  =1 1  R: 


(95) 


where ^fi  is  defined  as  the  RFI  frequency  relative  to  the  desired  carrier  frequency 

u 

o  From  Section  3.3.6,  Eqn.  (65),  the  value  of  Gxhreshoid  is  calculated  using  the 
following  equation: 

1 


<7r 


Radian 


(96) 


.jLSNREff 

where  LSNRes  is  calculated  from  (70). 

Step  2:  From  Section  3.3.6,  Eq.  (71),  using  the  value  of  GTracking,  evaluate  the  impacts  of 
the  carrier  PLL  tracking  error  on  BER  perfonnance  using  the  following  analytical 
models: 


BERrf,  =  \erfc\— - U  ^ 


1 


® Tracking 


'  Tracking-' 


Tracking  J  v  (^97^ 

where  a  is  defined  as  in  (69).  As  indicated  in  Figure  26,  the  loop-bandwidth-bit-time- 
product  is:  B\Tb  <0.01. 

Step  3 :  Detennine  the  threshold  values  of  the  BER^fi  and  7rfi  caused  by  the  RFI.  The 
threshold  values  of  the  BER^fi  and  7rFi  caused  by  the  RFI  are  defined  as  A£7?Threshoid  and 
F-rhRFi,  respectively. 
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o 


From  Eq.  (97),  the  threshold  value  BERThKsho\d  is  calculated  using  the  following 
equation: 


BER  Threshold  =  \ & fA  - " 

2  G 


Threshold  J 


1  a.(B,Th). 


® Threshold 


& Threshold 


(98) 


where  Gxhreshoid  is  defined  as  in  (96). 

o  The  threshold  of  the  interference  time  duration  Trfi  is  calculated  using  the 
following  relationship: 


Th  RFI  —  fCW,  5  (99 

where  T0  is  the  total  observation  time  within  the  ground  station  antenna  FOV 
(field-of-view),  and  /(%)  is  the  fraction  of  time  that  is  allowed  for  the  SATOPS 
data  to  be  corrupted  by  the  interference  sources.  For  NASA  missions,  the  fraction 
of  time  /(%)  of  1%  is  usually  selected,  and  the  threshold  of  the  interference  time 
duration  Trfi  is  given  by: 


Thitn  -  (0'01).Fo  .  (2Q( 

•  Step  4:  RFI  is  declared  when  BERg$i  and  Trfi  exceed  threshold  values,  i.e.,: 

O  BERrfi  >  5E7? Threshold  when  Glracking  =  Gjhreshold  for  TrFI  ^  F'lliRIT. 

•  Step  5:  In  conjunction  with  the  above  steps  (activities),  the  lock  detector  also  monitors 
the  PLL  tracking  error  Quacking. 

o  The  loop  is  declared  to  drop  lock  when  oxracking  ^  Gxhcar- 

o  Gxhcar  is  the  carrier  threshold  of  the  PLL,  and  it  is  derived  from  (65).  Note  that  in 


practice.  Oxhcar  —  ^Threshold- 

o  If  the  drop-lock  duration  is  greater  than  IxhRFi,  then  evaluate  the  BER. 

■  RFI  is  declared  when  BER  >  RfiRxhreshoid- 

In  practice,  to  ensure  meeting  the  required  command  QoS  and  un-interrupted  PLL  operation,  we 
want  to  set: 


C7 Thcar 


>  G 


Threshold 


(101) 
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Figure  27.  Proposed  Framework  and  Algorithm  for  RFI  Detection 


3.3.9  USB  Signal  Estimation  Models 
3.3.9.1  CW  RFI  Case 

Recall  that  the  SATOPS  command  signal  S(t)  and  RFI  signal  I(t)  received  at  the  satellite 
are  given  by: 


S(t )  =  ciy]2Pr  cos (27u(fc  +  fd  )t  +  md(t)  sin(2 7ifscf)  +  Or)) 
1(f)  =  ^2F)  cos(2 7r(fe  +  A fc )t  +  0j) 

Using  the  identities: 


(102) 

(103) 


oo 

cos  (m  sin(x))  =  J0  ( m )  +  2^  J2n  ( m )  cos(2  nx) 

72=1 


(104) 


oo 

sin(/nsin(x))  =  Z2J2 -  1(m)sin((2/7-l)x) 

«='  ,  (105) 
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the  signal  S(t)  can  be  expanded  as: 

S  (0  =  a^2P^[j0  (m)  cos(2 n(fc  +fd)t  +  Oc)~  2J,  (777)  d  (t)  sin(2  n(fc  +fd)t  +  9C )  sin(2^;,c0  +  h.o.t\ 
h.o.t  =  Higher  order  terms  that  contained  signal  components  at  twice  of  frequencies  or  higher  (106) 

The  goal  here  is  to  estimate  the  “carrier  component”  of  the  USB  signal,  i.e.,  to  estimate  the 
following  signal: 

S(t)  =  aj2P^J0(m)  cos(2tt(/c  +fd)t  +  Gc) 

The  channel  gain  factor,  a,  is  calculated  from  the  estimated  a  (dB)  as: 


(107) 


a(dB ) 

fl  =  10^ 

The  estimated  channel  gain  factor,  a,  expressed  in  dB  is: 

a  a(dB) 

a  =  10^ 


(108) 


(109) 


where 


a(dB)  -  L  (dB)  +  Gr (dB)  +  GR (dB)  +  LAtm  (dB)  +  LCaT(dB)  +  LCaR(dB) 
Ls  (dB)-  Space  loss  =  1 0  log 


'  c  ^ 


v  4Hhfc  j 


4  nA  f2 

Gr(dB)  =  T ransmitter  antenna  gain  in  dB  =  I  Ologl  T  c 


Gr  (dB)  =  Receiver  antenna  gain  in  dB  =  lOlog 


c" 


and 


lAtm  (dB)  =  Loss  due  to  atmospheric  in  dB 
LCaT(dB)  =  Loss  due  to  transmitter' s  cable  in  dB 
LCaR(dB )  =  Loss  due  to  receiver' s  cable  in  dB 

where 

c  =  Speed  of  light  in  m/sec 

D  =  The  distant  from  ground  station  tracking  to  the  satellite  in  meters 
Ax  =  Transmitter  antenna  diameter  in  m2 
Ar  =  Receiver  antenna  diameter  in  m2 

The  Doppler  frequency, /d,  is  estimated  using  the  following  equation: 


(110) 
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fd=fc 


C  +  V. 


c-v„ 


-1 


(111) 


where 

fc  =  Transmitted  carrier  frequency  in  Hz 
c  =  Speed  of  light  in  m/sec 
vs  =  Satellite's  speed  in  m/sec 

vg  =  Ground  station's  sp  eed  in  m/sec.  It  is  assumed  to  be  zero. 

It  should  be  pointed  out  that  the  channel  gain,  a,  can  also  be  derived  from  the  link  budget  model 
described  in  Section  3.3.1 1. 


3.3. 9.2  Estimating  the  Received  Carrier  Phase  6C  -  Case  1:  PLL  in  Tracking  Mode 

For  small  phase  tracking  errors,  Figure  28  illustrated  the  linearized  model  of  the 
advanced  NASA  PLL  shown  in  Figure  18.  As  shown  in  Figure  28,  during  tracking  mode,  the 
phase  error  is  small  and  given  by: 

Oe(t)  =  0c{t)-0vco{t) 

where  the  VCO’s  phase  is  defined  as: 


(112) 


where 


t 

dycoit)  =  2/r.Kvco.jo  J e(r )v(s  -  r)cl r.ds 


pco 

e(t)  =  AK.\o  Oe(r)g(t  -r)dr 


(113) 


(114) 


Figure  29  describes  a  mathematical  framework  for  characterizing  the  estimate  of  carrier  tracking 
phase  error  of  the  advanced  NASA  second-order  PLL  shown  in  Figure  18,  and  how  it  is  used  in 
the  USB  signal  estimation  model.  Figure  30  shows  where  the  carrier  tracking  phase  error  is 
incorporated  in  the  simulation  of  a  PLL  in  the  presence  of  RFI.  Simulation  models  for  the  NASA 
advanced  second-order  PLL  and  typical  Type-I  imperfect  integrator  second-order  PLL  employed 
by  many  SATOPS  systems  will  be  presented  in  the  following  sub-sections.  Both  CW  and  wide¬ 
band  (WB)  RFI  signal  types  will  be  investigated.  For  Type-I  imperfect  integrator,  discrete 
simulation  models  will  be  presented  for  implementation  in  Matlab.  Finally,  a  framework  to  show 


Approved  for  public  release;  distribution  is  unlimited. 


55 


how  these  simulation  models  of  the  PLL  tracking  phase  error  will  be  integrated  into  the  USB 
estimator  model. 


0c(O 


e 


sin(^(0)*^(0 

Oeit) 


AK.B(s)F(S) 


e{t) 


dvco(t)  G(s)  =  AK.B(s)F(s)  ->  g)t) 
V(s)  ->  v(t) 


KvcoK(s)V(S) 


PLL  Closed  loop  transfer  function: 

B(s)F(s)K(s)V(s ) 


Hpll(s)  = 


1  +  B(s)F(s)K(s)V(s) 


Figure  28.  Linearized  Model  of  Advanced  NASA  PLL  Shown  in  Figure  18  for  Tracking 

Mode 


Oe(t)  =  0c(t)-Qvco(t) 


t 


6,  =  6, 


YM)+^vC0(t)  =  ^A(t) 

dt  dt  dt 


7  00  7 

0,(0  +  2x-Kvcoj  v(r)e(r  -z)dz  =  —  0(0 


t 


dt 


vco 

Phase  of  the  PLL  VCO 
Provides  an  estimate  of 
the  received  carrier 
phase 


0pco(  0  —  J*  £?(t)v(v  z)dT.ds 


e(t)  =  AK.[d, 


d_ 

dt 


co  |- 

de{t)  +  2n.AK.KVC0.\v(T)  £  0e(s)g(t  - s)ds 

1 


LAlm  ( dB )  =  Loss  due  to  atmospheric  in  dB 
LCaT  {dB)  —  Loss  due  to  transmitter's  cable  in  dB 
LCaR{dB)  —  Loss  due  to  receiver's  cable  in  dB 
c  =  Speed  of  light  in  m/sec 

D  =  The  distant  from  ground  station  tracking  to  the  satellite  in  meters 
At  —  Transmitter  antenna  diameter  in  m2 
Ar  =  Receiver  antenna  diameter  in  m2 
fc  =  Transmitted  carrier  frequency  in  Hz 
c  =  Speed  of  light  in  m/sec 
vs  =  Satellite's  speed  in  m/sec 

v  =  Ground  station' s  speed  in  m/sec.  This  can  assumed  to  be  zero 


Tracking 

Mode 


PLL 


A 


0c  —  0|  'co 


input 


USB  Signal 
Estimation 
Model 


Output 
S(t) 

S(t )  =  a  j2r\J0(m)  cos(2  *(/c  +fd)t  +  0c) 


Figure  29.  Mathematical  Framework  and  Algorithm  to  Incorporate  PLL  Tracking  Error 

into  USB  Signal  Estimation  Model 
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S(t) 


0,0 


7(0  ”(0 

44- 


0,(0 

* 

<?(0 


^(s)  =  B(s)F(j)K(j) 


AK.B(s)F(s)V(s) 


r(?) 


^FCo(  0 


KvcoK(s) 


Figure  30.  Simulation  Model  for  Incorporating  Carrier  Tracking  Error  into  NASA 

Advanced  Second-Order  PLL 


3.3. 9.3  Mathematical  Model  of  NASA  Advanced  Second-Order  PLL  Tracking  Error 
3. 3.9. 3.1  Case  h  CWRFI 

Recall  from  previous  sections,  the  signals  S(t),  I(t),  n(t)  and  r(t)  shown  in  Figure  30,  are 
defined  as  follow: 


where, 


and 


Sit)  =  a  j2FrJ 0  (m)  cos(2;r  (fc  +fd)t  +  0c) 

I {t)  =  V273  cos(2/r(/;;  +  Afc)t  +  O, ) 

n(t)  =  AWGN  with  two-  sided  powerspectral density  of  N0  /  2 


(115) 

(116) 


r{t)  =  V2  sin(2/z/;/  +  0vco(t)) 
Let’s  define  the  Heavy-side  operator  “p”  as: 

d  M 
P  =  —  (•) 
dt 


(117) 


(118) 


For  typical  NASA,  civil  and  commercial  SATOPS  receiver,  the  loop-gain  and  filters  shown  in 
Figure  30  are  given  below: 

AK  —  loop  gain  =  2.4  x  107 ,  (119) 

(,20> 
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(121) 


trc  —  1.6  x  10  5  . 

For  a  Type-I  imperfect  integrator  loop  filter,  F(s)  is  given  by: 


Ks)  = 


1+  TjS 
1+  TiS 


where 


Ti  =  4707,  t2  =  0.0442 
T y  —  1.0  X  10-6  . 

The  phase  smoother  filter,  V(s),  is  defined  as: 


v<*>  = 


The  VCO  filter  is: 


K(s)  =  -  . 


We  define  the  following  system  parameters: 


/SR  =  F\ 


Jl(m)PT 


K  =  AK.Kvco 

dc  =  9C  (O  =  co0t  -e0=27f0t-  0{) 

0,  (0  =  2nA fct  +  2  rfj_osct  +  0,  =  2  rfjt  +  0 J 
where 

fl  =  Afc+fl-OSC 

Am0  =27Z-(/7+/o) 


AG,  =6,-0, 


(122) 


(123) 

(124) 


(125) 


(126) 


(127) 

(128) 
(129) 


(130) 

(131) 

(132) 


Oe(i)  =  Oc(t)  -Ovco(t)  ( 

The  carrier  tracking  phase  error  of  a  typical  second-order  PLL  in  the  presence  of  RFI  can  be 
characterized  by  a  differential  equation  shown  below: 


(133) 


- - Q  (/)  = - ~Fa(p)  sin(6>e (0)  +  VM sin(0  (?)  +  Amot  +  A0o)+ - ^-=  n(f) 

KJ,(ni)JPT  dt  KJ0(m )JFT  [ 

For  a  typical  advanced  NASA  Type-I  second-order  PLL,  FA(p)  is  given  by: 


.  (134) 
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fAp)  = 
fAp)  = 
fAp)  = 
fAp)  = 


1  +  z2p 

i 

1 

1 

isT 

+ 

_ i 

_1  trcP  _ 

A+tvP_ 

1  +  r2p 


(l  +  TJ)f\  +  T 

RC  p\l  +  TvP) 

\  +  t2p 


\  + T\P  +  TrcP  +  TrCT{P  +Tvp  +  T{Uvp  +  TRCTvp  +TrcT1TvP 

1  +  r2p 

1  +  (O  +  Trc  +  Tv)p  +  ( TrcTx  +  T{CV  +  TRCTv)p2  +  TRCTlTvpi 
1  +  r2p 

1  +  Tclp  +  TC2P  +  tC2p 

where 

^C1  =  T\  +  F RC  +  O 

TC2  =  GcO  +  00  +  GcO 
=  GcOO 


(135) 


(136) 


Thus,  from  (135),  a  differential  equation  characterizes  the  behavior  of  the  PPL  tracking  phase 
error  in  the  presence  of  narrowband  RFI  and  AWGN.  It  can  be  shown  to  have  the  following 
form: 

vr,  \  rrr  <3  +  TciP  +  TciP2  +  TcyP' )  1  A(0  = 

KJ0(m)^]PT  dt 


(1  +  rclp  +  rC2p2  +  rc3p3)cj0 


-(t  +  r2p) 


sin(  0e  (?))  +  V I  SR  sin(  0(,  (t)  +  A  m0t  +  A6?0)  4 - -  . — 

J0(m)  yjPT 


KJ0(m)-JPr 

One  can  rewrite  the  above  equation  and  obtain  the  following  equation: 

d  Q  . ,  d2  _  . ,  c/3  _  . ,  o'4  _  .  ’ 

(/?  dt  dt  dt 


77(0 


(137) 


KJ0(m)Jp 


1 


sin(£(,  (0)  +  a/ /Si?  sin(6(  (?)  +  A  i u0t  +  A0O)+  , — 

J0(m)^Pr 


77(0 


—  sin(#e (t ))  +  V /Si?  —  sin(6( (?)  +  A m0t  +  A 90 )  H - ^  , —  — 

dt  dt  J0(myPT  dt 


«(0 


(138) 


33.9.3.2  Case  2;  Wideband  (WB)  RFI 

For  WB  RFI,  recall  that  the  interference  signal  has  the  following  fonn: 

7(0  =  a/T0^/  (0  cos(2  jr(fc  +Afc)t  +  0j)  (|3  9) 

From  (137),  the  differential  equation  characterizes  the  behavior  of  the  PPL  tracking  phase  error 
in  the  presence  of  wideband  RFI  and  AWGN.  It  can  be  shown  to  have  the  following  form: 


Approved  for  public  release;  distribution  is  unlimited. 


59 


- — ==  —6(t)+Tcx^d(t)+TC2^e(t)  +  TC3^e(t)  = 

KJ0(m)^[dt  eW  Cldt2  e  C2dt2  e  C3dt4  cWJ 

-  sin(6((f))  +  V ISRdj ( t ) sin($, ( t )  +  Am 0t  +  A0O) - ^  , —  n(t) 

d0(myPT 

-  *jt  sin(^ (0)  +  yfdSRd, ( t ) sin(0e (0  +  Am 0t  +  AGQ)  +  y  ( ~ft  n ^ 


(140) 


33.9.4  Discrete  Simulation  Model  of  Type  I  Imperfect  Integrator  PLL  Tracking  Error 

This  sub-section  presents  the  simulation  models  for  a  second-order  PLL  using  Type-I 
imperfect  integrator.  Figure  31  depicts  a  block  diagram  of  this  PLL. 


Figure  31.  Simulation  Model  for  Incorporating  Carrier  Tracking  Error  into  Typical  Type-I 

Imperfect  Integrator  Second-Order  PLL 


33.9.4.1  Case  1:  CWRFI 


For  Type-I  imperfect  integrator  loop  filter  with  the  worst  case  CW  RFI,  Eq.  (134) 


4r 0e (0  =  —  -  sin(0,(f))  -  flSR  sin (0,(0)  -  ^ 


becomes: 


-  —  fl  +  K'  r,  cos(0  (0)  +  K'  t2  flSR  cos(/9  (?))1—  0e  (/)  +  —  —  nt  (f ) 
r,  "  dt  A't1  dt 

where  we  have  defined  the  worst  case  scenario  when  6 ft)  =  0,  A aT0  —  0,  and  A 9  —  0.  The 
parameters  A\  K ’  and  interference  to  signal  ratio  (I SR)  are  defined  as: 

A'=  J0(m)^2Pr 


(141) 


(142) 
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K'_k_  ak.kvco 

A'  J0(m)42PT 


(143) 


p 

ISR  =  — 


Note  that  m(t)  is  the  defined  as: 


nx(t)  =  nc(t)cos(0e(t))  -  ns(t)sm(6e(t)) 


(144) 


(145) 


where  nc(t)  and  ns(t)  are  the  in-phase  and  quadrature  components  of  the  AWGN,  n(t).  Viterbi  [7] 
has  shown  that  nj(t)  is  essentially  white  Gaussian  noise  with  spectral  density  No/2  when  the 
carrier  loop  bandwidth  is  smaller  than  the  input  noise  bandwidth.  Using,  Eqn.  (138)  can  be 
transformed  into  two  first-order  differential  equations  that  are  suitable  for  simulation  using 
Matlab  program.  If  we  define  Ts  as  the  sampling  interval  and  “/?  ”  as  the  number  of  sample  points, 
then  (138)  can  be  rewritten  as: 

9Jin  +  \)Ts)  =  9JnTs)  +  Tsd'jnTs) 

0e((n  +  1)7))  =  —  -  sin(0e  (nTJ)  -  sin (9e  ( nTs ))  - 

i JiViv 

-  —  \  +  K't2  cos(0e(nTs))  +  K\2yfMms(0e(nTs))^(nTs)  +  -^—nlN(nTs) 
h  h  SNR  (146) 


where 


e'M)  =  -ee{t) 

at 


(147) 


Oe(nTs)  = 


0e((n+\)Ts)-de((n)Ts) 


>h  N(.nTs)  = 


n\(nTs) 

n\  (nTs)\ 


When  07(t)  =  9I0,  A <n0  =  0,  and  A0O  =  60,  Eqn.  (146)  becomes: 


(148) 


(149) 
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0e((n+l)Ts)  =  0e(nTs)+Ts0'e(nTs) 

e'e{{n+\)Ts)=T-^ 


-sin {Oe{nTs)  +  0a)-yflSR  sm(0e  (nTs)+0IO)- 


SNR 


-TN^  +  K'T2cos(0e(nTs)+0o)+K'T24lSRcos(0e(nTs)+eio)\p'e(nTs)+-^—n1N(nTs) 
r,  t1SNR 


(150) 


33.9.4.2  Case  2:  WB  RFI 

For  Type-I  imperfect  integrator  loop  filter  with  the  worst  case  WB  (wideband)  RFI,  Eqn. 
(134)  becomes: 

^0Xt)  =  -\-sHm)-^dI(t)sm(0Xt))-^ 

dr  r,  A 


1 


|l  +  K't2  cos(0  (t))  +  K'r2 -JlSRdjit) cos {0 At))  —0 At) +  — — n, (t) 

dt  .IT  di 


(151) 


Similarly,  this  equation  can  be  transformed  into  two  first-order  differential  equations  that  are 
suitable  for  simulation  using  Matlab  program.  It  is  found  to  be: 

0e((n  +  l)Ts)  =  0e(nTs)  +  Ts0'e(nTs) 

0'e((n  +  l)Ts)  =  7^-  -sin {0AnTs))-^dI{nTs)Sm(nTs))-'hj^ 


I  +  K't. 


cos(0e(nTs))  +  K32yfmdI(nTs)cos(0e(nTs))^'e(nTs)  +  ^—nlN(nTs) 


t{SNR 


(152) 


When  0/(t)  =  9I0,  An>0  =  0,  and  A0O  =  0O,  Eqn.  (152)  becomes: 

9e((n  +  l)Ts)  =  0e(nTs)  +  Ts0'e(nTs) 

^((n  +  l)r,)  =  ^[-sin(^(«r5)  +  0J-VMd/(^)sin(^(«rs)  +  0/o)-^^ 

a/V/v 


M  L 

^\  +  K32Cos(03nTs)  +  0o)  +  K32^dI(nT3CoS(0A,lTs)  +  0I3}):(nT3  +  ^-nlN(nTs) 

T'  t>SNR  .  (153) 


The  characteristics  of  the  second-order  Type-I  imperfect  integrator  employed  by  a  typical  NASA 
PLL  is  shown  in  Table  4.  The  loop  parameters  xi  and  i2  in  (145)  and  (149)  are  calculated  based 
the  PLL  tracking  loop  bandwidth  selected.  Table  4  shows  the  values  for  the  two-sided  loop  noise 
bandwidth,  2BL,  two-sided  pre-detection  noise  bandwidth,  22?Lpd,  normalized  gain  constant^’, 
loop  parameters  xi  and  X2. 
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Table  4.  Typical  Type-I  Imperfect  Integrator  Employed  By  NASA  PLL 


Tracking 

Mode 

2Bl 

(Hz) 

2f?LPD 

(Hz) 

r 

(Sec) 

(Sec) 

Narrow 

1.0 

200 

1340.9 

33940.0 

1.50 

3.0 

200 

1340.9 

6510.0 

0.50 

10.0 

2000 

1340.7 

339.4 

0.15 

30.0 

2000 

1340.7 

65.1 

0.05 

Wide 

10.0 

2000 

13409.3 

3394.0 

0.15 

30.0 

2000 

13409.3 

651.0 

0.05 

100.0 

20000 

13406.7 

33.9 

0.015 

300.3 

20000 

13406.7 

6.5 

0.005 

3.3. 9.5  Lock  Point  and  Definition  of  Loss  of  Lock  of  Type-I  Imperfect  Integrator  PLL 

It  is  assumed  that  the  PLL  is  initially  in  lock  with  the  phase  of  the  desired  signal,  and  the 
RFI  signal  is  subsequently  injected  into  the  receiver.  Since  the  Type-I  imperfect  integrator  is  of  a 
lag-lead  type  filter,  the  PLL  is  a  second-order  loop,  and  consequently  only  one  stable  point  in  the 
phase  error  region  of  (-n,  +n).  In  this  same  region,  there  also  exists  a  single  saddle  point.  Thus, 
using  the  same  argument  in,  the  selected  threshold  for  detennining  the  loss  of  lock  of  a  PLL  is 
Till.  This  threshold  will  be  used  in  analyzing  the  loss  of  lock  of  a  PLL  using  Type-I  imperfect 
integrator. 

Eqn.  (150)  is  simulated  in  Matlab  and  the  results  are  presented  in  Figure  32,  Figure  33, 
and  Figure  34.  These  figures  show  the  simulation  of  the  second-order  Type-I  imperfect  integrator 
PLL  with  SNR  =16  dB,  ISR  =  -15  dB,  and  initial  phase  offset  6q  between  the  desired  signal  and 
the  reference  VCO  phase  is  set  at  tc/4  radians.  Figure  32  shows  the  tracking  jitter  settled  down  to 
7t/4  as  expected.  Figure  33  shows  the  difference  of  tracking  error  fluctuates  around  0  radians,  and 
Figure  34  shows  the  PLL’s  VCO  phase  starts  at  0  degrees  and  it  eventually  locks  onto  the 
incoming  phase  at  7t/4  radians.  When  the  initial  phase  offset,  between  the  desired  signal  and 
the  reference  VCO  phase  is  set  at  0  rad,  and  the  phase  offset  between  the  RFI  signal  and  the 
reference  VCO  phase  is  set  at  tc/  1 8  rad,  Figure  35  shows  that  the  PLL  drops  lock  on  the  desired 
signal  phase  at  0  degree,  and  it  eventually  locks  on  to  the  incoming  RFI  phase  at  tc/  1 8  radians. 
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Tracking  Phase  Error  in  Radians 


Time,  seconds 


Figure  32.  Tracking  Phase  Error  as  a  Function  of  Time 


Figure  33.  Difference  of  Tracking  Error  as  a  function  of  Time 
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Figure  34.  Phase  Plane  Plot  for  a  Typical  PLL  Using  Type-I  Imperfect  Integrator  for  ISR 

-15  dB 


Tracking  Phase  Error  in  Radians 


Figure  35.  Phase  Plane  Plot  for  a  Typical  PLL  Using  Type-I  Imperfect  Integrator  for  ISR 

25  dB 
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33.9.6  USB  Carrier  Signal  Component  Estimator 

Figure  36  shows  how  to  incorporate  the  simulation  models  for  the  carrier  phase  error, 
described  by  equations  (137),  (139),  (145),  and  (149),  into  the  USB  carrier  estimator  model  to 
generate  an  estimate  of  the  received  SATOPS  command  signal,  S(t). 


Figure  36.  Block  Diagram  of  the  USB  Carrier  Signal  Component  Estimator 


3.3.10  RFI  Signal  Estimation  and  Prediction  Models 

3.3.10.1  Case  1:  CW  RFI  Signal  Estimation  and  Prediction  Model 

This  section  describes  the  CW  RFI  prediction  models  using  the  USB  carrier  signal 
component  estimator  presented  previously.  Figure  37  presents  a  high-level  block  diagram  for 
extracting  the  CW  RFI  signal  from  the  received  signal,  R(t).  Recall  the  received  signal  at  the 
SATOPS  satellite  is: 

R{t)  =  S(t)  +  I(t)  +  n{f) 

S(t)  =  a  j2P^J0(m)  cos(2  n(Jc  +  fd)t  +  Oc)  ( 
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I(t)  =  JlP,  cos(2*r(/c  +  A/; )t  +  0, )  5  (155) 

where 

n(t)  =  AWGN  with  two-  sided  power  spectral  density  of  N0  /  2 

In  Figure  37,  Rj(t)  is  defined  as  the  estimate  of  the  CW  RFI  signal  at  RF  frequency  in  the 
presence  of  noise  and  error  signal  caused  by  inaccurate  estimate  of  the  USB  carrier  signal,  i.e., 

R,  (t)  =  I(t )  +  n(t )  +  s(t) 

where  s(t)  is  the  error  signal  caused  by  in  accurate  carrier  signal  estimator  due  to  RFI 

s(t)  =  S ( t )  -  S (t)  =  a  yf2P^J0 (m) cos(2 n{Jc  +  fd )t  +  0c)~  a -J 2P. J 0 (m) cos(2 n(fc  +fd)t  +  0c)  _  ^  1 5 6) 

Due  to  the  RFI,  the  PLL  may  not  provide  an  accurate  estimate  of  the  incoming  carrier  phase 
signal.  However,  in  practice,  one  knows  the  transmitted  power,  and  can  estimate  the  Doppler 
frequency  and  channel  gain  factor  fairly  accurate,  thus  without  loss  of  generality  one  can  set: 


a  =  a  =  \^’=fd=fd 


(157) 


The  error  signal  becomes: 


£(t)  =  ^2PRJJm) 


cos(2  n{fc  +  fd)t  +  6C )  -  cos(2  n(Jc  +  fd)t  +  6C) 


(158) 


Rewriting  the  error  signal,  one  gets: 


A  O+O 

sm(2 n(fc+fd)t  +  (  c2  c))sm 

(  A  > 

0-0 

c  c 

2 

\2PRJ0(m) 

.V  ) 

(159) 


Down-converting  the  received  signal,  Ri(t),  signal  to  baseband  signal,  Ribb(0,  one  obtains: 


RlBB  (0  — 


R, (t).2cos(2x:(  fc  +  fd)t  =  JlPj  cos(2 TiAfjt  +  9,)  -  ^2PR  J0(m)[siii(6»;)siii(6»;)]+  N'(t) 


,  (160) 


where 

Afd=Afc-fd 
0+  = 


0-  = 


2 

0-6, 


(161) 
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and 


N'(t)  =  j2N0.N(t) 


(162) 


where  N(t)  is  the  nonnalized  AWGN  with  zero-mean  and  unit  variance.  For  small  tracking  error, 
one  can  approximate  (160)  as  follows: 


ribb (0  ~  V2 ~p,  cos(2 nAfjt  +  6,)- 


2PRJ0(m) 


■&c(t)  +  N\t) 


(163) 


where,  for  small  PLL  carrier  phase  tracking  jitter,  one  can  write: 

Qc(t)  =  [sin(6>;)sin(6>;)] ~  0e  0).sin(6> );  0C  <  U- 

2  .  (164) 

Here,  one  has  assumed  that  6C  is  a  constant  phase  offset  between  the  desired  signal  and  the  PLL 
VCO’s  phase,  and  that  the  PLL  phase  tracking  error  is  small,  such  that  9e  —  8C  —  9C  «  0  and 
9C  *  8C. 


USB  Signal 
Estimation 
Model 


S(/)  =  a^lP^J„(m)cos{ln(fI+fd)t  +  ec) 


.  Ri(J) 

m — >(V) - ^ 

\  T  t 


I  R1BB  (0 


42  cos(2»(/c  +/,)/) 


Low  Pass  Filter 
Bandwidth  =  B, 


CW  RFI 
Estimation, 
Prediction 
and  Detection 
Model 


I(t)  =  V 2 P,  cos(2 7i(fc  +  A/, )t  +  0,) 


R(t)  =  5(0  +  /(0  +  «(0 


Figure  37.  Block  Diagram  for  Extracting  the  CW  RFI  Signal  from  the  Received  Signal,  R(t) 


Figure  38  presents  a  proposed  approach  for  CW  RFI  estimation,  prediction  and  detection 
model.  The  proposed  approach  estimates  the  CW  RFI  signal  power  first  and  uses  its  estimate  for 
estimating  the  CW  RFI  frequency.  The  estimated  CW  RFI  frequency,  along  with  the  estimated 
CW  RFI  power,  is  then  used  for  the  CW  RFI  phase  estimation  and  detection.  The  details  of  these 
steps  will  be  described  in  the  following  sections. 
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CW  RFI  Estimation,  Prediction  and  Detection  Model 


Reset  B, , 
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Ribb  (n) 


t  ^Tsamp 

Tsamp  =  Sampling  time 
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CW  RFI 
Detector 


9, 


CW  RFI 
Phase 
Estimator 


^IBB  (0 


CW  RFI  Signal 


IIP,  cos(2 n(fc+  Afc )f  +  6t ) 


Figure  38.  Block  Diagram  for  CW  RFI  Estimation,  Prediction  and  Detection  Model 


Before  getting  into  the  detailed  discussion  of  the  CW  RFI  estimation,  prediction  and 
detection  models,  it  is  noted  that  from  (160),  one  observes  that  if  the  CW  RFI  frequency  can  be 
estimated  accurately,  then  the  3-dB  cutoff  bandwidth  of  the  LPF  is  set  exactly  at  A/j  such  that  the 
CW  RFI  signal  is  passed  through  it  without  distortion.  If  the  3-dB  cutoff  bandwidth  of  the  LPF  is 
not  known,  the  LPF’s  bandwidth  is  set  larger  than  the  RFI  frequency,  A/i.  Then,  more  noise  is 
allowed  to  pass  through  the  filter,  causing  performance  degradation  in  the  “CW  RFI  estimation, 
prediction  and  detection”  model.  The  proposed  approach  to  resolve  this  problem  is  described  as 
follow: 

•  Step  1:  Set  the  initial  LPF’s  bandwidth  BWlpf,  =  BL  (PLL  tracking  loop  bandwidth).  Note 
that  in  order  to  cause  potential  interference  to  SATOPS  receiver,  the  RFI  frequency  A/J, 
must  be  less  than  the  PLL  tracking  loop  bandwidth. 

•  Step  2:  Perform  the  CW  RFI  signal  power  and  frequency  estimations  to  obtain  the  initial 
estimate  of  the  CW  RFI  signal  power  and  frequency. 

•  Step  3 :  Use  the  estimate  of  the  CW  RFI  signal  frequency  to  reset  LPF’s  bandwidth 
BWlpf- 
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•  Step  4:  One  can  repeat  Steps  2  and  3  until  the  estimate  of  the  CW  RFI  frequency  A/i, 
settles  to  a  constant  frequency  and  uses  that  frequency  to  set  the  LPF’s  bandwidth  . 

3.3.10.1.1  Estimating  the  CW  RFI  Signal  Power 

A  combined  CW  RFI  signal  power  estimator  and  the  RFI  signal  power-to-noise  power 
ratio  estimator  (INR)  is  shown  in  Figure  39.  The  proposed  approach  shown  in  Figure  39  has 
assumed  the  signal  properties  do  not  change  during  the  observation  time.  It  calculates  the 
average  power  of  all  points  received  since  the  last  reset.  When  signal  properties  change  during 
the  observation  time,  the  sliding-block  discussed  in  Ref.  [8]  would  be  used  to  estimate  the  CW 
RFI  signal  power. 


Figure  39.  Block  Diagram  for  CW  RFI  Signal  Power  Estimation 


If  one  defines  the  carrier  signal  power-to-noise  power  spectral  density  ratio,  CSNR, 
interference  signal  power-to-carrier  signal  power  ratio,  ICSR,  and  the  normalized  baseband 
signal  Rum(t),  Rmm(t)  as,  respectively: 


CSNR  =  Pr 

NoBlpf  ,  (165) 

where  the  7?Lpf  is  initially  set  at  BL.  As  described  above,  the  optimum  setting  of  RLpf  happens 
when  it  approaches  the  true  value  of  the  CW  RFI  signal  frequency: 

ICSR  =  — ^ — 

PRJl(m)  (166) 
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^IBBN  6) 


R IBB  (0 


\2  PRJ0(m) 

Eqn.  (163)  can  be  rewritten  in  terms  of  CSNR  and  ICSR  as  follows: 


^IBB  (t)<zyj2PRJ0(m)\ 


V/GW.cos(2M/;f  +  6,)  - (jj.0c(f)  + 


m 


J2CSNR 


(167) 


(168) 


or 


Ribbn  (t )  -  JICSR. cos(2^Af ,t  +  9,)-  [|]-@e(0  +  ^== 

Note  that  in  (165),  PLpf  is  defined  as  the  bandwidth  of  the  LPF  that  is  dependent  on  the  RFI  CW 
signal  frequency.  As  mentioned  above,  the  initial  LPF  bandwidth  will  be  set  at  the  PLL  tracking 
bandwidth  B\  . 

The  unknown  power  Pi  of  the  CW  RFI  can  be  estimated  by  using  power  estimator  or  by 
taking  the  expected  value  of  the  square  of  the  baseband  signal  Pibbn(0,  (168).  From  Figure  39, 
the  output  of  the  power  estimator  is  found  to  be: 


i{Ram(t)Y]*ICSR.E[coS(2xAfIt  +  0I)]2  +  (±W(oF  + 


Here 


E[cos(2^Afjt  +  6,  )]2  =  i 

E[®c(t)]2=al  |sm(0)f 


and 


J  2  CSNR  _ 

(169) 

2  CSNR 

(170) 

(171) 

(172) 

(173) 

p[^vfi)]2=i 

Note  that  we  have  assumed  that  0C  is  a  constant  phase  offset  between  the  desired  signal  and  the 
PLL  VCO’s  phase.  Hence  sin{0c )  is  a  constant  between  0  and  1.  Similarly,  from  Figure  39,  the 
output  of  the  variance  estimator  is  found  to  be: 


Va>%RmN  (t).cos(2rfHt)f  \PF  ] »  ^ jcr;  {sin(0  )f  + 

The  estimate  of  the  CW  RFI  power  is  given  by,  from  Figure  39: 
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p,= 


2  PRJ\m) 


\e\Ribhn  (0)2]-  Var\Rmw  (t).cos(27ifHt)f\ 


=  P 

,PF  J  I 


(175) 


As  shown  in  Figure  39,  the  sampling  frequency  is  Fsamp  and  the  observation  time  T0 b  in  seconds 
is  defined  as: 


T  =N  T  = 

±Ob  ly  oh 1  Samp 


N 


ob 


Fc 


(176) 

where  Na b  is  the  number  of  samples  during  the  observation  time.  The  perfonnance  of  the  CW 
RFI  signal  power  described  in  (175)  is  given  by: 


Var\  PI  = 


2 


N„ 


Wcr2  {sm(0  )f 


+  - 


1 


2CSNR 


[nobl] 


N„ 


(177) 


Note  that  the  factor  is  calculated  using  (79). 


3.3.10.1.2  Estimating  the  CW  RFI  Signal  Frequency 


Since  the  CW  RFI  frequency  is  unknown,  the  proposed  approach  here  is  to  estimate  the 
RFI  frequency  before  detection.  Using  (168),  the  nonnalized  input  signal  of  R\\m(t)  to  the  “CW 
RFI  Frequency  Estimator”  shown  in  Figure  38  can  be  approximated  as,  taking  into  account  of  the 
CW  RFI  power  estimate: 


R„ 


\2P,  ■  cos(2 7^f,t  +  9 ,)  +  NmN  (0 


where  the  frequency,  A /J,  and  phase,  0\,  of  the  CW  RFI  are  unknown,  and 


(178) 


,oc<- 

c  2 


(179) 


Mbbn(0  can  be  characterized  as  AWGN  with  zero-mean  and  variance  is  given  by: 


Vat  \N IBBN  (t)]  —  cIBBN  —2 PR  «/q  (m)\ 
Note  that  the  CSNR  for  this  case  is:  CSNR 


-sin2(0f  ).cr 

4  y  c  J  ee 

PrJo  (m) 


1 


N°BLPF  'Blpf=Bl 

interference  signal-to-noise  power  ratio  (ISNR)  is  given  by: 


2  CSNR 

_  PrIoQu ) 
N0BL 


;  0C  <- 
2 


(180) 


.  Therefore,  the 
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I  SNR  = 


P, 


NA 


p, 

1  +  A<:'v  sin2(i9c) 


(181) 


This  subsection  presents  an  approach  to  estimate  the  CW  RFI  frequency,  A /J,  at  low  ISNR 
using  iterative  weighted  phase  averaging  (iWPA)  [9].  From  Figure  39,  the  sampled  signal 
Rbbn(«)  is  given  by: 


Pbbn  («)  =  V2  P,  .cos(2  7f0n  +  0j)  +  Nam  (n) 


(182) 


The  corresponding  unknown  analog  frequency  A/J  that  is  being  estimated  can  be  obtained  from 
the  following  relationship: 


Samp 


(183) 


where  F^.dmv  is  the  sampling  frequency  that  is  defined  as  in  (176).  The  discrete  frequency  fo  is 
being  estimated  and  can  be  written  as: 


u  = 


(I  +  A)_ 


N, 


ob 


N. 


ob 


A<- 
2  2 


(184) 


where  N0 b  is  related  to  the  observation  time,  Tob,  which  is  defined  as  in  (176),  L  is  the  integer 
part,  and  A  is  the  fractional  part  of  Xq,  where  X(}  denotes  the  number  of  signal  cycles  contained  in 
the  observation  window.  To  understand  the  above  notations,  suppose  that  one  observes  a  CW 
RFI  signal  of  frequency  32.5  Hz,  i.e.,  A/J  =  32.5  Hz,  that  is  sampled  at  128  Hz  during  1  second 
(N=  128  samples).  Then ,fo  =  32.5/128.  In  this  case,  X0  =  32.5  (the  frequency  cycles),  L  =  32, 
and  A  =  0.5. 

Using  the  above  notations,  the  Discrete  Fourier  Transform  (DFT),  X(k),  of  Rmis(n)  at 
spectral  line  k  is  given  by,  assuming  the  noise  term  is  negligible: 


X(k)  = 


2  P, 


LejeiW\ 


Uk-k 

vv  Nob  j) 


ff 


w 


k  +  A0 


J) 


(185) 


where  W(f)  is  the  Fourier  transfonn  of  the  selected  time-domain  window.  For  instance, 
considering  a  rectangular  window  of  length  N,  we  have: 
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W(  f)  =  sirKW^X)  e-j!f(Nob-X) 

sin  (4) 


(186) 


Since  the  frequency  A/i  of  the  CW  RFI  is  far  from  the  origin, /=  0,  such  that  the  leakage  coming 
from  the  negative  part  of  the  spectrum  can  be  neglected,  i.e., 


X(k)  =  ±^eJd'W\ 


Uk-A 

VV  Not  jj 


(187) 


The  coarse  estimate  of  the  CW  RFI  frequency  is  performed  by  selecting  its  largest  peak,  i.e., 


rc, 

/o 


-Max\ 


gdw(C)|) 


N„ 


(188) 


The  fine  estimate  of  the  CW  RFI  frequency  is  perfonned  by  estimating  the  fractional  part  A  of  Xo 
using  iWPA  [9].  The  DFT  signal  X(k)  is  divided  into  M non-overlapping  segments  of  length  P: 

<>)  =  <>  +  ^)’°^n<P-l  (189) 

The  spectrum  of  each  segment  is  evaluated  at  the  previously  estimated /o: 


P- 1 


Wo)=iy „(«)« 


-j2xf0(n+lP) 


n= 0 


(190) 


Substituting  (182),  assuming  that  the  noise  term  is  negligible  and  rectangular  window,  and  taking 
into  account  that  f0  —  (Nob  —  A)  / Nob  ,  one  can  show  that: 


Wo)  = 


2  P, 


sin  (nAP/Nob) 
sin (nA/Nob) 


J0,  e-j(2  API  Ncbt\+l) 


(191) 


The  phase  of  Xt  (/0)  is  given  by: 


arg|  Xt  (/0 )  \  -6i-  (l  +  /) 


N 


ob 


(192) 


The  frequency  error  can  be  estimated  by: 


A  N 
A  =  ob 


2  nL 


arg 


A 


Wo) 


-arg 


A 


X2(f«) 


v  7  v  7J.  (193) 

In  the  above  equation,  we  have  split  the  original  register  into  two  segments  ( P  =  Nt. 2).  The  fine 
frequency  estimator  is  obtained  using  the  following  relationship: 


Approved  for  public  release;  distribution  is  unlimited. 


74 


(194) 


A 


r  fine 
'  0 


L  +  A 


TV 


The  block  diagram  for  CW  RFI  frequency  estimation  is  shown  in  Figure  40.  The  frequency 
estimator  described  above  approaches  the  Cramer-Rao  lower  bound  at  high  SNRs,  which  is  the 
lowest  variance  that  can  be  attained  by  any  unbiased  estimator.  At  high  ISNR,  the  variance  of  the 
above  frequency  estimator  can  be  approximated  as: 


var 


(  A  ) 

f  fine 

Jo 

« 

V  / 

12  [N0Bl] 


2^P/AoA(7Vo6--1)_ 


1  +  ^ 


[sin(0  )]2 


(195) 


P, 


^  IBB  ) 


Figure  40.  Block  Diagram  for  CW  RFI  Frequency  Estimation 
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3.3.10.1.3  Estimating  the  CW RFI  Signal  Phase  and  Detecting  the  CW RFI  Signal 


This  subsection  describes  a  Maximum  Likelihood  CW  RFI  phase  estimator  and  CW  RFI 
signal  detector  based  on  the  approach  presented  in  [10].  Figure  41  illustrates  the  proposed 
approach  by  incorporating  the  CW  RFI  power  and  frequency  estimates.  From  Figure  41,  the 
objective  of  the  CW  RFI  phase  estimator  is  to  estimate  the  phase  of  the  following  CW  RFI  signal 
by  incorporating  the  power  and  frequency  estimates: 


RB(t) 


2P,  .cos(2;r  A f,  t  +  d,)  +  Nam (t) 


(196) 


where  the  phase,  0\,  of  the  CW  RFI  are  unknown,  and  the  noise  term  Mebn^is  defined  as  in 
(179).  It  can  be  characterized  as  AWGN  with  zero-mean  and  variance  as  defined  as  in  (180). 
Note  that  (180)  can  be  rewritten  as: 


O' mm,  =  Var[Nmm  (1)]  =  [N0Bl  ]  1  +  sin2  (dc ) 


K<^ 


(197) 


The  CW  RFI  phase  estimator  described  in  Figure  41  can  be  implemented  in  digital  domain  using 
a  sampling  frequency  of  Fsamp  and  an  observation  time  Fob  in  seconds.  The  Fob  is  defined  as  in 
(176).  For  this  case,  since  the  CW  RFI  frequency  must  be  within  the  PLL  tracking  loop 
bandwidth  BL,  the  sampling  frequency  must  be  sampled  at  the  Nyquist  rate,  i.e., 


^ Sump  2.Bl 


N 

rp  _  V  ob 

1  Ob  ~ 


2  Bt 


(198) 


(199) 


Again,  N0 b  is  the  number  of  samples  during  the  observation  time. 

At  high  ISNR,  the  performance  of  the  CW  RFI  phase  estimator  shown  in  Figure  41 
approaches  Cramer-Rao  lower  bound,  which  can  be  shown  to  have  the  following  form: 

6  '  ^  [Un2Noh+\2nNob(Nob -l)  +  2Nob(Nob  -l)(2Aoi  -1)] 


var 


0, 

v  J 


'  IBBN  | 


PjNob\Nob2- 1) 


(200) 


where  no  is  the  time  to  at  which  the  first  sample  is  taken  for  the  estimate  of  the  CW  RFI  phase. 
Here  to  is  defined  as: 


■  n  T  =  -l^L 

"O1  Samp 


(201) 
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Figure  41.  Block  Diagram  for  CW  RFI  Phase  Estimation  and  RFI  Detection 


From  [10],  the  false  alarm  probability  as  a  function  of  the  threshold  detection  value,  y,  can  be 
shown  to  have  the  following  fonn: 


PF=e 


7 


2a 


2 

IBBN 


Solving  (198)  for  the  threshold  value: 

Y  =  ] 


where  ofBBN  is  given  by  (197).  Similarly,  the  probability  of  detection  as  a  function  of  the 
threshold  detection  value,  y,  can  be  shown  to  have  the  following  fonn  [10]: 


Pd=Q 


ISNR. 


<T„ 


(202) 

(203) 


(204) 


Approved  for  public  release;  distribution  is  unlimited. 


77 


where  ISNR  is  defined  as  in  (181),  and  Q(a,  b )  is  the  Marcum  Q-function  defined  by: 

- 

Q(a,b)=j  ze-  2  l0(az)dz 


(205) 

By  combining  (198)  and  (199),  one  obtains  the  receiver  operating  characteristic  (ROC)  given  by: 

PD(r)  =  Q{lSNR^j-2\n(PF))  (206) 


Plot  of  Eq.  (202)  is  shown  in  Figure  42  with  the  detection  interference  SNR  (ISNR)  as  a 
parameter.  For  a  given  false  alarm,  the  threshold  detection  value  y  is  chosen,  using  (203),  to 
achieve  the  desired  probability  of  detection  at  a  given  ISNR. 


Figure  42.  Receiver  Operating  Curve  for  the  CW  RFI  Signal  Detector  Shown  in  Figure  41 


3.3.10.2  Case  2:  WB  RFI  Signal  Estimation  and  Prediction  Model 

This  section  describes  the  WB  RFI  prediction  models  using  the  USB  carrier  signal 
component  estimator  presented  in  the  previous  section.  Similar  to  the  CW  RFI  case  shown  in 
Figure  37,  Figure  43  shows  a  high-level  block  diagram  for  extracting  the  WB  RFI  signal  from 
the  received  signal,  R(t). 
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Vwb  (0  =  V2  pi  di  (0  cos(2 7T(  fc  +  A/;  )t  +  9j ) 


Rivb  0)  =  S it)  +  IWB  it)  +  nit) 

Figure  43.  Block  Diagram  for  Extracting  the  WB  RFI  Signal  from  the  Received  Signal, 

R  unit) 

The  received  signal  at  the  SATOPS  satellite  in  the  presence  of  WB  RFI  is  given  by: 

Rjvb (0  =  n(t)  ^0 

where  S(t)  and  n(t)  are  defined  as  before,  and  the  WB  RFI  is  defined  as: 


hvR  (0  =  yf2 ~Pldl  ( t )  cos(2 n{fc  +  Afc  )l  +  O, ) 


(208) 


Again,  in  Figure  43,  RwyvJt)  is  defined  as  the  estimate  of  the  WB  RFI  signal  at  RF  frequency  in 
the  presence  of  noise  and  error  signal  caused  by  inaccurate  estimate  of  the  USB  carrier  signal, 


R-IWB^f)  —  Iwb( 0  +  «(0  e(t) 


(209) 


where  s(t)  is  the  error  signal  caused  by  inaccurate  carrier  signal  estimator  due  to  RFI,  which  is 
defined  in  (159).  Down-convert  the  received  signal  R\wB(t)  signal  to  baseband  signal  RiwBB(t), 
one  obtains: 


RlWBB  (0  R IWB  (t).2cos(2n(fc  +  fd)t\ 


LPF 


(210) 


Similarly,  for  small  tracking  errors,  one  can  approximate  (209)  as  follows: 


v 


y 


(211) 


where  Afh  6c(t)  and  N’(t)  are  defined  as  in  (161),  (162)  and  (164). 
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Figure  44  shows  an  approach  for  WB  RFI  estimation,  prediction  and  detection  modeling. 
The  proposed  approach  estimates  the  WB  RFI  signal  power  first  and  uses  its  estimate  for 
estimating  the  WB  RFI  data  rate  (hence  bandwidth  of  the  WB  RFI  signal).  These  estimates  are 
then  used  for  the  estimation  of  WB  RFI  frequency.  The  estimated  WB  RFI  frequency  along  with 
the  other  estimates  is  then  used  for  the  CW  RFI  phase  estimation  and  detection. 

Similar  to  CW  RFI  case,  the  3-dB  cutoff  bandwidth  of  the  low  pass  filter  (LPF)  must  set 
according  to  the  RFI  signal  bandwidth  (or  RFI  data  rate  Ri )  such  that  the  WB  RFI  signal  can  pass 
through  without  distortion.  Since  the  3-dB  cutoff  bandwidth  of  the  LPF  is  not  known  at  the 
beginning,  the  following  steps  are  proposed  to  address  this  issue: 

•  Step  1:  Set  the  initial  LPF’s  bandwidth  BWlpf,  =  Brcx,  where  Rrcx  is  the  SATOPS 
receiver  bandwidth.  This  bandwidth  is  set  based  on  the  99%  power  containment 
bandwidth.  As  an  example,  for  uncoded  bit  rate  of  2  kbps,  the  minimum  subcarrier 
frequency  recommended  by  CCSDS  is  8  kHz.  For  this  case,  the  subcarrier  frequency-to- 
bit  rate  ratio  is  4.  The  CCSDS  Recommendation  4.2.2  recommends  that  the  99% 
power  containment  bandwidth  for  this  case  is  16f?b  for  modulation  index  “m”  from  0.8 
radians  to  1 .2  radians.  For  this  case,  Rrcx  is  set  at  1 6/C  and  that  the  CW  RFI  bandwidth 
can  be  assumed  to  be  less  than  16Rb  for  uncoded  system.  For  coded  system,  it  is  16RS, 
where  Rs  is  the  symbol  rate. 

•  Step  2:  Perfonn  the  WB  RFI  signal  power  and  WB  RFI  data  date  estimations  to  obtain 
the  initial  estimate  of  the  WB  RFI  signal  power  and  bandwidth. 

•  Step  3 :  Use  the  estimate  of  the  WB  RFI  signal  bandwidth  to  reset  LPF’s  bandwidth 
B  Rlpf- 

•  Step  4:  One  can  repeat  Steps  2  and  3  until  the  estimate  of  the  WB  RFI  data  rate,  R\, 
settles  to  a  constant  data  rate  and  uses  that  data  rate  to  set  the  LPF’s  bandwidth  BWlpf- 
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Initial  Bandwidth 
Setting: 

Blpf  =  SATOPS 

Receiver 

Bandwidth 


WB  RFI  Estimation,  Prediction  and  Detection  Model 

R-IWBB  (0 


WB  RFI  Signal 


2  Pt  dj  ( t )  cos(2/r(  fc  +  Afc  )t  +  6, ) 


Figure  44.  Block  Diagram  for  the  WB  RFI  Estimation,  Prediction  and  Detection  Model 


3.3.10.2.1  Estimating  the  WB  RFI  Signal  Power 


A  combined  WB  RFI  signal  power  estimator  and  the  WB  RFI  signal  power-to-noise 
power  ratio  estimator  (INR)  is  shown  in  Figure  45.  Similar  to  the  CW  RFI  case,  the  proposed 
approach  shown  in  Figure  45  has  assumed  the  signal  properties  do  not  change  during  the 
observation  time.  It  calculates  the  average  power  of  all  points  received  since  the  last  reset.  When 
signal  properties  change  during  the  observation  time,  the  sliding-block  discussed  in  Ref.  [8]  is 
used  to  estimate  the  WB  RFI  signal  power. 

The  carrier  signal  power-to-noise  power  spectral  density  ratio,  CSNR,  interference  signal 
power-to-carrier  signal  power  ratio,  ICSR,  and  the  normalized  baseband  signal  Rmm(t), 

^iwbbn (t)  as,  respectively: 


CSNRwb 


PrJIQk) 

N0BWlpf 

? 


where  the  BWLPF  is  initially  set  at  Brcx-  The  optimum  setting  of  BWlpf  happens  when  it 
approaches  the  true  value  of  the  WB  RFI  signal  bandwidth: 


(212) 
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and 


ICSRwb  = 


Prvbbn  (0 


p, 

prJ  <>0) 

PflVBB  O 

■^2PRJ0(m) 


R-IWBBN  (0  *  4 1 CSRnB.d , (0cos(2M//  +  e,)~ Q j.0t ._WB (0  + 

? 

where  0C  wsCO  denotes  the  tracking  phase  error  in  the  presence  of  wideband  RFI. 


(213) 

(214) 


(215) 


Figure  45.  Block  Diagram  for  the  WB  RFI  Power  Estimation 

Note  that  in  (212),  BWlpf  is  defined  as  the  bandwidth  of  the  LPF  that  is  dependent  on  the 
RFI  CW  signal  bandwidth.  As  mentioned  above,  the  initial  LPF  bandwidth  will  be  set  at  the 
SATOPS  receiver  bandwidth,  BRCx-  Note  that  BRCX  is  initially  set  at  16Rb  for  uncoded  SATOPS 
systems.  For  this  case,  the  CW  RFI  bandwidth  is  assumed  to  be  less  than  1 6/C. 

The  unknown  power  Pi  of  the  WB  RFI  can  be  estimated  by  using  power  estimator  or  by 
taking  the  expected  value  of  the  square  of  the  baseband  signal  Riwbbn(0,  Eqn.  (215).  From 
Figure  45,  the  output  of  the  power  estimator  is  found  to  be: 
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and 


Arisen  0)f  ]*  E  J ICSRWB  .[cos(2M//f  +  e, )]-^00c_TO (*)]+ 


(O)2 ]»  /CW(),.E[^,(f)coS(2M//f  +  0, )]2  +  fi W_ra(of  + 


(216) 


(217) 


Here, 


E[d  t  ( t )  cos(2  nbfjt  +  0,  )]2  =  ^ 

4®«(0F  =  <4,  _,,-[sin(0.  )f 

4w(op=i 

? 

where  <JwB_ee  is  calculated  using  (81).  Note  that  one  has  assumed  that  0C  is  a  constant  phase 
offset  between  the  desired  signal  and  the  PLL  VCO’s  phase.  Hence,  sin(Qc)  is  a  constant  between 


(218) 

(219) 

(220) 


0  and  1.  Similarly,  from  Figure  45,  the  output  of  the  variance  estimator  is  found  to  be: 


Var%R,WBBN  (t).  coi(2  jf Ht)f  \PF  ]  = 


-  \.<7u, 


.[sin(6»c  i]2 


1 


2CSNRn 


The  estimate  of  the  CW  RFI  power  is  given  by,  from  Figure  45: 


P,= 


2Pr  J2(m) 


' 0))  (f).cos(2;zfwf))  ^ 


P 

,PF  1  1 


(221) 


(222) 


From  Figure  45,  when  the  observation  time  is  /V0b7samP,  the  performance  of  the  WB  RFI  signal 
power  estimator  described  in  (222)  can  be  shown  to  be: 


Var  Pj  = 


2PRJHm) 


N 


.[sin(^)]2 


+  - 


1 

2CSNR,, 


[n0bl] 


N„ 


1  +  A'";  [sin(0  )f 


.  (223) 


Note  that  the  factor,  AWb,  is  calculated  using  (82). 


3.3.10.2.2  Estimating  the  WB  RFI  Signal  Frequency 

Similar  to  CW  RFI  case,  the  WB  RFI  frequency  is  unknown  and  the  proposed  approach 
is  the  identical  the  CW  RFI  case  with  a  slight  modification.  Figure  46  shows  a  block  diagram  for 
the  WB  RFI  signal  frequency  estimator.  The  input  signal  Riwbb(«)  to  the  “CW  RFI  Frequency 
Estimator”  shown  in  Figure  38  can  be  approximated  as: 
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R-IWBBN  («)  ~  V2 Pi  -di  (0 C0S(2 7f0n  +  6,)  +  Nwbbn (n) 

where  the  frequency,  A/J,  d\(n)  and  phase,  0i,  of  the  CW  RFI  are  unknown,  and 


N  ib bn  (0  —  pPRJ0(m). 


From  Figure  45,  the  minimum  sampling  time,  rSamp,  at  Nyquist  must  be  set  at: 

T  = _ 1 _ =  1 

Samp  2(16  Rb)  32  Rh 

where  Mwbisn^  is  AWGN  with  zero-mean  and  variance  is  given  by: 


^ar\N IWBBN  (O]  —  < 


'  IWBBN  (m) 


1  . 


sin  (Oc  )-(JWB  + 


1 


2CSNR,, 


=  '•<! 


The  CSNRwb  for  this  case  is: 


CSNRm  = 


PRJ;{m) 


PRJ;{m) 


N0BWlpf  'BW^=m6R^  16 N0Rb 


(224) 


(225) 


(226) 


(227) 


Therefore,  for  WB  RFI  case,  the  Interference  Signal-to-Noise  Power  Ratio  (ISNR)  becomes: 


ISNRura  = 


WB  2 

°  IWBBN  16  MR 


l  +  ^p-sin2(<9c) 


(228) 


Again,  A0b  is  the  number  of  samples  selected  during  the  observation  time  Toh  (see  Eqn.  (176)). 
The  factor  AWb  is  calculated  using  (82). 
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p, 


Figure  46.  Block  Diagram  for  the  WB  RFI  Signal  Frequency  Estimator 


If  one  ignores  the  noise  term  in  (224)  for  a  moment,  and  process  only  the  WB  RFI  signal 
using  the  proposed  frequency  estimator  approach  shown  in  Figure  46,  one  can  show  that: 


Rwbbn  («)  =  V 2  P; .  cos(2;r(2  f0  )n  +  6, ) 


(229) 


The  corresponding  unknown  analog  frequency  A/i  that  is  being  estimated  can  be  obtained  from 
the  following  relationship: 

=  f 

yr  SampJ  0 

Samp  ,  (230) 


2A f,=- 


where  Fsamp  is  the  sampling  frequency  that  is  defined  as  in  (226).  The  discrete  frequency /o  is 
being  estimated  can  be  written  as: 


2/a=2(^)=2^° 


N„, 


1  .  1 

— <A<  — 
2  2 


(231) 
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Similar  to  CW  RFI  case,  the  value  of  N0 b  is  related  to  the  observation  time,  T0 b,  which  is  defined 
as  in  (176),  L  is  the  integer  part,  A  is  the  fractional  part  of  Ao,  and  Ao  denotes  the  number  of 
signal  cycles  contained  in  the  observation  window.  One  can  go  through  the  remaining  signal 
processing  steps  shown  in  Figure  45  and  show  that  the  coarse  estimate  of  the  WB  RFI  frequency 
is  perfonned  by  selecting  its  largest  peak,  i.e., 


fo 


coarse  WB 


-Maxi 


arg 


(kwi)' 


N, 


N. 


(232) 


Similarly,  the  fine  estimate  of  the  WB  RFI  frequency  is  performed  by  estimating  the  fractional 
part  AWb  of  A().  It  can  be  shown  to  have  the  following  form: 


A  wb  =2. 


N. 


ob 


2  tzL 


arg 


Wo) 


arg 


Wo) 


=  2  A 


(233) 


From  Figure  45,  the  fine  frequency  estimator  can  be  shown  to  be  identical  to  the  case  for  CW 
RFI,  i.e., 


r  fine_WB 

/  0 


_  L  _  A 
2 - +  2- 


N„ 


Nn, 


L  +  A 


N, 


ob 


(234) 


By  substituting  (234)  into  (230),  one  obtains  the  estimated  WB  RFI  frequency  as  follow: 


A  f  -  F  fflneJWB  _  r? 

Lt>  ,  —  1  Samp  Jo  ~  1  Samp- 

The  performance  of  the  proposed  WB  RFI  frequency  estimator  shown  in  Figure  46  is  expected  to 
be  worse  than  the  CW  RFI  frequency  estimator  due  to  the  increase  of  the  noise  power  caused  by 
the  squaring  processing  of  the  signal. 


L  +  A 


N 


ob 


3.3.10.2.3  Estimating  the  WB  RFI  Signal  Phase  and  Detecting  the  WB  RFI  Signal 

This  subsection  modifies  the  Maximum  Likelihood  CW  RFI  phase  estimator  and  CW 
RFI  signal  detector  presented  in  the  previous  section  for  CW  RFI  case.  Figure  47  presents  the 
modified  approach.  From  Figure  46,  the  signal  objective  of  the  CW  RFI  phase  estimator  is  to 
estimate  the  phase  of  the  following  CW  RFI  signal  by  incorporating  the  power  and  frequency 
estimates: 


Approved  for  public  release;  distribution  is  unlimited. 


86 


(236) 


Rwbb  (0  *  i 2  P, .  d,  (t)  dj  (t)  cos(2;r  Af  t  +  e, )  +  Nwb  (0 
where  the  phase,  8h  of  the  CW  RFI  are  unknown  and  the  noise  tenn  Nweft)  is  defined  as: 

A 

^wr  (0  =  di(t)  NmBN  (t)  (237) 

The  noise  tenn  expressed  in  (236)  can  be  characterized  as  AWGN  with  zero-mean  and  variance 
is  defined  as  in  (197).  If  one  assumes  that  the  estimated  wideband  data  d/(t)  is  synchronized 
with  d\(t),  then  one  can  approximate  the  product  d/(t)d/(t)  as: 

(t)  ( t )  «  dj  ( t )  =  1 

Using  the  above  approximation,  one  can  write  (235)  as: 


(238) 


rwbb  (f)  ~  V2  P1  ■  cos(2;r  Aflt  +  0,)  +  NWB  (t) 


(239) 


Using  (238),  one  can  show  that  the  signal  Rb(0  is  identical  to  CW  RFI  case.  The  perfonnance  of 
the  proposed  WB  RFI  phase  estimator  and  WB  RFI  detector  shown  in  Figure  47  is  expected  to 
be  worse  than  the  CW  RFI  frequency  estimator  due  to  the  increase  of  the  noise  power  caused  by 
the  timing  synchronization  error  between  the  incoming  WB  RFI  and  the  estimated  WB  RFI  of 
the  signal.  When  the  timing  synchronization  error  is  negligible,  one  can  characterize  the 
perfonnance  of  the  phase  estimator  and  signal  detector  using  the  CW  RFI  analytical  models 
presented  above. 
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Figure  47.  Block  Diagram  for  WB  RFI  Phase  Estimation  and  RFI  Detection 


3.3.10.2.4  Estimating  the  WB  RFI  d,(t)  and  Data  Rate 


Figure  48  describes  the  proposed  estimator  for  estimating  the  WB  RFI  data  stream  d\(t) 
and  data  rate  R\.  From  the  figure,  the  signal  Rd/(t)  is  found  to  be: 

R-d,  (0  =  PM  (f)-R-IWBB  (0 


Rd[  (t)  -2PId1  ( t ) cos(2;r  A/',  t  +  6j) cos(2^  Af)  t  +  0l)  +  Noise T erm 


(240) 


Equation  (240)  has  assumed  that  Pj  «  Pt.  If  one  assumes  that  A/)  «  Af,  and  9,  ~  6,  then  (240) 
becomes: 


Rd/  ( t )  =  Pl  d, ( t)+  P,  dj ( t ) cos(4;r  AfIt  +  20I)  +  Noise T erm 


(241) 
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From  Figure  48,  the  output  of  the  LPF  is  found  to  be: 

Rdj  (?)  =  Pj  dj  (?)  +  Noise  T erm 

The  input  to  the  hard-limiter  shown  in  Figure  48  becomes: 

_  .  .  ...  Noise  Term 

RdN(t)  =  d,(t)  + - ^ - 

P, 

The  noise-tenn  in  (143)  is  characterized  by  (225).  It  can  be  shown  that: 

[sin(<9  )f 


(242) 


(243) 


Var 


Noise  Term 

_  \.N0BrCX. 

l  P,  J 

p,  L 

1  +  -A"»l 


(244) 


where  BRCx  is  the  SATOPS’  receiver  bandwidth,  and  it  is  initially  set  at  1 6/?h  for  uncoded 
SATOPS  systems.  From  Figure  48,  the  output  of  the  integrator  is  RdlN(Tj )  given  by,  assuming 
NRZ  data  format  with  estimated  bit  duration,  7),  for  the  WB  RFI  data  stream: 


Rdmm= ±\+ni(ti) 

where  one  can  show  that  1^(7})  is  an  AWGN  with  zero-mean  and  variance: 

fall 


(245) 


(J  = 

Ni 


p  T  L 
1 1 1 1 


i+^ 


[sin(^)]2 


(246) 


Again,  the  factor  AWb  is  calculated  using  (82). 
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A/; 


Figure  48.  Block  Diagram  for  WB  RFI  Data  Stream  and  Rata  Rate  Estimation 


From  Figure  48,  the  estimate  of  the  WB  RFI  data  stream  d/(t)  is  derived  from  Eqn.  (245), 


i.e., 


dI(t)  = 


RdJTI)>0^+l 
JL  (£)<()-> -1 


(247) 


The  error  of  the  estimate  is  characterized  by  the  bit  error  rate  (BER)  of  the  estimated  WB  RFI 
data  stream.  From  (245),  Rd;jv(7))  is  a  Gaussian  random  variable  with  zero-mean  and  variance  is 
given  by  (246).  One  can  show  that  the  BER  for  the  estimated  WB  RFI  data  stream,  BERwb  rfi, 
has  the  following  fonn  [11]: 


BER,, 


=  Q 


d2. 

min 

I 4.cr2 

ni  y 


(248) 


where  cr^;  is  the  variance  of  the  noise  IV/ ( 7)),  and  dmm  is  the  minimum  Euclidian  distance 
between  the  binary  bits  +1  and  -1,  which  is  2.  Substituting  dmm  =  2  into  (248),  one  gets: 
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BER„ 


=  Q\ 


v%  J 


(249) 


where  Q(x)  is  defined  as: 


Q(x )  =  J —  |  e  2  dt  =  —  erfc{x) 


^ Jv  2  .  (250) 

As  shown  in  Figure  48,  the  integration  time  of  the  integrator  is  controlled  by  the  WB  RFI  timing 
synchronizer,  and  the  timing  synchronization  error  can  cause  potential  degradation  in  the 
estimate  of  the  WB  RFI  data  stream.  If  one  assumes  that  the  digital  transition  tracking  loop 
(DTTL),  shown  in  Figure  49,  is  used  to  track  the  WB  RFI  timing  clock,  then  the  WB  RFI  timing 
clock  error  for  a  first-order  loop  (i.e.,  F(s)  =1)  is  found  to  have  the  following  fonn  [12]: 
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iA+-L 

4  2  E 


f 


2p, 


1-^ 

v  4 


(251) 


where 


s  e )  11 
A  =  — — x — —  =  Normalize  timing  error; —  <  A  <  — 

rri  2  2 


C  =  DTTL  window  width;  0  <  c  <1 
E,  =PITI  =  WB  RFI  Bit  Energy 


(252) 

(253) 

(254) 


Pi  = 


N  R 


=  DTTLLoop  SNR;  BDTTL  =  DTTLloop  bandwidth  in  Hz 


The  timing  clock  jitter  is  given  by: 


cr  =TjJa- 


(255) 


(256) 


where  is  given  by  (25 1)  and  7)  is  the  estimate  of  the  WB  RFI  bit  duration  that  is  generated  by 
the  WB  RFI  data  rate  estimator.  As  shown  in  Figure  48,  the  data  rate  estimator  estimates  the  WB 
RFI  data  rate,  and  hence  the  WB  RFI  bit  duration  7),  from  the  power  spectral  density  (PSD) 
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estimator.  Following  is  the  discussion  on  the  WB  RFI  PSD  and  data  rate  estimators  (and  hence 
WB  RFI  bit  duration  7)). 


Tj 


T, 


Figure  49.  Block  Diagram  for  WB  RFI  Synchronizer 


Figure  50  illustrates  the  block  diagram  for  the  WB  RFI  PSD  and  data  rate  estimators. 
Similar  to  previous  sections,  the  signal  RdlN(t)  is  sampled  at  sampling  Tsam p  and  the  observation 
time  rob  is  defined  by  (176).  Auto-correlation  function,  rdjdj(k)  ,  is  computed  using  the  sampled 
signal  RdlN(n).  The  input  signal  to  the  DFT  (Discrete  Fourier  Transform)  block,  shown  in  Figure 
50,  is  calculated  using  the  following  relationship: 


r 

drd, 


d, 

J2+\ 


(*)  =  IX  {n)Rdm  {n  -  k);  k  =  0,  ±  1,  ±  2,  ±  3,  ...,±^ 


(257) 


The  PSD  function  of  the  signal,  RdlN(n),  is  computed  by  taking  the  DFT  of  the  auto-correlation 
function  given  by  (257),  and  it  is  given  by: 
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PSD(f)  \  =  PSD(k)  =  DFT{d  di  (*)} 


Tob 

Not- 1 


PSD{k)=  X  rdidi(n)e 


N„ 


;k  =  0+l,±2,...±^ 


(258) 


For  NRZ  data  format,  the  discrete  PSD  of  a  WB  RFI  is  found  to  be: 


PSD(k)  =  -k- 


srn 


irkT, 


nkT, 


N 

•  ■  k  =  0  +1  +2  +— 2k 

,  ,  IV  - J.  , - ^ 


(259) 


The  plot  of  (259)  is  shown  on  the  bottom  left  comer  of  Figure  50.  As  shown  in  Figure  50,  the 
spectral  spacing  of  PSD(k)  or  frequency  resolution  is  defined  in  terms  of  the  observation  time  in 
sec,  Fob,  as: 
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Figure  50.  Block  Diagram  for  WB  RFI  PSD  and  Data 
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It  is  important  to  note  that  the  sampling  time  rsamp  must  be  fast  enough  to  accurately 
capture  the  WB  RFI  spectrum.  For  the  case  of  interest,  the  CW  RFI  bandwidth  is  assumed  to  be 
less  than  16f?b-  Thus,  the  minimum  sampling  frequency,  Fs,  must  be  at  the  Nyquist  rate,  i.e., 

FSmp=32Rob.  (261) 


From  (261),  the  required  number  of  samples  of  PSD(k)  or  the  observation  data  point,  /V„b,  within 
the  observation  time,  rob,  is  found  to  be: 


N0b  ~  ^’samp^ob  ~  T°'  ~  32RbTob 

Sa'"p  .  (262) 

The  WB  RFI  data  rate  can  be  estimated  using  the  nulls  (zero-points)  of  the  PSD  curve.  Since  the 

WB  RFI  signal  is  corrupted  by  the  additive  noise,  sampling  noise  and  quantization  noise,  the 

zero-points  might  be  shifted  and  raised  above  zero.  The  proposed  algorithm  to  estimate  the  WB 

RFI  data  rate  for  this  case  is  illustrated  in  Figure  50.  The  algorithm  searches  for  the  minimum 

points,  say  within  the  observation  time  T0b.  There  would  be  Nmm  points.  If  one  assumes  that  the 

7th  zero-point  is  shifted  by  £,  and  (/-l)th  zero-point  is  shifted  by  £j_i,  then  the  distance  between  the 

7th  minimum  point  and  the  0'-l)th  minimum  point  is  calculated  using  the  following  equation: 


«0  =  ^-t^M.=0,±l,±2 . 


(263) 


Note  that  in  the  absence  of  noise,  i.e.,  a,  =  Si_i  =  0,  and  (263)  represents  a  true  estimate  of  the  WB 
RFI  data  rate.  Due  to  the  presence  of  noise,  8;  ^  Sj_i  ^  0.  The  WB  RFI  data  rate  can  be  estimated 
using  the  following  equation: 


R'  (W„  -1)  M 


i  ou-n 

£*/(»') 


The  estimate  of  WB  RFI  data  bit  duration  is  given  by: 

T  _  1 
1  A 

Rr  0 


(264) 


(265) 


As  mentioned  earlier,  the  timing  synchronization  error  can  cause  potential  degradation  in  the 
estimate  of  the  WB  RFI  data  stream.  It  is  of  interest  to  assess  the  impacts  of  the  timing  error  on 
the  estimate  of  WB  RFI  data  stream.  The  effect  of  timing  error  on  the  BER  perfonnance  of  the 
estimate  of  WB  RFI  signal  can  be  modeled  by  observing  the  following  cases: 
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•  Case  1:  If  two  successive  (rectangular)  bits  are  identical,  an  incorrect  bit  reference  will 
have  no  effect  on  the  bit  error  probability;  and 

•  Case  2:  If  two  successive  bits  differ,  the  magnitude  of  the  expected  correlator  output  is 
reduced  by  a  factor:  (1  —  2\A\),  where  X  is  the  nonnalize  timing  error,  i.e., 


A  =  - 


T  iming  Error 
Data  Bit  Duration 


(  A-) 

£-£ 

V _ )_ 

T: 


(266) 


Based  on  the  above  observation,  the  conditional  probability  of  WB  RFI  bit  error  in  the  presence 
of  timing  error  can  be  shown  to  have  the  form: 


PWB  x{E!X)  =  \-Q 


1 


[1_2lAl  1 


(267) 


The  Probability  Density  Function  (PDF)  of  the  WB  RFI  timing  error  X  in  the  presence  of  AWGN 
for  the  proposed  DTTF  synchronizer  can  be  characterized  by  a  Tikhonov  distribution,  namely, 


p  _  exP  {Pwb-rfi  cos(2^L)}  2^1  <  y 

^o(PwB-RFl) 


(268) 


where 


PwB- 


RFI 


1 


cr; 


'*  .  (269) 

Here,  er^  is  given  by  (25 1).  The  average  BER  for  the  estimate  of  the  WB  RFI  data  stream  in  the 
presence  of  imperfect  timing  recovery  becomes: 


BER 


WB-RFI  _  Ti  min  g 


//>„_,(£/ UP, 


(270) 


where  Pwb~.a(E/A)  and  (A)  are  given  by  (267)  and  (268),  respectively. 


3.3.11  USB  Link  Budget  Models  for  Assessing  SATOPS  Performance 

This  section  provides  a  detailed  description  for  the  USB  link  budget  model,  shown  in 
Figure  9,  for  estimating  the  received  signal  power  PR  at  the  satellite  receiver  and  assessing 
impacts  of  RFI  on  the  USB  command  system  perfonnance.  The  link  budget  model  described 
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below  is  for  the  uplink  command  SATOPS  models.  These  models  consist  of  two  key 
components: 

•  Link  margin  model  for  carrier  tracking,  and 

•  Link  margin  model  for  command  demodulation. 


3.3.11.1  Link  Margin  Model  for  Carrier  Tracking 


The  link  margin  model,  ZMcar rier,  for  the  carrier  tracking  is  given  by: 

C ' 


LM  ^  (dB)  =  EIRP(dBm)  +  Gr  (dBi)  -  -  (dB) 

VNjreqd 

~KTfTc  (dBm)  -  Ls  (dB)  -  L0  (dB)  -  MLc  (dB)  (271) 

where  EIRP  is  the  transmitted  signal  power  expressed  in  Effective  Isotropic  Radiated  Power 
(EIRP)  expressed  in  decibel-milli-Watts  (dBm);  antenna  gain,  Gr,  is  in  decibels  referenced  to 
isotropic  gain  (dBi);  (C/N)reqd  is  the  required  carrier-to-noise  power  for  a  specified  tracking  loop 
bandwidth,  BL ;  k TsWc,  the  maximum  thermal  noise  power,  is  in  decibel-miliwatts  (dBm)  and  is 
given  in  terms  of  the  Boltzmann’s  constant,  k,  the  system  temperature,  Ts,  and  the  carrier  loop 
bandwidth,  Wc;  Ls(dB)  is  the  space  loss  in  dB;  L(}  is  the  circuit  loss  plus  weather  loss  plus 
scintillation  loss  plus  degradation  due  to  RFI  plus  other  losses;  and  MLC( dB)  is  the  carrier 
modulation  loss  expressed  in  decibels. 


3.3.11.2  Link  Margin  Model  for  Command  Demodulation 


The  link  margin  model,  LMcm d,  for  the  command  demodulation  is  given  by: 


LMcmd  (dB)  =  EIRP(  dBm)  +  Gr  (dBi)  - 


<  E  ^ 

ft 

V^o  Jreqd 


(dB)  -  At dB-bits/s) 


-kTs  (dBm/Hz)  -  L  (dB )  -  L0  ( dB)  -  ML.  (dB) 


(272) 


where  all  the  parameters  are  defined  as  before  except  (Eb/No)reqd  is  the  required  bit  SNR  for  a 
specified  command  BER  (bit  error  rate);  MLcm( dB)  is  the  command  modulation  loss  expressed  in 
decibels;  and 

R(dB  -  bits  /  sec)  =  10  log[/?/;  ]  (273) 

Recall  that  Rb  is  the  command  bit  rate  in  bits/second. 
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For  both  link  budget  models,  the  space  loss,  Ts(dB),  can  be  calculated  using  (110),  and  the  EIRP 
can  be  calculated  using  the  following  equation: 


EIRP(dB)  =  PT  ( dB )  +  Gt  (dB)  -  Feed  _  Loss(dB ) 


(274) 


where  7\(dB)  is  the  transmitted  power  in  dB,  FeedLoss(dB)  is  the  antenna  feeder  loss  in  dB, 
and  G'  KdB)  is  the  transmitted  antenna  gain  in  dB,  and  it  is  given  by: 


GT(dB)  =  101og(?7rv)+  201ogf  nDl^'' 


(275) 


V  c  j 

where  r\ Tx  is  the  transmitted  antenna  efficiency,  DTx  is  the  transmitted  antenna  diameter,  and  fc  is 
the  transmitted  carrier  frequency,  and  c  is  the  speed  of  light.  In  the  following  subsections,  the 
weather  loss,  scintillation  loss,  degradation  due  to  RFI,  and  the  carrier  modulation  loss  will  be 
discussed  in  detail. 


3.3.11.3  Weather  Loss  in  USB  SATOPS  Systems 

The  weather  effect  at  S-band  is  almost  negligible.  Specifically,  the  attenuation  due  to  dry 
atmosphere  and  water  vapor  is  almost  0  dB.  In  addition,  Table  5  shows  the  rain  attenuation 
computations  for  2  GHz  and  3  GHz  at  various  values  of  annual  probability.  The  values  in  the 
table  were  computed  assuming  that  the  polarization  is  circular,  elevation  angle  is  15°,  the 
satellite  orbit  is  180°  W,  the  city  is  Tampa,  Florida,  which  has  the  highest  rain  precipitation  in 
the  U.S.,  and  the  Crane  rain  region  is  E. 
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Table  5.  Rain  Attenuation  Loss  in  dB  Calculated  for  2GHz  and  3  GHz 


Annual  Percent 

Rain  Attenuation  (dB) 

Rain  Attenuation  (dB) 

Of  Time 

Frequency  =  2  GHz 

Frequency  =  3  GHz 

0.1 

0.02 

0.04 

0.8 

0.02 

0.05 

0.6 

0.02 

0.06 

0.5 

0.02 

0.07 

0.4 

0.03 

0.07 

3.3.11.4  Scintillation  Effects  in  USB  SATOPS  Systems 

SATOPS  communication  links  at  S-band  can  be  subject  to  the  effects  of  ionospheric 
scintillations,  which  are  principally  related  to  the  occurrence  of  F- layer  irregularities. 
Scintillations  cause  both  enhancements  and  fading  of  signal  about  the  medium  level  as  the  signal 
transits  the  disturbed  ionospheric  region.  When  scintillations  occur  that  exceed  the  fade  margin, 
perfonnance  of  the  communications  link  will  be  degraded.  This  degradation  is  most  serious  for 
propagation  paths  that  transit  in  the  auroral  and  equatorial  ionospheres.  The  degree  of 
degradation  will  depend  on  how  far  the  signal  fades  below  the  margin,  the  duration  of  the  fade, 
the  type  of  modulation,  and  the  criteria  for  acceptability.  Table  6  and  Table  7  show  the  signal 
losses  computed  at  15°  and  10°  elevation  angles  for  transmitted  frequencies  of  2.06  GHz  and  1.8 
GHz. 


Table  6.  Signal  Losses  Due  to  Scintillation  for  15°  and  10°  Elevation  Angle  at  2.06  GHz 


Annual  Percent 
of  Time,  p  (%) 

Signal  Loss  Due  to  Ionospheric  Scintillations  (dB) 

Elevation  Angle  0E|V  =  1 5° 

Elevation  Angle  0Eiv  =10° 

0.05 

22.9 

36.5 

0.10 

19.5 

31.6 

0.50 

11.5 

19.5 

Table  7.  Signal  Losses  Due  to  Scintillation  for  15°  and  10°  Elevation  Angle  at  1.8  GHz 


Annual  Percent 
of  Time,  p  (%) 

Signal  Loss  Due  to  Ionospheric  Scintillations  (dB) 

Elevation  Angle  0Eiv  =  1 5° 

Elevation  Angle  0Eiv  =10° 

0.05 

29.1 

45.8 

0.10 

25.1 

39.8 

0.50 

15.2 

25.0 
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3.3.11.5  Degradation  Due  to  RFI  for  USB  SATOPS  Systems 


Sections  3.3.5,  Section  3.3.6,  Section  3.3.7,  Section  3.3.8  and  Section  3.3.9  above  have 
devoted  the  modeling  of  impacts  of  RFI  on  the  second-order  PLL.  The  carrier  SNR  (C/No)  and 
bit  SNR  (Eb/N0)  degradations  can  be  calculated  using  the  analytical  models  presented  in  those 
sections. 


3.3.11.6  Modulation  Losses  for  USB  SATOPS  Waveforms 


From  (1),  the  modulation  losses  for  USB  standard  waveforms,  which  use  PCM/PSK/PM 
commanding  signals  and  PCM/PM-Biphase  commanding  signals  modulated  on  the  carrier,  can 
be  shown  to  have  the  following  fonns,  respectively: 

(1)  PCM/PSK/PM  USB  SATOPS  wavefonn: 

MLc  =  P^  =  Jl(m) 

M  (276) 

M  (277) 

(2)  PCM/PM-Biphase  USB  SATOPS  waveform: 


ML,.  = 


Carrier 


Jo  (m) 


(278) 


ML 


CMD 


'cmd 


2  j\  (in) 


(279) 


where  m  is  the  modulation  index,  Reamer  (or  Pc)  is  the  power  in  the  carrier  component,  and  Rcmd 
is  the  power  in  the  command  signal  component. 

As  mentioned  at  the  beginning  of  this  section,  the  channel  gain  “a”  also  can  be  estimated 
from  the  link  budget  model  described  above.  From  (162),  the  channel  gain  “a”  can  be  estimated 
from  the  received  signal  power,  PR,  at  the  satellite  transponder,  and  PR  can  be  estimated  using 
the  following  relationship: 


PR  =  EIRP(dBm)  +  GfdBi)  -  KTfVc(dBm)  -  Ls(dB)  -  L0(dB) 

where 


(280) 
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L0(dB)  =  LAtm  ( dB )  +  LCaR(dB)  +  ARF,  (dB) 

LAtm  (dB)  =  Loss  due  to  atmospheric  (including  scintillation)  in  dB 
LCaR(dB)  =  Loss  due  to  receiver' s  cable  in  dB 

Am  (dB)  =  Degradation  due  to  RFI  in  dB  ^8 ) 

The  degradation  due  to  RFI,  ARFi(dB),  can  be  calculated  using  the  analytical  models  presented  in 
Section  3.3.7. 


3.3.11.7  Optimum  Transmitted  Signal  Power  for  USB  SATOPS  Waveforms 


The  communications  designer  can  use  the  framework  described  in  Section  3.2  and  the 
algorithms  presented  in  this  section  and  Section  3.3  to  estimate  the  optimum  transmitted  signal 
power,  Pj  opt’,  to  maintain  a  required  USB  SATOPS  Quality-of-Service  (QoS)  in  the  presence  of 
both  friendly  and  unfriendly  RFI  sources.  The  optimum  transmitted  power  can  be  calculated 
using  the  following  relationship: 


(  p  A 

P.,  „  (dB)  =  LM  .(dB)  +  6 


VN0P-eqd 


+  R(  dB  -  bits/s)  +  ML.  (dB)  -  GfdBi) 


T_Opt ' 

+  kTsWc (dBm)  +  Ls (dB)  +  L0 (dB)  -  Gr  (dB)  +  Feed  _ Loss(dB) 


(282) 


where  the  USB  SATOPS  QoS  is  specified  through  the  required  bit  SNR  (Eb/No)reqd  is  through  the 
parameter  Z0(db)  described  in  (281). 


4  RESULTS  AND  DISCUSSION 

4.1  RFI  Detection  and  Prediction  Model  Integration 

Our  approach  was  to  leverage  existing  RFI  analysis  tools  that  were  available  in  public 
domain.  The  team  developed  a  web-based  graphical  user  interface  (GUI)  for  the  RFI  modeling 
and  prediction  tool  with  the  following  capabilities: 

•  Predict  USB  SATOPS  carrier  acquisition  time  consisting  of  frequency  and  phase 
acquisition, 

•  Predict  USB  SATOPS  carrier  tracking  perfonnance,  and 

•  Predict  USB  SATOP  command  BER  performance  in  the  presence  of  AWGN,  continuous 
wave  or  wideband  RFI,  and  imperfect  carrier  tracking  caused  by  AWGN  and  continuous 
wave  or  wideband  RFI. 
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To  develop  the  web-based  GUI,  the  team  employed  the  following  programming  technologies: 

•  Java  scripting  for  front-end  data  handling; 

•  C-programming  to  handle  the  common  gateway  interface  (CGI)  interactions  (including 
data  entry  checking  for  increased  security); 

•  Develop  and  implement  Matlab  codes  for  evaluating: 

o  Carrier  acquisition  time  consisting  of  frequency  and  phase  acquisition, 
o  Carrier  tracking  performance, 

o  Command  BER  perfonnance  in  the  presence  of  AWGN,  continuous  wave  or 
wideband  RFI,  and  imperfect  carrier  tracking  caused  by  AWGN  and  continuous 
wave  or  wideband  RFI. 

The  web-based  GUI  for  the  RFI  modeling  and  prediction  tools  include  features  such  as  double 
precision,  data  verification  routines  to  catch  user  errors  on  input,  and  a  default  feature  case  that 
requires  minimal  user  input.  Figure  5 1  shows  a  screen-shot  of  the  web-based  GUI  for  the  RFI 
modeling  and  prediction  tool. 

4.2  IFT  SATCOM  Tool  Integration 

This  section  describes  the  integration  of  the  RFI  prediction  and  detection  models 
presented  in  the  previous  sections  into  the  IFT  SATCOMM  tool.  The  IFT  tool  has  been 
developed  based  on  the  open  source  Jsattrak  with  Java.  The  integration  consists  of  the  following 
steps: 

•  Design  the  menu  items  to  call  the  display  functions  of  various  analysis  and  simulation 
results  from  various  CW  RFI  and  WB  RFI  models,  including  RFI  prediction,  SATOPS 
command  carrier  acquisition  and  tracking  performances,  and  SATOPS  command  BER 
perfonnance. 

•  Design  a  graphical  user  interface,  which  is  user  friendly  and  functional,  including  a 
display  of  SATOPS  system  setting,  and  corresponding  analytical  and  simulation  results. 

The  programming  languages  required  to  perform  the  integration  task  include  Java  for  front-end 
data  handling  in  Netbeans. 
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Command  Link  Margin  Analysis 
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Figure  51.  RFI  Modeling  and  Prediction  Tool  Integration 

Figure  52  shows  the  overview  of  IFT  SatComm  tool  integrated  with  the  RFI  prediction 
and  detection  models.  One  can  see  that  the  menu  item  “IFT  SatComm”  includes  the  “BER  and 
Tracking  error”  and  “Acquisition  Time  analysis”,  which  are  the  key  components  in  the  RFI 
detection  and  prediction  models  described  earlier. 
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Figure  52.  IFT  SatComm  Tool  Integrated  With  RFI  Prediction  Results 

Figure  53  demonstrates  the  SATOPS  command  BER  performance  in  the  presence  of 
AWGN,  RFI  signal  and  imperfect  carrier  tracking  loop.  On  the  other  hand,  Figure  54  shows  the 
PLL  carrier  tracking  error  performances  under  no  RFI  and  RFI  scenarios.  The  RFI  scenarios 
include  continuous  wave  RFI  and  wideband  RFI.  One  can  see  that  in  the  display  window,  it  has  a 
test  area  displaying  the  information  of  system  setting.  One  will  be  able  to  switch  between  BER 
results  and  the  carrier  PLL  tracking  error  by  the  buttons. 
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Figure  53.  BER  Performance  With  CW  RFI  or  WB  RFI  and  Imperfect  Carrier  Tracking 
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Figure  54.  Tracking  Error  Results  Under  No  RFI  and  CW  RFI  or  WB  RFI  Scenario 
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Figure  55  and  Figure  56  show  the  display  for  the  acquisition  time  perfonnances  under 
CW  RFI  and  WB  RFI,  respectively.  The  display  window  includes  a  test  area  displaying  system 
setting  information  and  the  switch  buttons.  For  the  acquisition  time  performance  analysis,  there 
are  four  figures  in  the  windows:  (1)  relationship  between  SNR  and  loop  SNR,  (2)  frequency 
acquisition  time  as  a  function  of  loop  SNR,  (3)  phase  acquisition  time  as  a  function  of  loop  SNR, 
and  (4)  the  total  acquisition  time  of  the  carrier  tracking  loop  as  a  function  of  loop  SNR. 


Figure  55.  Acquisition  Time  Result  Under  CW  RFI 
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Figure  56.  Acquisition  Time  Result  Under  WB  RFI 


4.3  Verification  and  Validation  of  RFI  Detection  and  Prediction  Tools 

The  analytical  models  for  RFI  detection  and  prediction,  which  were  derived  and 
discussed  in  Section  3.3,  were  implemented  in  MATLAB  for  verification  and  validation  (V&V) 
purpose.  These  models  were  tested  and  verified  for  various  operating  scenarios  to  ensure  the 
accuracy  of  the  performance  prediction  of  a  typical  USB  SATOPS  command  system.  Figure  57- 
Figure  60  depict  the  impacts  of  CW  and  wideband  RFI  signals  on  the  second  order  PLL 
acquisition  perfonnance,  respectively.  In  particular,  they  compare  the  carrier  acquisition  times 
without  RFI  and  in  the  presence  of  RFI  as  functions  of  loop  SNR  for  the  following  cases: 

•  Signal  carrier  frequency  =fc  =  2075  MHz, 

•  m  =  command  modulation  index  =1.1  rad, 

•  The  threshold  carrier  jitter  a  =n/2  rad, 

•  Command  bit  rate  =  R\,  =  2  kbps,  and 

•  Interference-to-signal  power  ratio  =  ISR  =  -10  dB  and  -40  dB. 
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Figure  57  shows  the  assessment  of  PLL  acquisition  performance  in  the  presence  of  CW  RFI  at 
ISR  =  -10  dB.  The  plots  of  loop  SNR  versus  received  SNR  without  RFI  and  with  CW  RFI  are 
shown  on  the  top  left  of  Figure  57.  The  plots  of  carrier  frequency  acquisition  times  as  functions 
of  loop  SNR  without  RFI  and  with  CW  RFI  are  shown  on  the  top  right  of  Figure  57.  The  plots  of 
carrier  phase  acquisition  times  as  functions  of  loop  SNR  without  RFI  and  with  CW  RFI  are 
shown  on  the  bottom  left  of  Figure  57.  As  expected,  the  acquisition  time  increases  due  to  the 
presence  of  CW  RFI. 


t  =2075  MHz.  m=l.l  rad.  a  =  i/2  fad.  R  =  2000 dps.  ISR  = -ID  dB 

G  D 


Figure  57.  Plots  of  Total  Carrier  Acquisition  Times  Versus  Loop  SNR  Without  RFI  and 

With  CW  RFI 


Approved  for  public  release;  distribution  is  unlimited. 


107 


Figure  58  provides  results  for  the  assessment  of  the  PLL  acquisition  performance  in  the 
presence  of  WB  RFI  at  ISR  =  -10  dB.  The  plots  of  loop  SNR  versus  SNR  without  RFI  and  with 
WB  RFI  are  shown  on  the  top  left  of  Figure  58.  The  plots  of  carrier  frequency  acquisition  times 
as  functions  of  loop  SNR  without  RFI  and  with  WB  RFI  are  shown  on  the  top  right  of  Figure  58 
The  plots  of  carrier  phase  acquisition  times  as  functions  of  loop  SNR  without  RFI  and  with  WB 
RFI  are  shown  on  the  bottom  left  of  Figure  58.  The  results  show  that  the  acquisition  time 
increases  due  to  the  presence  of  WB  RFI. 


1  =2075  MHz.  m=1.i  rad,  *  =  ,i/2rad,  R  =2000  Dps.  ISR  =  -l0dB.  R  =1  Dps 

c  b  I 


LoopSNR)V  ;|  (dB) 


Figure  58.  Plots  of  Total  Carrier  Acquisition  Times  as  Functions  of  Loop  SNR  Without 

RFI  and  With  WB  RFI  at  ISR  =  -10  dB 
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Figure  59  provides  the  assessment  of  PLL  acquisition  performance  in  the  presence  of  CW 
RFI  at  ISR  =  -40  dB.  The  plots  of  loop  SNR  versus  SNR  without  RFI  and  with  CW  RFI  are 
shown  on  the  top  left  of  Figure  59.  The  plots  of  carrier  frequency  acquisition  times  as  functions 
of  loop  SNR  without  RFI  and  with  CW  RFI  are  shown  on  the  top  right  of  Figure  59.  The  plots  of 
carrier  phase  acquisition  times  as  functions  of  loop  SNR  without  RFI  and  with  CW  RFI  are 
shown  on  the  bottom  left  of  Figure  59.  The  results  show  that  the  total  acquisition  time 
degradation  is  negligible  when  the  CW  RFI  interference  -to-signal  power  ratio,  ISR,  is  at  -40  dB. 


f  =2075  MHz.  m=i.t  rad.  a  -  i/2rad.R„  =2000 Dps.  ISR  =  -40dB 

a  D 


Figure  59.  Plots  of  Total  Carrier  Acquisition  Times  as  Functions  of  Loop  SNR  Without 

RFI  and  With  CW  RFI 
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Figure  60  presents  the  assessment  of  PLL  acquisition  performance  in  the  presence  of  WB 
RFI  at  ISR  =  -40  dB.  The  plots  of  loop  SNR  versus  SNR  without  RFI  and  with  WB  RFI  are 
shown  on  the  top  left  of  Figure  60.  The  plots  of  carrier  frequency  acquisition  times  as  functions 
of  loop  SNR  without  RFI  and  with  WB  RFI  are  shown  on  the  top  right  of  Figure  60.  The  plots  of 
carrier  phase  acquisition  times  as  functions  of  loop  SNR  without  RFI  and  with  WB  RFI  are 
shown  on  the  bottom  left  of  Figure  60.  Similarly,  for  WB  RFI,  the  numerical  results  show  that 
the  total  acquisition  time  degradation  is  negligible  when  the  ISR  is  at  -40  dB. 


t.  =2075  MHz.  mal.l  rad,  *  =  ,i/2rad,  R;  =2000 Dps,  ISR  =  -40dB.  R,  =  1  Dps 


Figure  60.  Plots  of  Total  Carrier  Acquisition  Times  as  Functions  of  Loop  SNR  Without 

RFI  and  With  WB  RFI  at  ISR  =  -40  dB 


Approved  for  public  release;  distribution  is  unlimited. 


110 


Figure  61  shows  plots  of  PLL  tracking  jitter  in  the  absence  of  RFI  and  with  both  CW  and 
WB  RFI  signals  as  functions  of  loop  SNR  with  A/rh  =  5  Hz,  in  =1.1  rads,  amargin  =  2  deg,  Rb.=  2 
kbps,  ISR  =  -KMB,  R\.=  10  bps,  B\~  10  Hz.  The  plots  show  that  the  tracking  performance  of  the 
PLL  in  the  presence  of  CW  RFI  is  worse  than  WB  RFI  under  the  specified  operating  conditions. 


A  t-_.  =  5  Hz,  m  =  i.i  rad,"  =2ctog.  ft  =2000 Dps. I5R--10 dB.R=  loops.  B.  =10 hz 

RFI  '  margn  3  b  ^  I  L 


Figure  61.  Plots  of  PLL  Tracking  Jitter  in  the  Absence  of  RFI  and  Presence  of  CW  and  WB 
RFI  Signals  as  Functions  of  Loop  SNR  With  A/rfi  =  5  Hz 

Figure  62  shows  plots  of  PLL  tracking  jitter  in  the  absence  of  RFI  and  with  both  CW  and 
WB  RFI  signals  as  functions  of  loop  SNR  with  A/rh  =10  Hz,  in  =  1.1  rads,  amargin  =  2  degs,  Rb.= 
2  kbps,  ISR  =  -10  dB,  Rh=  10  bps,  B\  =  10  Hz.  The  plots  show  that  the  tracking  performance  of 
the  PLL  in  the  presence  of  CW  RFI  is  still  worse  than  WB  RFI  under  the  same  operating 
conditions.  For  this  case,  the  plots  show  that  the  carrier  tracking  performance  of  the  WB  RFI  is 
better  than  the  previous  case  for  A/rFi  =  5  Hz.  This  is  expected,  since  the  RFI  is  moving  away 
from  the  center  carrier  frequency. 
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Figure  62.  Plots  of  PLL  Tracking  Jitter  in  the  Absence  of  RFI  and  Presence  of  CW  and  WB 
RFI  Signals  as  Functions  of  Loop  SNR  With  A/'rm  =  10  Hz 

Figure  63  and  Figure  64  depict  the  BER  perfonnance  due  to  AWGN  and  RFI  signal  as 
functions  of  loop  SNR  taking  into  account  carrier  synchronization  loop.  Figure  63  shows  plots  of 
BER  performance  in  the  absence  of  RFI  and  with  both  CW  and  WB  RFI  signals  as  functions  of 
loop  SNR  with  A/R|.|  =  5  Hz,  m  =  1.1  rads,  amargin  =  2  degs,  Rb  =  2  kbps,  ISR  =  -10  dB,  R\  =  1 0 
bps,  Bi  =  10  Hz.  The  plots  show  that  the  BER  performance  of  the  PLL  in  the  presence  of  CW 
RFI  is  worse  than  WB  RFI  under  the  same  operating  conditions. 

Figure  64  shows  plots  of  BER  performance  in  the  absence  of  RFI  and  with  both  CW  and 
WB  RFI  signals  as  functions  of  loop  SNR  with  the  same  operating  conditions  shown  in  Figure 
63,  except  with  A/rfi  =10  Hz.  The  plots  show  that  the  BER  performance  of  the  PLL  in  the 
presence  of  CW  RFI  is  still  worse  than  WB  RFI  under  the  same  operating  conditions.  But  when 
compared  to  this  case  with  Figure  63,  the  plots  show  that  the  BER  performance  of  the  WB  RFI  is 
better  than  the  previous  case  for  A/Rfi  =  5  Hz.  As  pointed  out  earlier,  this  is  as  expected,  since 
the  RFI  is  moving  away  from  the  center  carrier  frequency. 
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A  tRR  =  5  Kz.  m  =  M  rad,  anm^n  -2  dog.  Rfc  =2000  Dps.  !5R=-10dB.  R(  =  10  Dps.  Bl  =  10  HZ 


Figure  63.  Plots  of  BER  due  to  RFI  for  Both  CW  RFI  and  WB  RFI  Signals  as  Functions  of 

Foop  SNR  With/RFi  =  5  Hz 


A  l__.  =  10  Hz.  m  =  i.i  rad.  a  =  2  deg.  R.  =2000  Dps,  I5R  =  -io  dB.  R,  =  to  Dps.  B  =iohz 
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Figure  64.  Plots  of  BER  due  to  RFI  for  Both  CW  RFI  and  WB  RFI  Signals  as  Functions  of 

Foop  SNR  With/RFi  =  10  Hz 
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5  CONCLUSIONS 


This  report  presented  the  work  done  by  the  Catholic  University  of  America  (CUA) 
research  team,  which  included  North  Carolina  State  University  (NCSU)  and  Intelligent  Fusion 
Technology  (IFT)  on  radio  frequency  interference  (RFI)  modeling  and  prediction.  The  team 
accomplished  research  objectives  and  completed  the  six  tasks  described  in  our  proposal. 

The  survey  results  showed  that  the  proposed  RFI  detection  approach  using  a  SATOPS 
receiver’s  carrier  synchronizer  presented  in  this  report  is  original.  The  CUA  team  conducted 
research  on  more  than  6  RFI  tools.  We  leveraged  existing  tools  for  incorporating  the  RFI 
detection  and  prediction.  We  developed  a  framework  and  models  to  assess  the  impacts  of  RFI  on 
SATOPS  USB  systems. 

An  innovative  framework  was  developed  by  the  CUA  team  with  an  intention  to  provide  a 
unified  approach  to  address  the  RFI  challenges  for  USB  SATOPS.  We  described  the  framework 
for  the  RFI  tool  development,  RFI  detection  and  assessment,  and  RFI  prediction.  The  RFI  tool 
development  framework  provided  detailed  description  of  the  tools  and  models  required  for 
detecting,  predicting  and  assessing  RFI  signals  on  USB  SATOPS  command  systems.  The  RFI 
detection  framework  described  how  to  incorporate  the  SATOPS  receiver’s  carrier  synchronizer 
in  the  detection  of  CW  and  WB  RFI  events.  The  RFI  assessment  framework  presented  a  unified 
approach  for  evaluating  the  impacts  of  CW  and  WB  RFI  sources  on  USB  SATOPS  command 
systems.  The  RFI  prediction  framework  addressed  a  unified  approach  for  estimating  and 
predicting  the  characteristics  of  CW  and  WB  RFI  sources.  It  is  worth  mentioning  that  the 
framework  is  applicable  to  both  satellite  communications  and  satellite  operations. 

The  team  developed  analytical  and  simulation  models  to  evaluate  the  carrier  acquisition 
and  tracking  performances  of  practical  USB  SATOPS  command  systems  in  the  presence  of  CW 
and  WB  RFI  signals.  The  impacts  of  the  RFI  signals  on  the  carrier  synchronizer  were  used  in  the 
detection  of  the  RFI  events.  We  described  an  approach  to  accurately  predict  the  RFI  interfering 
time  duration  using  the  carrier  synchronizer  and  assessed  the  impacts  of  the  carrier  perfonnance 
degradation  of  the  synchronizer  on  the  command  BER  performance.  Advanced  signal  processing 
algorithms  to  estimate,  predict  and  characterize  the  CW  and  WB  RFI  signals  were  described  in 
detail  along  with  analytical  models  to  describe  the  estimators’  perfonnance.  The  characteristics 
of  the  RFI  signals  estimated  and  predicted  by  our  team’s  proposed  algorithms  included  RFI 
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power,  RFI  carrier  frequency,  RFI  phase,  and  RFI  bandwidth.  In  summary,  the  key  RFI 
analytical  and  simulation  models  that  were  developed  under  this  task  included: 

•  Models  to  predict  USB  SATOPS  carrier  acquisition  time  in  the  presence  of  CW/WB  RFI. 
The  acquisition  time  consists  of  frequency  and  phase  acquisition. 

•  Models  to  predict  USB  SATOPS  carrier  tracking  perfonnance. 

•  Model  to  predict  USB  SATOP  command  BER  perfonnance  in  the  presence  of  AWGN, 
CW/WB  RFI,  and  imperfect  earner  tracking  caused  by  AWGN  and  CW/WB  RFI. 

•  Models  to  predict  RFI  duration  by  monitoring  canier  synchronizer. 

•  Model  to  detect  RFI  events  by  monitoring  carrier  synchronizer. 

•  Models  to  estimate,  predict,  and  characterize  CW  RFI  signals. 

•  Models  to  estimate,  predict,  and  characterize  WB  RFI  signals. 

•  Models  to  estimate  the  optimum  transmitted  USB  SATOPS  command  signal  power  to 
maintain  a  required  USB  SATOPS  quality-of-service  (QoS)  in  the  presence  of  CW  and 
WB  RFI  signals. 

This  report  described  the  team’s  approach  to  integrate  the  analytical  and  simulation 
models  developed  in  Task  #3  into  existing  RFI  tool  developed  by  The  Aerospace  Corporation. 
The  newly  integrated  RFI  tools  have  new  capabilities  including:  (i)  predict  USB  SATOPS  carrier 
acquisition  time  consisting  of  frequency  and  phase  acquisition,  (ii)  Predict  USB  SATOPS  carrier 
tracking  performance,  (iii)  predict  USB  SATOP  command  BER  perfonnance  in  the  presence  of 
AWGN,  CW/WB  RFI,  and  imperfect  carrier  tracking  caused  by  AWGN  and  CW/WB  RFI,  (iv) 
predict  RFI  duration  by  monitoring  carrier  synchronizer,  and  (v)  detect  RFI  events  by  monitoring 
canier  synchronizer.  Currently,  only  capabilities  (i),  (ii)  and  (iii)  were  verified  and  validated. 

We  presented  the  key  software  integration  results  obtained  by  the  team  for  Task  #5.  We 
discussed  how  the  integrated  RFI  Tool  developed  in  Task  #4  were  incorporated  into  existing 
Intelligent  Fusion  Technology  SATCOM  tools  to  display  the  RFI  detection  and  prediction 
results.  Screen-shots  of  how  the  tools  display  the  RFI  results  for  USB  SATOPS  command 
systems  were  provided.  Cunently,  the  team  only  incorporated  three  capabilities  into  the 
Intelligent  Fusion  Technology  SATCOM  tool,  including  prediction  of  USB  SATOPS  carrier 
acquisition  time  consisting  of  frequency  and  phase  acquisition,  prediction  of  USB  SATOPS 
carrier  tracking  perfonnance,  and  prediction  USB  SATOP  command  BER  perfonnance  in  the 
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presence  of  AWGN,  CW/WB  RFI,  and  imperfect  carrier  tracking  caused  by  AWGN  and 
CW/WB  RFI. 

We  provided  a  summary  of  the  team’s  approach  for  Task  #6  on  the  verification  and 
validation  of  the  analytical  and  simulation  models  developed  in  Task  #3.  The  following  models 
were  verified  and  validated  by  the  team:  (i)  models  to  predict  USB  SATOPS  carrier  acquisition 
time  consisting  of  frequency  and  phase  acquisition,  (ii)  models  to  predict  USB  SATOPS  carrier 
tracking  performance,  and  (iii)  models  to  predict  USB  SATOPS  command  BER  performance  in 
the  presence  of  AWGN,  CW/WB  RFI,  and  imperfect  carrier  tracking  caused  by  AWGN  and 
CW/WB  RFI. 

In  conclusion,  the  team  successfully  developed  a  unified  framework  to  address  the  RFI 
modeling  and  prediction  challenges  by  incorporating  effects  of  RFI  on  the  satellite  receiver’s 
carrier  synchronizers,  and  sophisticated  analytical  and  simulation  models  to  thoroughly  evaluate 
the  impacts  of  CW  and  WB  RFI  signals  on  the  USB  SATOPS  command  system.  The  team  also 
developed  advanced  signal  processing  algorithms  to  estimate/predict/characterize  CW  and  WB 
RFI  signals.  Using  RFI  modeling  and  prediction  models  developed  by  this  research  effort,  a 
communications  designer  can: 

•  Evaluate  effectiveness  of  the  existing  USB  SATOPS  command  waveforms  employed  by 

civil,  commercial  and  military  SATOPS  ground  stations; 

•  Predict  the  impacts  of  RFI  on  USB  SATOPS  command  systems; 

•  Characterize  both  friendly  and  unfriendly  RFI  sources;  and 

•  Estimate  the  optimum  transmitted  signal  power  to  maintain  a  required  USB  SATOPS 

Quality-of-Service  (QoS)  in  the  presence  of  both  friendly  and  unfriendly  RFI  signals. 

6  RECOMMENDATIONS 

Currently,  we  have  not  yet  fully  incorporated  all  of  the  analytical  and  simulation  models 
that  have  been  developed  into  the  existing  RFI  tool  and  IFT  SATCOM  tool.  We  have  only 
incorporated  the  models  to  estimate  and  predict: 

(i)  USB  SATOPS  carrier  acquisition  time  in  the  presence  of  CW/WB  RFI, 

(ii)  USB  SATOPS  carrier  tracking  performance,  and 
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(iii)  USB  SATOPS  command  BER  perfonnance  in  the  presence  of  AWGN,  CW/WB 
RFI,  and  imperfect  carrier  tracking  caused  by  AWGN  and  CW/WB  RFI. 

To  effectively  address  the  RFI  challenges  on  the  USB  SATOPS  systems,  we  recommend 
incorporating  the  following  models  into  the  existing  RFI  tool  and  IFT  SATCOM  tool: 

•  Models  to  predict  RFI  duration  by  monitoring  carrier  synchronizer, 

•  Models  to  detect  RFI  events  by  monitoring  carrier  synchronizer, 

•  Models  to  estimate,  predict  and  characterize  CW  RFI  signal, 

•  Models  to  estimate,  predict  and  characterize  WB  RFI  signal,  and 

•  Models  to  estimate  the  optimum  transmitted  USB  SATOPS  command  signal  power  to 
maintain  a  required  USB  SATOPS  Quality-of-Service  (QoS)  in  the  presence  of  CW  and 
WB  RFI  signals. 

The  focus  of  this  research  project  was  on  the  USB  SATOPS  carrier  synchronizer.  It  is  anticipated 
that  the  influences  of  the  RFI  signals  on  the  command  subcarrier  synchronizer  and  symbol 
timing  synchronizer  are  non-trivial,  and  they  should  also  be  incorporated  into  the  detection  of  the 
RFI  events  and  prediction  of  the  RFI  signals.  We  recommend  investigating  the  following 
research  topics: 

•  USB  SATOPS  subcarrier  and  timing  acquisition  time  in  the  presence  of  CW/WB  RFI, 

•  USB  SATOPS  subcarrier  and  timing  tracking  performance, 

•  USB  SATOPS  command  BER  perfonnance  in  the  presence  of  imperfect 
subcarrier/timing  tracking  caused  by  AWGN  and  CW/WB  RFI, 

•  Prediction  of  RFI  duration  by  monitoring  subcarrier/timing  synchronizers,  and 

•  Detection  of  RFI  events  by  monitoring  subcarrier/timing  synchronizers. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


3G 

Third  Generation  Cellular 

4G 

Fourth  Generation  Cellular 

8PSK 

8-Phase  Shift  Keying 

16QAM 

16-Quadrature  Amplitude  Modulation 

a 

Channel  Gain 

ACI 

Adjacent  Channel  Interference 

AEHF 

Advanced  Extremely  High  Frequency 

AFSCN 

Air  Force  Satellite  Control  Network 

AGI 

Analytical  Graphics,  Incorporated 

AWGN 

Additive  White  Gaussian  Noise 

-^Acq-DL 

Acquisition  Mode  Foop  Bandwidth 

Foop  Bandwidth 

BER 

Bit  Error  Rate 

Bl 

Foop  Bandwidth 

B  LAcq 

Acquisition  Foop  Bandwidth 

Blp 

Foop  Bandwidth  Perfect  Loop  Filter 

BPF 

Band  Pass  Filter 

bps 

Bits  per  Second 

B  rcx 

SATOPS  Receiver  Bandwidth 

c 

Speed  of  Light 

CCI 

Co-Channel  Interference 

CCSDS 

Consultative  Committee  for  Space  Data  System 

CDMA 

Code  Division  Multiple  Access 

CDU 

Command  Detector  Unit 

CGI 

Common  Gateway  Interface 

C/I 

Carrier  to  Interference  Ratio 

CMM 

Carrier  Modulation  Mode 
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C/No 

Carrier  to  Noise  Ratio 

cos 

Cosine 

CSNR 

Carrier  Signal  Power-to-Noise  Power  Spectral  Density  Ratio 

CSNRwb 

Wideband  Carrier  Signal-to-Noise  Power  Spectral  Density  Ratio 

CUA 

Catholic  University  of  America 

CW 

Continuous  Wave 

dB 

Decibel 

dBi 

Decibel  Relative  to  Isotropic  Gain 

dBm 

Decibel  relative  to  a  rnilli  Watt 

DC 

Direct  Current,  0  Hz 

DFT 

Discrete  Fourier  Transform 

DTTL 

Digital  Transition  Tracking  Loop 

Eb 

Energy  per  Bit 

EIRP 

Effective  Isotropically  Radiated  Power 

Eq.,  Eqn. 

Equation 

fc 

Carrier  Frequency 

FCC 

Federal  Communications  Commission 

h 

Doppler  Frequency  Shift 

FOV 

Field  of  View 

/rfi 

RFI  Frequency 

Fs 

Sampling  Frequency 

1  samp 

Sampling  Frequency 

/sc 

Sub-carrier  Frequency 

fj n 

Frequency  Uncertainty 

GBS 

Global  Broadcast  Service 

GCN 

Ground  Communications  Network 

GMSK 

Gaussian  Minimum  Shift  Keying 

Gr 

Antenna  Gain 
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GTS 

Ground  Tracking  Station 

GUI 

Graphical  User  Interface 

Hz 

Hertz 

ICSR 

Interference  Signal  Power-to-Carrier  Signal  Power  Ratio 

IF 

Intennediate  Frequency 

IFT 

Intelligent  Fusion  Technology 

INR 

Interference  Power  to  Noise  Ratio 

IPC 

Interference  Protection  Criteria 

ISR 

Interference  to  Signal  Ratio 

ITU 

International  Telecommunication  Union 

Avu(t) 

Wideband  Interference  Signal 

iWPA 

Iterative  Weighted  Phase  Averaging 

JtoJl 

Bessel  Function,  First  Kind 

JPL 

Jet  Propulsion  Laboratory 

kbps 

Kilo  Bits  per  Second 

kHz 

Kilo  Hertz 

LMcarrier 

Link  Margin  Model  for  Carrier 

LNA 

Low  Noise  Amplifier 

LPF 

Low  Pass  Filter 

LSNR 

Loop  Signal  to  Noise  Ratio 

m 

Modulation  index 

MHz 

Mega  Hertz 

ML 

Maximum  Likelihood 

n 

Bit  Rate  Ratio 

NASA 

National  Aeronautics  and  Space  Administration 

NCC 

Network  Control  Center 

NCSU 

North  Carolina  State  University 

No 

Noise  Power 
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Nob 

Observation  Samples 

NRZ 

Non  Return  to  Zero 

NTIA 

National  Telecommunications  and  Infonnation  Administration 

P  P 

1  cn1  earner 

Power  in  Carrier 

PCM 

Pulse  Code  Modulation 

P  (Ml) 

Power  in  Command  Signal 

PCS 

Pulse  Coded  Signals 

PFD 

Power  Flux  Density 

Pi 

Interference  Power 

PLL 

Phased  Locked  Loop 

PM 

Phase  Modulation 

PRF 

Pulse  Repetition  Frequency 

PSD 

Power  Spectral  Density 

PSK 

Phase  Shift  Keying 

Pt 

Transmit  Power 

Pt 

Total  Transmitted  Power 

Q(a,b) 

Marcum  Q-function 

QoS 

Quality  of  Service 

QPSK 

Quadrature  Phase  Shifted  Keying 

Rad,  rad 

radians 

Rb 

Bit  rate  in  bits  per  second 

RF 

Radio  Frequency 

RFI 

Radio  Frequency  Interference 

Ri(t) 

Received  Signal 

R\\m(t) 

Received  Signal,  Baseband 

^ibbn(0 

Received  Signal,  Baseband 

^iwb(0 

Received  Wideband  Interference  Signal 

^iwbb(I) 

Received  Baseband  Wideband  Interference  Signal 
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RL 

Radio  Loss 

ROC 

Receiver  Operating  Characteristic 

Rs 

Bit  rate  in  bits  per  second 

s 

Laplace  Variable, yco 

SATCOM 

Satellite  Communications 

SATOPS 

Satellite  Operations 

see 

Satellite  Control  Center 

SCN 

Satellite  Control  Network 

SGLS 

Space-to-Ground  Link  Subsystem 

sin 

Sine 

SNR 

Signal  to  Noise  Ratio 

SPW 

Signal  Processing  Workstation 

STK 

Satellite  Tool  Kit 

Tacq 

Carrier  Frequency  Acquisition  Time 

Th 

Bit  Duration  in  Seconds 

TDRSS 

Tracking  and  Data  Relay  Satellite  System 

TF 

Time  Factor 

Fob 

Observation  Time,  in  Seconds 

FpAcq 

Carrier  Phase  Acquisition  Time 

Ts 

Bit  Duration  in  Seconds 

Ts 

Sample  Period 

TT&C 

Tracking,  Telemetry,  and  Command 

USB 

Unified  S  Band 

VCO 

Voltage  Controlled  Oscillator 

WB 

Wideband 

W-CDMA 

Wideband  Code  Division  Multiple  Access 

WiFi 

Local  Area  Wireless  Computer  Networking  Technology 

WiMAX 

Worldwide  Interoperability  for  Microwave  Access 
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A/c 

RFI  frequency  offset  from  the  desired  carrier  frequency 

s(t) 

Error  Signal 

n 

Antenna  Efficiency 

K 

Boltzmann’s  Constant 
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